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FOREWORD 


This report was submitted by Ultrasystems, Inc., Environmental and Applied 
Sciences Division, 2400 Mlchelson Drive, Irvine, California 92664, under Contract 
No. F04611-73-C-0G38, Job Order No. 305909LZ with the Air Force Rocket Propul- 
sion Laboratory, Edwards, CA 93523. 

This report consists of three volumes. Volume I describes a computer 
program for the prediction of Solid Propellant Rocket Motor Performance. Tne com- 
puter program described herein will be referred to as the SPP program, and describes 
the engineering analysis which was used in developing this computer program and 
the results obtained to date. 


Volume II of this report is a programming document of the computer program 
which was developed under this contract. It Includes a subroutine-by-subroutine 
description of all of the elements of the SPP program. 

Volume III of ".his report is a Program User's Manual which describes the 
input necessary to execute the SPP computer program and the information required 
to interpret the output. A sample case is also included. 

This report has been reviewed by the Information Offlce/DOZ and is re- 
leasable to the National Technical Information Service (NTIS). At NTIS it will be 
available to the general public, including foreign nations. 


This report is unclassified and suitable for public release. 
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computed as a series of independent efficiencies. The program currently treats 
the following losses: two-dimensional/two-phase (coupled), nozzle erosion, 
kinetics, boundary layer, combustion efficiency, submergence. The program 
predicts average delivered performance, as well as mass flow, pressure, thrust, 
impulse, and specific impulse as functions of time and trajectory. 

In order to assess the validity of the SPP program, calculated results were com- 
pared to firing data for four different types of motors. 'Vhile conclusive state- 
ments regarding the accuracy and range of validity of the SPP program cannot be 
made until additional verification efforts are conducted, the results of these four 
test cases were encouraging. These calculations also served to demonstrate 
i the desirability of eliminating some of the present limitations of the program. 
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PREFACE 


This is Volume in of a three-part manual which describes a computer pro- 
gram for the prediction of Solid Propellant Rocket Motor Performance. The com- 
puter program described herein will be referred to as the SPP program. 

Volume I of this report describes the engineering analysis which was used 
in developing thU computer program. 

Volume II of this report is a programming document of the computer program 
which was developed under this contract. It includes a subroutine-by-subroutine 
description of all of the elements of the SPP program. 

Volume HI of this report is a Program User's Manual which describes the 
input necessary to execute the SPP computer program and the information required 
to interpret the output. A sample case is also included. 
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1 . METHOD OF APPROACH 


The goal of the present effort was to develop a computerized analysis tech- 
nique that could predict delivered specific impulse to within ±2%, and delivered 
thrust coefficient and total impulse to within ±S%. The method of approach used 
in the development of the Solid Rocket Performance Program (SPP) stressed modu- 
larity and flexibility. The program was designed so that the various elements 
which comprise the total performance prediction methodology could be utilized as 
an integrated set, individually, or in any user determined combination. In most 
cases both analytical and empirical methods have been incorporated into the pro- 
gram to account for the phenomena which causes motor performance to deviate 
from a specified “ideal.” The user, depending on his needs, can select either method. The 
option is phenomena oriented, so some effects may be treated empirically (simple, 
short computation time) while others can be evaluated analytically (more calcula- 
tion time, greater accuracy and generality). 


The program has been structured in a modular fashion to allow future im- 
provements and mod-flcatlons to be easily incorporated. The program modules 
(with the exception ox the master control module) arc pre-existing computer 
programs which have been adapted for the present purpose. The following pre- 
existing programs were used as the basis for the SPP program: 

One Dimensional Equilibrium Program (ODE)^ ^ 

One Dimensional Kinetics Program (ODKr ~ 

Two Dimensional Two Phase Program (TD2P)^ ^ 

Turbulent Boundary Layer Program (TBL)^ 4 ^ 

Three Dimensional Gram Design Program (64101 5 ^ 

Motor Ballistics Program (637 ® ^ 


These programs have been combined with a control inrvjuir* which permits them to 
be run together automatically. The extensive use of internal communication be- 
tween modules has Just about eliminated the need for redundant inputs by the 
user. 


The next two subsections will define the “losses* which have been con- 
sidered and the manner in which delivered performance is calculated. The balance 
of the section briefly iescribcs ihe various modules, what they do, how they do 
it, and their inter-relctionships. 
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1 . 1 Termlnoloc v and Definition of Losses 

The basic method of approach used in the development of the SPP computer 
program is to calculate the maximum or ideal ‘theoretical* performance of a rocket 
motor and to subtract from that maximum performance the individual losses or de- 
viations from the “theoretical ** performance which ere known to occur in a real 
rocket motor. This method of approach is neither new to the field of rocket per- 
formance prediction or engineering analysis in general. It is, however, founded 
on our understanding of the basic physical processes involved in a real rocket 
motor and the interrelation of these processes with each other. 

The basic assumption involved in this type of approach is that certain "losses” 
can be treated as essentially Independent from the other losses which occur in 
the rocket motor. As noted in Section 1.2, all of the losses should not be treated 
independently. Ther<= are some Interactions strong enough to warrant considera- 
tion. With these exceptions in mind, the concept of treating the losses indepen- 
dently appears to be valid, considering the success of this approach and others 
of a similar nature. 

The term “loss" as used here is really a misnomer, since it is used to 
describe a known physical mechanism or phenomena which causes the performance 
of a “real" solid rocket motor to be less than the so-called “theoretical" perfor- 
mance. In the same sense the terminology "theoietical, “ or '‘ideal, “ if. also a 
misnomer, because the whole methodology Incorporated herein is based on our 
theoretical understanding of the processes which are pre^r:* in real rocket motors. 
However, this terminology has gained wide spread acceptance in the engineering 
community snd will be used In this report. 

In order to further define the terminology used here the following defini- 
tions will be used: 


• Performance 


unless otherwise specified, performance 
will bo In terms of specific impulse, I , 
delivered to a vacuum. 5 " 

•Theoretical I 

sp 


tho maximum possible I gp which can be de- 
livered at the initial area ratio by a propel- 
lant at a given chamber pressure and total 
enthalpy. Referred to as 

• lose 


the decrement from the theoretical I #p which 
can be attributed to a physical phenomena 
net included in the calculation of I 5p . 



• Two phase tlow loss 

•Two dimensional or 
divergence loss 

•Finite Rate Kinetics 
loss 

•Boundary Layer Loss - 
•Combustion Efficiency - 


•Submergence Loss 
•Erosion Loss 


•Delivered Performance - 


the decrement in performance due to finite 
velocity and temperature differences be- 
tween the gas and condensed phase. 

the decrement in performance which is due 
to the momentum of the rocket exhaust not 
being totally aligned with the axis of the 
motor. 

the decrement in performance which is due 
to incomplete transfer of latent heat to sen- 
sible heat caused by the finite time required 
for gas phase chemical reactions to occur. 

the decrement in performance due to viscous 
forces adjacent to the nozzle wall and heat 
transfer to the nozzle wall. * 

as used here, the term combustion efficiency 
will refer only to the degradation of I sp due 
to a departure of chamber total temperature 
from the theoretically calculated total temper- 
ature. 

the decrement in performance due to propel- 
lant grain extending pass the nozzle inlet. 

the decrement in performance due to the ero- 
sion of the nozzle throat. This loss encom- 
passes both area ratio effects (considered 
only in the TD2P module) and throat roughness 
effects (not cons idered in the present analysis) 

the calculated performance of the rocket 
motor which includes all of the losses con- 
sidered herein. 


With our terminology now defined, we can now describe the rational and 
methods used by the SPP computer program to calculate the delivered performance 
of a solid rocket motor. 


1.1.1 Calculation of Delivered Performance 

Within the framework already discussed (Le., separation of individual 

losses), there are trute basic methods whibh can be used to calculate delivered 

I . The first method is to calculate discrete decrements in performance, aI , 
sp ®P 

ami to subtract these losses from the theoretical specific impulse, I • 

a.... 


,# Pd" I,p 


th 


-E Al 


rp 


lots 


(1-1) 
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The second method Is to calculate a series of efficiencies, n’s, to be applied 
multiplicative!/, i.e. 


Isp D* I,p th ^ 10,8 * 


( 1 - 2 ) 


Both of the above methods are the same through first order, thus, if the losses 
are small they yield the same result. 


The third method was selected by the JANNAF Committee on Performance 
Standardization for liquid rocket engines and is a combination of the first two 
methods. That is, all of the losses except the boundary layer loss are treated as 
efflciences while the boundary layer loss is treated as a decrement. Hence, 


Isp D* ^Pth 1 ^ 10 ” * * Isp BL 


(1-3) 


This latter method of apr. roach was selected for incorporation in the SPP code for 
two reasons. The first reason being that there was no compelling reason to de- 
viate from what has been established as a "standardized" method in the Industry, 
and secondly, that any other combination of the above methods yields substantially 
the samo end result. 

The main p-colem associated with the separation technique is, of course, 
the selection of a tase or reference value of performance to compare against the 
model which incorporates the loss. Ideally, this reference calculation should be 
made by the same computer code which is selected to sake the loss calculation 
since this would minimize deviations due to numerical techniques, truncation and 
roundoff errors. However, as is generally the case, su«“’i an approach is not 
practical. For examnle, both the equations for the two so flow and chemical 
kinetics losses become "stiff" in the equilibrium limit result in unacceptably 
long and difficult calculation^ procedures. Even if sophisticated numerical tech- 
niques are used to '-elculete the "near* equilibrium case, extrapolation to the full 
equilibrium solution introduces errors which negate any advantage gained by using 
tho same computer code. Hence, the following criteria was adopted in -aiculatlng 
* the reference performance for each of the losses (except hie boundary layer loss 
which is calculated directly). 

The reference method should be the most exact method of cal- 
culation available. The conceptual differences between the 
reference and os a models should Include only those mechan- 
isms which ^use the loss. 


1-4 


1 


TABLE l-l INTERACTION OF PHYSICAL PHENOMENA WITH 
PERFORMANCE LOSS CALCULATIONS 


^^PERFORMANCE 

^vLOSSES 

PHENOMEN^V^ 

Divergence 

Loss 

Boundary 

Layer 

Loss 

Kinetic 

Losses 

Two- Phase 
Flow 
Losses 

Combustion 

Inefficiency 

Loss 

Non 

One- Dimensional 
Flow 

* 


1 

2 

2 

3 

Viscous And 

Heat Transfer 

3 

a 

3 

2 

3 

Finite Rate 
Chemistry 

3 

2 

a 

3 

3 

Multiphase Flow 

1 

2 

2 

IK 

3 

Incomplete 

Combustion 

l 

3 

2 

2 

3 

3 


Log end: 

1 . Primary Importance (could be ■» 0.2* effect on I #p ) 

2. Secondary Importance (probaLiy <0.2* effect on l #p ) 

3. Generally Not Important 


l-S 


r 


































Table 1-1 (from Ref. 9 ), shows the interaction of the physical phenomena 
considered here with each of the performance losses. The table should be read 
down and then to the left. For example, the main interaction between the diver- 
gence loss and the various physical phenomena is that associated with multi- 
phased flow, while the boundary layer loss is chiefly affected by the nature of 
the specified boundary layer edge conditions, (i.e., from a one or two-dimensional 
flow field calculation). Table 1-2 Is a condensation of Table 1-1, and indicates 
only those cross coupling effects that were deemed strong enough to warrant con- 
sideration. 


Table 1 -2 Cross Coupling of Performance Losses 


Loss Mechanism 

Cross Coupling 

Finite Rate Chemistry {Kinetics) 

none 

Divergence Loss 

Two Phase Flow 

Two Phase Flow 

none 

Boundary Layer Loss 

Non One-Dimensional Flow 

Combustion Efficiency Loss 

none 


Table 1-2 shows that while the kinetics, two phase flow and combustion 
efficiency losses can be treated independently, the divergence and boundary layer 
losses are coupled to non-one-dlmensional and two phase Low phenomena, re- 
spectively. Hence, it was decided that the best course of action to take in the 
SPP code was to couple the two phase flow and divergence losses together and to 
u e the results of that calculation to obtain the edge conditions for the boundary 
layer calculation. 

With the above in mind, the following sub-sections describe the basic modules 
and assumptions used in the SPP code. 

1.2 Master Control Module (MC) 

The master control module (MC) controls the execution of the SPP computer 
code. The MC module selects, via user input, which of the five basic calculation 
modules (0DE, BAL, 0L'iC, TD2P, and TBL) are to be executed, and whether calcu- 
lated, input, or empirical, losses are to be selected for the calculation of delivered 
specific impulse, thrust, and total impulse. 
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This module handles or specifies all of the linkage communications bet- 
ween modules. The MC module selects which parameters (either calculated or 
input) are transmitted from one module to another and tries to determine the best 
overall delivered performance consistant with the input to the SPP computer pro- 
gram. 

Figure 1-1 shows the basic input and output of the MC module, while Table 
1-3 shows the variables used to calculate delivered thrust, I , and total impulse 
and indicates the module in which each is calculated. 


Table 1-3 Variables Used for Performance 
Prediction by the MC Module 


Module 

Variable 

0DE/One Dimensional Equilibrium 

Isp th.' Is Preo' f 

BAL/Graln Deslgi end Ballistics 

m, r*. L*, ry c# as functions of t 

0D&/One Dimensional Kinetics loss 

^ 

sp «N 

TD2P/Two Dimensional Two Phase flow loss 

I?p TD2P ' T ’ 7X,2P a * funct,ons of «,<’ and C D 

TBL/Turbulent Boundary layer loss 



The vacuum delivered specific impulse is calculated as in equation 1-3 
and is written in full here. 


! sp D = *sp th n KIN ' "TD2P ' • 77 SUB " M sp m (: ' 4) 

The theoretical specific Impulse, I , is rase-' on the initial nozzle 

exit area ratio. If a motor experiences throat e rosion during the course of a firing 
the nozzle exit area ratio, and hence, *'.c specific impulse, varies with time, 

Ihls gives rise to an additional loss mechanism; one which was only mentioned 
but not discussed In Section 1.1. The performance decrement associated with this 
phenomena is usually called the erosion loss. While this loss can be estimated 
based on one-dimensional calculations, it was felt that a two-dimensional esti- 
mate of this effect wculd be superior. Thus, tor eroding nozzles, Hj^p in 0( l ua ' 
tlon (1-4) is modified as follows: 

<<> 

— * : TD2P 

*>TD2P “ ^TD2P I T71 (1 " 5) 

sp TD2P * 
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FIGURE I -I MASTER CONTROL MODULE (MC) 



where 



r? TD2P 



When defined In this way r? Tn2 p includes the "erosion” loss in addition to 
the coupled two-phase and divergence losses. The erosion loss, as defined here- 
in, is equal to the ratio of the 2-D, 2-phase I at the time averaged expansion 

sp 

ratio, I (7), o the specific impulse evaluated at the initial expansion ratio, 

sp TD2P 

* sp TD2P ^ c i^ (Future modifications to the program should break the erosion loss 
out as a separate efficiency.) 


or 


The combustion efficiency, is defined as follows: 

^CE = C D ’ ^c* 

tj CE * ^ the previous relation yields a value 

greater than unity. 


This definition prevents the mass flow effect related to the nozzle discharge coefficient, 
C n , from being coun + ed twice since it is an integral pert of the 2D--2 phase loss. In 
the empirical performance prediction procedure is taken to be equal to tj c *. The 
empirical determination of c* efficiency is discussed ir. Volume 1 of this report. 


The loss due to finite rate chemical kinetics is treated as an efficiency 
factor, The Kiretlc efficiency is calculated ac the »at:o of an I calcula- 

ted by the QDK module, to an I calculated by a sp'Vu-1 option in the ODE module. 

sp 

This special ODE option, referred to as the "Restricteu Equilibrium" option causes 
the ODE program to compute an based on the more restrictive physical assump- 
tions employed in the ODK calculation. This procedure is accessary if certain 
losses arc not to be counted twice. The aforementioned assumptions are describ- 
ed in Section 1.5. 

The submergence efficiency, rjg^p/ rased on an empirical correlation 
is a function of the throat radius, mass fraction of condensed phase, length of sub- 
mergence to length of the Internal motor, ar.d chamber pressure. This correlation 
is discussed in more detail in Volume I of this report. 
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The decrement in performance due to the formation of a turbulent boundary 
layer in the motor nozzle is transmitted directly to the MC module and is reported 
as a function of area ratio. The decrement in I gp is taken Lo be the value corres- 
ponding to the last (l.e. the maximum) area ratio input by the user. This decrement 
is assumed to be independent of erosion induced area ratio changes (it is more 
a function of the wettsd length of the nozzle) and throat roughness. 

The equations used to calculate the "theoretical” performance losses, and 
the formulas used in the "empirical", or simplified, performance calculations are 
presented in Volume I and Volume II. 

1.3 Theoretical I Module (0DE) 

sp 

The theoretical I en module (ODE) is an outgrowth of the NASA Lewis chemical 
sp 

equilibrium computer program described in NASA SP-273 (Reference 1) and in the 
JANNAF standard TDK computer program manual (Reference 2 )• 

This module calculates the theoretical, or maximum, I which can be de- 

sp 

livered by a propellant at a given area ratio. The ODE module uses a free energy 

minimization technique along with the following assumption* to calculate the 

vacuum I . 

sp 

One dimensional flow 

Thermal and velocity equilibrium of the condensed phase with the gas 
phase flow 

This version o* the ODE program has also been modified to calculate the 
transport properties needed by the ballistics (SAL) module, two phase flow loss 
(TD2P) module, and the boundary layer loss (TBL) module. 

In addition, the program has been modified to calculate the reference con- 
dition for the loss in performance due to finite rate chemical kinetics. This option 
is referred to as the "restricted equilibrium" option and is discussed in more de- 
tail in section 1.5. 

Figure 1-2 illustrates the major functions of the ODE module as well as 
the Inputs and outputs of the module. 


1.4 Grain Design and Ballistics Module 

When utilized as part of the overall performance prediction methodology 



i I 


I I 


FIGURE 12 THEORETICAL ISP MODULE (0DEl 




the purpose of the Grain Design and Ballistics (GD&BC) module Is to provide values 
of pressure, mass flow rate, throat radius, L*, and C* efficiency, as functions of 
time. Time averaged values of these quantities are also computed and used in the 
other modules. Figuie 1-3 illustrates how the GD&BC module functions in this 
mode of operation. Various thermodynamic and transport properties are internally 
transmitted to the GD&BC module from the ODE module. Likewise, the nozzle 
geometry and ambient oressure history are obtained from the Master Control module. 
The final thrust versus time output of the programs in this integrated mode of oper- 
ation are based on the GD&BC calculated chamber pressure and mass flow rate; the 

delivered I is calculated from the combined results of the other modules, 
sp 

The GD&BC module can also be employed in a stand-alone mode. In this 
case, all of the parameters which would otherwise be supplied by the ODE and 
MC modules must be input directly into the Ballistics module. The thrust versus 
time curve in this mode of operation is based on the calculated mass flow rate, 
an empirical C* efficiency, and an input value of nozzle (C^) efficiency. 

The essential features of the GD&BC module are: 

The Hercules Inc. 64101 grain geometry analysis program, omitting non- 
essential elements for considerable size reduction, and modified to be com- 
patible with the Lockheed Propulsion Company interior ballistics analysis 
program. 

An improved version of the Lockheed Propulsion Company 637 interior ballis- 
tics analysis program. 

A correlation of C* efficiency. 

A correlation of nozzle throat erosion 

1.4.1 Gra in Geometry Program 

TheGrain Design Program examines a grain in three dimensions and solves 
for the surface versus bum distance relationship. The grains are not limited to 
families of grains such as star, slotted tube, etc., but will handle these and more 
complex grains. A detailed users manual for the original Heicules 64101 program 
exists as Rof. (5). The following extracts are essential portions for present purposes. 

1 . 4 . 1 . 1 Method of Approach 

The Grain Design Program simulates crafting techniques used in develop- 
ing the surface area versus bum distance curve of a solid propellant motor, given 
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an arbitrary number of symmetrical angular sectors. The program usee a right- 
handed orthogonal ;oordinate system in X, Y, Z where X is the longitudinal axis 
passing along the center lire of the motor and the Y and Z axes pass from the center 
of the motor to the outer boundary. 

Three basic input figures ar used in various combinations to describe the 
Initial propellant geometry. Thes' figures are: frustrum of a right circular cone, 
right triangular prism, and a sphere, A union of the input figures is performed by 
computing the intersection of a line Y*C and each of the inpuc figures. The over- 
lapping points are eliminated and those left represent the propellant geometry at 
the particular X, Y and bum coordinate consistent with the symmetry of the pro- 
pellant. 

Each figure may be input as consisting of void or propellant, where the 
void is burning outward and the propellant is burning Inward. Since the burning 
is assumed normal to the surface and the surface area is equal to the derivative of 
the volume of the propellant with respect to the normal, the surface is computed 
as the derivative o' the volume of the propellant with respect to the bum distance. 
Volume can be computed as the integral of the propellant areas at each X station 
with respect to X, and propellant area can be expressed as the Integral of the 
lengths across the propellant with respect to Y, 

1 • 4 • 1 . 2 General Description of the Sol ut ion 

The besic maus balance equation 

m * s r ? (1-7) 

is the basis for the solution. By eliminating time and density, equation (1-7) 
ry*comes 

s ■ dv/db - v/\b (1-8) 

and from the mean value theorem s will occur srmewhere in the Interval *b. The 
half-way point was atbif’erlly accepted, 
m * mass 
r » rate of bum 
o * mass density 
s * burning surface area 
v * volume of propellant 
b * bum distance 


1-14 


In the solution for b versus s the proc r am . >•-; not make any direct analytic 
solution for a burning surface, but rather makes c.d out ions for the volume of 
propellant at each bum distance specified In the ir;p , r cords the change in volume 
between two consecutive bum distances, and divu. : y the difference in the two 
burn distances. Thus, the user should keep in mlr* > t he is preparing input 
to simulate the chenges in volume which will om-i ; i the motor, not to ciiroctly 
compute burning surface areas. 

1.4. 1.3 Method o: Determining Volumes: Confic on Simulation 

The case is assumed full of propellant irb<; ; , The voids in the case 
are simulated by figures - cones, cylinders, prism.:- - r pheres - whlefc may over- 
lap and/or protrude outside the case, if needed. b ^ouid simulate the void 
as it is initially. The figures may be normal out ■ . : c ; cures , grain -fill inbum- 
ing figures, or nonourning figures and may have r corners and .edges, if needed. 

The order of input can be important. Each figure m. ; *y be placed in any orientation 
anywhere in space, inside or outside of the gr~i • nf these is described 

in detail later. Each figure is located in thrr- dt.r u d space by the follow .no 

means. 
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Spheres: 

V V 2 i 


Generally, the grain will be symmetric about the axis of the motor. If 
there are symmetries in the grain they should be reflected in the input in order to 
save work and computer time. See Figure 1-4. The symmetries are automatically 
accounted for in the whole motor. 


Figure 1-4: 


Sample symmetries: (End views of motors) 


SX «* 6.0 SX - 3.0 

6-60* A - 120* 

(Tic X axis is perpendicular to the paper) 


SX - 2.0 
A - 160* 



1,4*2 Interior Ba, listlcs 


The Interior Ballistics Program models the internal ballistics of solid rocket 
motors, including motor environmental effects on burning rate (erosion, radiation), 
effects of radial slots betweeu segments on the port flow port flow pressure drop, 
and delayed burning (ignition) in deep narrow slots. Semi-empirical correlations 
for combustion inefficiency and nozzle throat erosion also are included. This 
digital computer program was adapted from an existing model used to analyze nozzle- 
less solid rocket motors ^ ® K It was selected because it contains features lacking 
in many published internal ballistics programs. The Intent was to predict internal 
ballistics for motors ranging from small tactical motors , to large boosters , wherein 
aspects of all the aforementioned features are encountered. The major effort in 
the present program development involved updating certain model expressions, and 
rendering the computational procedure compatible with the Grain Design 
Module and the balance of the SPP Program Modules. 


1.5 Finite Rate Chemical Kinetics Loss Module (gPK) 

In most rocket motors there is a loss of performance due to the fact that 
a finite time is required for chemical reactions to proceed to the maximum amount 
of sensible heat release. Normally, this loss is causi i by the incomplete recom- 
bination in the rocket nozzle of the highly energetic species which are formed during 
combustion in the rocket motor. 

The 9DX module only considers the losses which are due to gas phase 
chemical kinetics. Other losses which include incomplete or partial solidification 
of condensed phases are ignored in the ODK module since adequate models of +hese 
processes have yet to be formulated. The ODK modulo also assumes that the flow 
is one dimensional, with the gat and condensed phases in dynamic and thermal 
equilibrium with each other. No transfer of mass between 'ms phase species and 
the condensed phase species are allowed. The liquid phase, however, 4s allowed 
to undergo solidification when the gas temperature reaches the melting temperature 
of the condensed phase. 

In order to properly calculate the chemical kir.ctlc performance efficiency, 
^KIN* a r *for*nce equilibrium I 0p which is baseo on the same assumptions as ODK 
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r 


should be used. In ordqr to provide this reference value; the 0DE module was modi- 

fled to ca’culate a "restricted equilibrium" specific impulse I . This modifl- 

sp RE 

cation to the ODE module allows the available gas phase species to remain In 
chemical equilibrium with themselves while restricting any mass transfer to the 
condensed phases. Hence, the kinetic performance efficiency, which Is defined as 

"KIN " '.Pqd/.Pre 

would approach unity as the chemical kinetic reaction rate constants approach In- 
finity. 

Another advantage of the use of the restricted equilibrium option for computing 
kinetic efficiencies Is that only thermodynamically Important species need be con- 
sidered In the ODK-ro calculations. Thus, trace species, and reactions Involv- 
ing them, need not be included in this calculation. 

i Figure 1-5 suows the basic Input, output and assumptions used In the ODK 

module. A comprehensive set at reactions for aluminized propellants Is also In- 
cluded In the volume. (See Table 2-13) Appendix A of Volume I Indicates how the 
reactions and species where selected for aluminized propellants. 


1.6 Two Dimensional -Two Ph ase Loss Module (TD2P) 

Almost all rocket nozzle* of practical Interest have a lo ss in performance 
due to a nonaxtal component of velocity at the nozzle exit plena. This loss Is 
usually referred to as a divergence, or two dimensional flow, loss. If the pro- 
pellant Is metal txed, there Is an additional loss when the particles cannot main- 
tain a state of dynamic and thermal equilibrium with the gas flow. This loss 
is referred to as the two phase flow loss. 

There are two distinct <-oupllng effects between the two dimensional and 
two phase flow losses. The first Is that If the gas and conurnsed phase are not 
In dynamic (velocity) equilibrium with each other then the gas and particle stream- 
lines are not coincident . This can have a noticeable effect on the divergence 
loss. Ihe second coupling effect is on the mess flux through the nozzle. It is 
now well known that the mess flow through e nozzle for single phase two dimen- 
sional flow Is less then the amount of mass flow that can be achieved et a given 
chamber pressure for an equivalent one dimensional nozzle Row. This effect results In 
nozzle discharge coefficients being less th in unity. However, it Is also well 
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known that in nonequilibrium two phase flow more mass flows through a nozzle 
at a given chamber pressure than would if the particles were in velocity and 
thermal equilibrium. By itself, this effect would produce values greater than 
unity. For most metalized solid propellant motors the net result of these two op- 
posing effects is & discharge coefficient greater than unity. If the combustion 
efficiency were 100% (complete combustion) the C* efficiency of the motor, 77 

c 

is given by 

Thus, in most cases, tj^would be less than unity for metalized propellants, and 
would always be greater than unity for non-metallzed propellants*. As a result of 
the relationship between t) c *nd C^, the C D effect must be corrected out of the 
empirically obtained c* efficiency if one is to obtain a value for combustion 
efficiency which does not count this factor twice. 

The TD2P module is a modified version of the computer program developed 
by Kllegel and Nickerson (Ref. 3 ). This module, which has been extensively 
modernized, calculates the losses which are usually the largest encountered in 
most solid rocket motors. These losses are, the two phase flow loss, erosion 
loss, and divergence loss. 

The basic inputs, characteristics, and outputs of the TD2P module are 
shown in Figure l-(*. Since this is the most important loss module in the SPP 
code it is felt that the following comments are in order. 


♦Most solid rocket motors have normalized radii of curvature ^2 and, In the ab- 
sence of particles Cjs tends to be .99 < C D < 1 , or, in other words, with perfect 
combustion ti c ^woulcroe bounded as follows, 1 < £ 1 .01 . This small effect 

tends to be masked by measurement and data reduction errors, and many times 
test engineers are loathe to report values of greater than unity due to the 
mistaken assumption that such values are physically impossible to achieve. 
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TWO DIMENSIONAL TWO PHASE FLOW CALCULATION MODULE (TD2P) 




1) The TD2P module assumes that the gas phase flow behaves as an ideal 
gas. The gas phase properties are calculated by matching ideal one di- 
mensional "equivalent" gas results to the output of the ODE module at 
the solidification point. In order that the condensed phase properties, 
and, hence, the results of the TD2P module, would not be artificially de- 
pendent on aiea ratio they were taken as the actual values obtained from 
the 0DE module. 


2) The transport properties which are transmitted from the 0DE module 
are selected to insure the greatest degree of accuracy in the throat region. 


3) The reference I sp (one dimensional flow with zero lags, I sp ), to 

which the TD2P I S p is ratioed, is calculated using the same idJJi^ks properties 
as are used in the TD2P analysis. By using a ratio calculated in this man- 
ner, Instead of using the absolute I sp number from TD2P, the effect of 
neglecting real gas effects is minimized. Also, extra care was taken to 
accurately caictriate the reference I during solidification so as to avoid 
the use of different assumptions in fne two modules. 


4) The fact that the two phase flow loss tends to decrease with increased 
throat diameter (erosion) was intentionally ignored, since it was felt that 
the added loss due to throat roughness would more than compensate for this 
effect. 


5) The transonic analysis currently incorporated int3 the TD2P module is 
restricted to running normalized throat radii of curvature of no less than 
approximately 1.5 and inlet angles less than about 45°. In addition the 
program requires the throat geometry to be circular with equal upstream 
and downstream radii of curvature. These restrictions on the geometries 
which may be realistically treated are a major shortcoming of the TD2P 
module. 

6) The edge conditions which are generated by the TD2P module for the 
turbulent boundary layer, THL, module have been adjusted so that the sub- 
sonic and transonic regions will faired into the supersonic regions in a 
smooth manner. The stagnation conditions transmitted to the TBL module 
have also been adjusted so as to allow the TBL module to reproduce the 
TD2P edge conditions (velocity, pressure) as closely as possible. 


1.7 Turbulent Boundary Laver Module (TBL) 

The losses associated with the presence of a boundary layer result from 
a transfer of heat to tne nozzle walls and from a reduction in momentum due to shear 
stress along the walls. These losses can be calculated by properly integrating the 
stress tensor along the wall (pressure plus skin friction), or by calculating the 
momentum deficit in the boundary layer at the nozzle exit plane. The second 
method is the easiest to apply, and is the JANNAF recommended procedure. 

The boundary layer development in the nozzle is currently computed using 
the Bartz-Silver integral method, as embodied in the JANNAF adopted THL 
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program^ 4 \ Thl? is not the most accurate boundary layer program available, 
however, Its virtues (speed and simplicity), coupled with the usually small 
magnitude of the loss (usually less tha 1% for all but small motors), make it a 
reasonable choice for the present purpose. 

The boundaiy layer solution requires a specification of the inviscld flow 
along the wall streamline. In the SPP program these edge conditions are supplied 
automatically by the TD2P module (if executed). This allows the interaction be- 
tween two dimensional and boundary layer effects to be properly accounted for. 

The TD2P module supplies the TBL module with the nozzle geometry, Mach number 
along the wall streamline, ideal gas properties, and stagnation pressure and temp- 
erature. This procedure allows the remaining edge properties to be computed in- 
ternally in the TBL module (pressure, velocity, etc.) to match, as closely as pos- 
sible, the respective TD2P computed counterparts. 

The boundary layer loss is computed directly as an I decrement, rather 

sp 

than as an efficiency, using the following formula 


tX 


sp 


BL 


(2ir r ^u a Q cos a ) & 
9 m XD2P 



i 


in which the various quantities are all evaluated at the nozzle exit plane. It 
should be pointed out that this equation gives the complete boundary layer loss. 

The reduction in fluid momentum due to shear and heat losses to the wall are both 
properly accounted for. 

In most cases the effect of the ir.vicid -viscous interaction between the TD2P 
and TBL modules is small and does not warrent the subsequent Iteration which re- 
quires the re-running of TD2P module using the new Invlcid wall contour supplied by 
the TBL module. However, in extrem e cases, the above procedure may be Implemented 
by using the invlcid wall contours supplied by the TBL module as input to second 
TD2P-TBL calculation. However, when boundary layer losses are of such importance 
that such effects should be considered, it is doubtful that the TBL analysis, and 
the assumptions inco&porated within, are valid. Therefore, is is suggested that in 
these cases, that a more accurate method or analysis of the turbulent boundary layer 
loss be used. 
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OUTPUTS 



FIGURE 1-7 TURBULENT BOUNDARY LAYER LOSS MODULE (TBL) 



That Is to S3y, when the losses are small ( u 1/2%) even a poor character- 
ization of the loss (± 50%) is acceptable. When the losses due to a turbulent 
boundary layer are large, then secondary effects (lnvlcid-vicid flow interactions) 
are small compared to the primary boundary layer analysis loss. Hence, if the 
boundary layer loss is significant corrections made to the THL module method of 
calculation are of minor Importance since the overall method selected is of 
doubtful accuracy. 

The input, Internal linkages, characteristics, and output of the TBL module 
are schematically portrayed in Figure 1-7. 
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2. INPUT DESCRIPTION AND LINKAGE 



The SPP Program is structured so that the various program elements which 
comprise the total performance prediction methodology can be utilized as an integra- 
ted set, as individual program^ in a stand-alone mode, or in any user specified 
combination. The individual program modules have been combined with a control 
module which provides a means for inter-module communication. This allows 
the complete methodology to be executed automatically, with a minimum of user 
prepared input. Minimum is used here in the relative sense, as the amount of 
input required to perform all the analyses is not small. 

Table 2-1 is essentially a block description and flow chart of the input 
data requirements of the SPP Program. In order to exercise the total program 
capability in a single run, one would have to consider tne input requirements of 
all of the data sets listed in Table 2-1. The notes column of the Table indicates 
which data sets can be skipped if one desires to execute only a limited number 
of the program options. The subsection in which each of the data sets is de- 
scribed is also indicated in Table 2-1. 

Table 2-1 should be used as a guide when preparing input for a given pro- 
blem. It lists the data sets in the order in which they should appear in the data 
deck, and also shows the special cards which must appear in each set (first card, 
last card, etc.) if the program is to function properly. Iabie 2-1 is basically 
self-explanatory when used together with the detailed input descriptions which 
follow. The following additional notes are parenthetical. 

1. Thermodynamic data will normally be on tape or on a mass storage 
permanent file. Card input will be required only when species not 
residing on the tape are considered. 

2. Thermodynamic data below 300 g K will rarely be required for solid rocket 
motor problems. 

3. The Geometry data set is not optional if one desires the post-module 
execution summary outputs, or the final performance summary tables. 

4. Figures 2* 1 and 2-2 of Volume II can be very useful in understanding 
how and in what order data is transmitted between modules. 

Section 2.11, following the input data descriptions, describes the Internal 
(and external, via punch cards, perm, files, etc.) linkages between the basic 
computational modules of the program. 
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Table 2-1 Input Data Set Description 


Data Sot 


Title c.trda 


Section EaplatnU* 
Date Set 


card col. no* 

)|2 3 4 sU 7 I 9 10 


first card t 

1 TLE title formation 

lust card 

nono required 


Tbermo.iyn.vnlc Data 


lirst card T 

I 1 E RK 10 

last card 1 

N D 


Thermo Jynamtc Data 
Selow 30 d*k 


first card L 

HW T CP HS 

last card l 

ND LQIW T CPHS 


Geometry card 


first cerd C 

t (BMC T* Y 

second card 

$ GE C M j 

last card 

SEND 


Trejociory card 

Urtt cerd "* 
second card 
last card 
Problem card 

lira! cerd **" 
second card 
laat card 

* 

Seeetents card 


T R A J C CTgRV 

Z mt i 

_ $CNC 

*P R >Bl E M 
"SPRCB 
“ ~nT 


first card R 

C ACT A N 1 S 

last card 

blank card required 


On It end Insert card* 










( 9 DE Namelist 


ftrat cerd 


second cerd 
last card 






























Table 2-1 (continued) 


0D6 Namelists 




first card 

E 

AMFL 

1 STS 


second card 




lest card 

$ E NO 


0DK option for input of 
Initial Species Concentre* 
tlons 

■ 




| first card 

3 

l» E Cf 

E S 



|iasi card 


none required 

enactions 





first card 

R 

QQES 

I ONS 



last card 

I 

OSH 

REAX 

Inert Species Option 

■ 

• 



(first card 

D 

ffiii 

S 



(lest card 


none required 

Third body Reaction 

Rate Ratios 

■ 




first card 

*< 

9 

irn r r 

90DY 



le»; card 

l 

ni 

CARD 

0DK Namelist 


S 




first card 

graM 




lest card 


S C NO 


102 F Namelist 


.... 




first card 


$ TD2 




lest card 

K NL 


TIL KamoMst | 





liras card 

i> ^ 9 t 




lost -aid 

j 

$ C 4 0 









Kflri 

I Opt tor il. 

I More then one 

) Mey eppeer Ir. any or.i *.■< before die PR09M.M r«rd. 

4 Oi'tiom.1 , U « r< .^ny: iMte dete already ^uti, 

% Xtundut^y. 

( fc-qtilM.I only «f t'.* c*: i* uj Mo* is select*-.!. 

i d it s*.f t*Ai option ia 

I l« qulre.1 inly If r- • ’n< *• 'I V'quitlhnufn option *»!ect**H. 
9 only If t* • .*P»: <h Mon is selected, 

10 fcquli'.tl only It •?.« l i */l* option Is select* d. 

11 Neqtiltcd only If »■*« TIU. option l» selteled. 

It It requested by irput In iSt WW9 eee*|ist. previously 
tenanted )u**q **y be Included after tHe PK» 

eaettltst set (see Section J JO, 


2.1.4 

2.5.5 


2.9.1 


2.9.2 


2.9.1 


1.9.4 


2.9.1 


1.9 


t.l* 
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TITLE CARDS 


This input permits labeling of runs with alphanumeric information. As 
many title cards as desired may be Input in sequence. Card fonnat is as follows: 
col 1-5 col 6-77 

TITLE any alphanumeric Information 


2-4 


2.2 


THERMODYNAMIC DATA 


This data set Is identical to the THERMO DATA described in Appendix 
D of NASA SP-273 (i.e. Reference 1 ). 


Using this data set, thermodynamic data curve fit coefficients may be 
read from cards. The curve fit coefficients are generated by the PAC computer 
program described in NASA TN D-4097 (i.e. Reference 7 ). 


The thermodynamic data (i.e. 



etc.) are expressed as functions of 


temperature using 5 least squares curve fit coefficients (a^ 5 ) and two integration 
constants (a g7 ) as follows: 




Rt 


a j + a 2 T + ajT* ♦ a^T* ♦ a g T 4 




( 2 . 2 - 1 ) 

( 2 . 2 - 2 ) 


S® a,T* a.T* a^T 4 

— jj— » aj lnT + a 2 T + -^~ +-|— +a 7 (2.2-3) 


For each species, two sets of coefficients (a^_ 7 and a j_ 7 ) are specified 
for two adjacent temperature intervals, lower and upper respectively. For the 
data available in Reference: 1 the lower temperature Interval is 300° te 1000°K 
and the upper temperature interval is 100(P K to S000°K. 

The input format required for this thermodynamic data Is defined in Table 
2-2. Data cards for the species H. H02, H2, H20CL), H20, H202, O, OH, 02. 
A1203(L)« and AL203(s) are listed in Table 2-3 as examples. Thermodynamic data 
coefficients for 42 1 species are available in Reference 1 and are supplied with 
the computer program. A list of these 421 species names is presented in Table 
2-4. 


Data Tape Generation and Usage: 

A computer run using thermodynamic data card Input will generate a data 
tape on logical unit JANAF. This tape may then be saved and used at a later time. 
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The program writes the THERMO data card images on unit JANAF as read but with 
two minor exceptions. The THERMO code card and the card numbers in card 
column 80 are omitted. 

If thermodynamic data cards are not input, the program assumes the 
thermodynamic data is on logical unit JANAF. Logical unit JANAF is currently 
assigned a value of 25. 
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TABLE 2-2 


FORMAT TOR THERMODYNAMIC DATA CARDS 


Card 

order 

Contents 

Format 

Cord 

column 

1 

THERMO 

3A4 

1 to 6 

2 

Temperature ranges for 2 sets of co- 
efficients: lowest T, common T, and 
highest T 

3F10.3 

1 to 30 

3 

Species name 

3A4 

1 to 12 


Date 

2 A3 

19 to 24 


Atomic symbols and formula 

4(A2,F3.0) 

25 to 44 


Phase of species (S,L, or G for solid, 
liquid, or gas, respectively 

A1 

45 


Tempi rature range 

2F10.3 

46 to 65 


Integer 1 

115 

80 

4 

Coefficients a Mi * 1 to 5) in equations 
(2*2-1) to (2.2 i -3) 

S(E1S.8) 

1 to 75 


(for upper temperature interval) 

Integer 2 

15 

80 

5 

Coefficients in equations (2*2-l}-(2«2-3! 
(ai, ei for upper temperature interval 
ana a^' am * °3 * or lower) 

S(C15.8) 

1 to 75 


Intcycr * 

15 

80 

6 

Coefficients In equations (2,2-l)-{2*2-3 
(a*,Oe#a-,a, for lower temperature In- 
terval! 6 7 

4(E15.8) 

1 to 60 


Integer 4 

* CO 

60 

(a) 

Repeat cards numbered 1 to 4 in cc 80 
for each species 



(Final 

card) 

END (Indicates end of thermodynamic 
data) 

3A4 

1 to 3 


*Gascou$ species end condensed species with only one condensed phaso can bo 
in any order* How over, the sets ter two or more condcnsod phases of the same 
species must be adjacent* If there arc more than two condonscd phases of a 
species their rets must be either ht increasing or decreasing order according 
to their temperature intervals. 
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TABLE 2-3 . SAMPLE THERMO DATA CARDS 


{Species Aft, fe, H 2 , HjO. N 2 , O, OH, 0 2 ) 


THERMO 

300,000 1000 

AR 

0,250000006 01 
•0,745375026 03 i 
0. i 

H 

0,250000006 01 i 
0,254710276 05- 
0, 

HR 

0.310019016 01 
-0,577300926 03- 
0,552103916-00- 

"f^l 676336 01 
-0,299050246 05 
-0,290374046-08 
N2 

0,209631946 01 
-0,905041046 03 
-0,432175596-09- 
0 

0,254205946 01- 
0,292300036 05 
-0,140204326-00 
OH 

0,29)044276 01 
0,39353015E 04 
0,107139726-09- 
02 

0.362195356 01 
-0,120190256 04 
*0,676351 376-00 
FNO 


,000 5000,1 
L 5/66AR 

0, 

0,436400046 

•• 

J 9/05H 

0. 

0,460117636 

0, 

J 3/61H 
0,511 194646- 
0,19629*216 
0,101227396' 

0,^94513746 
0,663054716 
0,8570?! C3C 
J 9/6SN 
0,15154**66 
0.616151406 
0,22577;.'51E 
J 6/626 
>0,27650096 
0. 49203ft 006 
0*309069646 
J 3/660 
0,959316506 
0,544234456 

•0,225710946 
J 9/650 
0.736102446 
0,361509606 
0,215559936 


000 

100 000 000 
0* 

01 0,250000006 
•0,745374906 
100 000 000 
0* 

00 0. 200000006 
0.2*4716276 
20 00 00 

•03 0.526442106' 
01 0.305744516 
-11-0,980904746 
20 100 000 

-02-0.002243746' 
01 0,407012756 
-12-0,35?797?2C 
20 00 00 

-02-0.572352776' 
: 0) 0*3674*2616 
-!?-0*106U5*AE 
140 000 000 

-04-0.310200336' 
: 01 0*294642*76 
;-12 0*291476446 
)H 100 000 

-03-0,194417026 
. 01 0,303759436 
-12 0*36412*236 
20 00 OO 

>03-0*196522206 
: 01 0*362559056 
>11*0,104752266 


06 300,000 5000.000 

0 . 0 . 

01 0, 0, 

03 0,430600066 01 

00 300.000 5000*000 

0 , 0 , 

01 0. 0, 

05*0,460117626 00 

06 300,000 5000.000 

•0 7-0 , 3490997 3E- 1 0 0, 369453456-14 
01 0,26765200£-02-0,56099i6?6-05 
03-0.22997056E 01 
06 300,000 5000,000 

•06 0,1 02266026-09-0,404721456-1 4 
01-0,1 10844096-02 0,415211006-05 
05-0,322700466 00 
06 300,000 5000,000 

•06 0,996073936-10-0,652235556-14 
01-0.120815006-02 0,232401026-05 

04 0.235004246 01 

06 300.000 5000.000 

•00 0,45S10674E-11*0,43600515E-15 
01-0*163016656-02 0*242103106-05 

05 0*296399496 01 

06 300*000 5000,000 

>06 0,137566466-10 9.142245426-15 
01-0,107780506-02 0,960303706-06 
04 0,493700096 00 
06 300*000 5000*000 

•06 0, 3620 1 5506- 1 0- 9 • 269456276-14 
01*0,107021046-02 0.705545446-05 
04 0,430527706 01 
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Table 2 -4. SPECIES WITH THERMODYNAMIC DATA PROVIDED 


AL(S) 

BCL 

BEO(S) 

AL(L) 

BCL+ 

BEO(L) 

AL 

BCLF 

BEO 

AL+ 

BCL 2 

BEOH 

ALB02 

BCL2+ 

BEOH+ 

ALCL 

BCL2- 

BE02H2 

ALCL+ 

BCL3 

BE 20 

ALCLF 

BF 

BE20F2 

ALCLF2 

BF2 . 

BE202 

ALCL 2 

BF2+ 

BE303 

ALCL 2+ 

* BF2- 

BE404 

ALCL 2- 

BF3 

BR 

ALCL2F 

BH 

BR2(L) 

ALGL3(S) 

BHF2 

BR2 

ALCL3(L) 

BH2 

C(S) 

ALCL 3 

BH3 

C 

ALF 

BN(S) 

C+ 

ALF+ 

BN 

C- 

ALF2 

BO 

CCL 

ALF2+ 

BOCL 

CCL2 

ALF2- 

BOF 

CCL3 

ALF3 (S) 

BOF2 

CCL4 

ALF3 (S) 

B02 

CF 

ALF3 

B02- 

CF2 

ALH 

BS 

CF3 

aln(g) 

B2 

Of H 

ALN 

B20 

CH 

ALO 

B202 

CH2 

ALO+ 

B203(L) 

CH20 

ALOCL 

B203 

CH3 

ALOF 

B303CL? 

CH4 

ALOH 

B303F3 

CN 

ALOH+ 

BE(S) 

CN+ 

ALOH- 

BE(L) 

CN- 

AL02 

BE 

CN2 

AL02- 

BE+ 

CO 

AL02H 

BEB02 

CCCL 

AL2CL6 

BECL 

CCCL2 

AL2F6 

BECL+ 

CCF 

AL20 

BECLF 

C0F2 

AL20+ 

BECL2 (S) 

COS 

AL202 

BECL 2 (L) 

C02 . 

AL202+ 

BECL2 

C02- 

AL203(S) 

BEF 

CP 

AL203(L) 

BEF2(o) 

CS 

AR 

BEF2(S) 

CS2 

AR+ 

BEF 2 (L) 

C2 

B(S) 

BEF2 

C2- 

B(L) 

BEH 

C2CL2 

B 

BEH+ 

cm 

B+ 

BEN 

C2F4 


C2H 

FECL 2 (S) 

KOH 

C2HF 

TECL2(L) 

KOH(S) 

C2H2 

FECL2 

KOH(L) 

C2H4 

FECL3 (S) 

K2 

C2N 

FECL J(L) 

K20(S) 

C2N2 

FECL 3 

LI(S) 

C20 

FEO(S) 

LI(L) 

C3 

FEO(L) 

LI 

C302 

FEO 

LI+ 

C4 

FE02H2(S) 

LICL(S) 

C5 

FE02H2 

LICL(L) 

CL 

FE03H3(S) 

LICL 

CL+ 

FE203(S) 

LIF(S) 

CL- 

FE304(S) 

LIF(L) 

CLCN 

H 

LIF 

CLF 

H+ 

LIF2- 

CLF3 

H- 

LIFO 

CLO 

HALO 

LIH(S) 

CL02 

HBO 

LIH(L) 

CL2 

HBO+ 

LIH 

CL20 

HB02 

UN 

CS(S) 

HCL 

LIO 

CS(L) 

HCN 

LIO- 

CS 

HCO 

LIOH(S) 

CS+ 

HCO+ 

LIOHOJ 

CSCL(S) 


LIOH 

CSCL(S) 

HF 

LION 

CSCL(L) 

HNO 

LI2 

CSCL 

H02 

LI2CL2 

CSF(S) 

H2 

U2F2 

csp(l) 

H20(S) 

LI20(S) 

CSF 

1120(1.) 

LI20a) 

CSO 

H>0 

LUO 

CS2 

H202 

LI202 

CS2CL2 

H2o 

I.I202H2 

CS2F2 

U3B306 

LI3CL3 

CS2C 

HE 

LI3F3 

E 

HE+ 

MG(S) 

r 

K(kj) 

MG(L) 

r- 

K(l.) 

MG 

FCN 

K 

MG+ 

FO 

K+ 

MGCL 

ro2 

KCL(S) 

MGCL+ 

rz 

KCL(L) 

MGCLF 

F20 

KCL 

MGCL2(S) 

FE(S) 

Kr(3) 

MGCL2(L) 

FE(S) 

K1'(I.) 

MGCL? 

FE(S) 

KF 

MGF 

FE(L) 

KF2- 

MGF2(S) 

FE 

K2F2 

MGF2(L) 

FECL 

KO 

MGF2 
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Table 2-4. (cont'd) 


MGH O 


MGN 

0+ 

MGO(S) 

0- 

MGO(L) 

OH 

MGO 

0H+ 

MGOH 

OH- 

MGOH+ 

02 . 

MG02H2 

02 - 

N 

P 

NF 

P(S) 

NF2 

P+ 

NF3 

PCL3 

NH 

PF3 

NH2 

PF5 

NH3 

PH 

(NO 

PH3 

NO+ 

PN 

NOCL 

PO 

NOF 

PS 

•NOF3 

P2 

N02 

P4 

N02- 

S(S) 

N02CL 

S(L) 

NO?F 

S 

N2 

S+ 

N2C 

SF4 

N2H4 

SF6 

N20 

SH 

N204 

SN 

NA(S) 

SO 

N/a) 

S0F2 

NA 

S02 

NA+ 

S02F2 

NACL(S) 

.S03 

NACL(L) 

S2 

NACL 

SI(S) 

NAF(S) 

SI(L) 

NAF(L) 

SI 

NAF 

SI+ 

NAF2- 

SIC 

NAH 

SIC2 

NAO 

SI CL 

NAO- 

SICL2 

NAOH(S) 

SI Cl. 3 

NAOH(L) 

sia4 

NAOH 

SIF 

NA2 

SIF2 

NA2CL2 

sirs 

NA2F2 

SIF4 

NE 

SIH 

NE+ 

SIH4 


SIN 

SIO 

SI02 (S) 
SI02 (S) 
SI 02 (S) 
S102 (L) 
SI02 
SIS 

512 
SI2C 
SI2N 

513 
XE 
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2.2.1 THERMODYN A MIC DATA BELOW 300°K 

i 

Tor low temperature calculations It may be necessary to extend the 
curve fit data In the Thermodynamic Data file (see Section 2.2)* The lower 
temperature limit, T^, in the Thermodynamic Data supplied with the program is 

T 300° K 
* 

Thermodynamic Data below the temperature, T # , may be input by data cards as | 

described below. 

col U 


card 1 L0W T CPHS I Directive for start of 

I low temperature CPHS 
tables (col 1 through 10), 

card 2 n 12 character species 


name, left justified 
followed by the integer, 
n, punched in column 21. 

The integer n must be such 
that 1 ^ns3 and represent 
the number of Thermodynamic 
Data points to be input for 
this species. 

card 3 T^K Cp 1 f irst xnermodynamic Data 

T point for the above species, 

• input 4F 10.0, 15. 

card n +2 T°K Cp H^, S^, 2 nth (l^n^3) nth Thermodynamic Data 

T point for the above species. 

Input 4JT 10.0, 15. 

. . . Repeat cards 2 through n + 2 above for each species to be input. 

Temperature must be T ^ < T 2 < T 3 < T 

(find END L01V T CPHS end directive (col 1 

card) through 14) 

An example of this input is given in Table 2-5 which shows a card listing extend- 
ing the Thermodynamic Data for an 0 2 /H 2 propellant to 100°K. Data in Table 2-5 
is taken directly from tne JANAF tables (Reference 8), except for Argon which is 
taken from NASA SP-3001. 

The quantity H^, is defined os 

Hi *» (H° -Ht 0 J - AHj , cal/molc 
T 290 *298 

Ci , cal/molc - deg K 
S!p , cal/mole - deg K 
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TABLE 2-5 . LOW TEMPERATURE C^, , 1^,, DATA FOR AN 0 2 /H 2 PROPELLANT 


LOW T CPHS 
AR 


2 


1 

100.0 

4.9661 

-984.5 

31.556 

200.0 

H 

4. 9681 

-487.7 

2 

34,999 

2 

100.0 

4.968 

51118. 

?1 .965 

1 

200.0 

M2 

4,968 

51614, 

2 

2S.40B 

2 

100.0 

5.393 

-1265. 

24.387 

1 

200.0 

H20 

6.518 

-562.0 

2 

?8.520 

2 

100.0 

7.961 

-59378,9 

36.396 

1 

200.0 

N2 

7,969 

-58581,9 

2 

41.916 

2 

100.0 

7,074 

-1387.0 

38,113 

1 

200,0 

o 

6.989 

•684. 

? 

*2.986 

2 

100.0 

5.666 

58479. 

32,466 

1 

200.0 

OH 

5.434 

59036. 

2 

36,340 

2 

100.0 

7.567 

7879, 

35.852 

1 

200,0 

02 

7,309 

8623. 

2 

41.021 

2 

100.0 

6.958 

-1381. 

41.395 

1 

200*0 
end l°>' 1 

6.961 

CPHS 

-685. 

*6.218 

2 
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2.3 Geometry Directive and Input 


Hie GE0METRY directive allows the user to input the nozzle geometry in 
the $GE0M namelist set and have that information transmitted to the appropriate 
calculation modules. Any or all of the values input via the GE0METRY directive 
may be overridden in the input to a specific calculation module. 

The format of the GE0METRY directive and description of the input is as 
follows: 

The word GE0METRY must start in card column 1 (Only the first 4 charact- 
ers are checked). The card following the GE0METRY card should contain $GE0M 
in card columns 2 through 6. The cards following the $GE0M card may be in any 
order up to the $END card which specifies the end of the geometry input data. 

All cards in the namelist input set must start in card column 2 or greater. The 
input items are 


Item 


Lo script ion 

Units 

Assumed 

Value(s) 

$GE0M 


namelist name 



ASUP(l) 

8 

supersonic area ratios at 
which information will be 
printed out 

none 

0.0 

NAS UP 

* 

number of entries in the 

AS UP array * 40 

none 

0 

ASUB(l) 

9 

subsonic area ratios at 
which Information will be 
printed out 

none 

0.0 

NASUB 

m 

number of entries in the 
ASUB array **10 

none 

0 

NN0Z 

a 

number of nozzles 

none 

1 

RSI 

a 

initial throat radius, r^*. 
normalized throat radius of 

in 

0.0 

RCURV 

a 

none 

0.0 

1 


curvature, r^r*. 



RSD0T 

•m 

nozzle erosion rate, r*. 

mils/sec 

0.0 

+ TBURN 

m 

motor bum time, used to 

sec 

0.0 


calculate r* as a functifn 
of time, i.e., r*(t)»r*+r*t. 


+ These values will be overridden by the Grain Design & Ballistics Module cal- 
culated values if BAL®1 or 2 in $PROB Namelist. 
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Item 


Description 

Units 

Assumed 
Value (s) 

IWALL 

= 

wall specification option 
flag 

none 

0 


1 cone option (input THETA, 

RCURV, THETAI, RSI, and 
EPS) 


2 parabolic nozzle contour 
option (input THETA, RCURV, 

THETAI, RSI, RMAX, ZMAX) 

3 circular arc nozzle contour 
option (input-same as in 
IWALL-2) 

4 nozzle contour (spline) option 
(input THETA, RCURV, THETAI, 

THE, RSI, RS, ZS, NWS) 

♦♦Note** - Certain of the calculation modules ((JDK and TD2P in particular) con- 
tain more extensive nozzle geometry specification options than are presented 
here. However, the options cited above are applicable to all of the calculation 
modules in the SPP code. 


The following items are illustrated in Figure 2.1 on the next page. 


THETA 

a 

nozzle attachment angle 

DEG. 

0.0 

THETAI 

sc 

nozzle inlet angle 

DEG. 

0.0 

THE 


nozzle exit angle (input 
if IWALL-4 only) 

DEG. 

0.0 

EPS 


nozzle expansion ratio (in- 
put if IWALL* 1 only) 

none 

0.0 

RMAX 


normalized radius at the 
nozzle exit plane (input if 

IWALL *2 or 3) 

none 

0.0 

ZMAX 


normalized axial position 
of the nozzle exit plane (in- 
put If IWALL*? or 3) 

none 

0.0 

RS(2) 

a 

table of normalized wall radii 
downstream of the nozzle 
tangency point. RS(1) is com- 
puted as r T * (Input if IWALL*4) 

none 

0.0 

ZS(2) 

m 

table of normalized axially 
positions downstream of the 
nozzle tangency point. ZS(1) 
is computed at Z. (input if 

IW ALL-4) r 

none 

0.0 
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isssa 


DescrlDtlon 

Units 

Assumed 
Value (s) 

NWS 

S 

number of entries in the RS, 

ZS tables. Includes the first 
entry. NWS * 20 . (Input if 
IWALL-4) 

none 

0 

RZN0RM 

S 

Optional normalizing factor 

none 

1.0 


for the RS, ZS table. For ex- 
ample, If RS & ZS were input 
as dimensional numbers, RZNORM 
would be the throat radius In those 
units, 

RI * normalized Inlet wall radius none 0,0 

of curvature. Only used If 
0DK calculation selected. 

It Is not necessary to Input 
RI even If an 0DK calcula- 
tion Is to be made. 


ECRAT 

m 

Contraction ratio at which 
the ballistics calculation 
will assume the nozzle In- 
let to start. 

none 

0 

FSUBM 

m 

length of submergence/length 
of Internal motor, (used In 
calculating the empirical sub- 
mergence loss for the nozzle). 

none 

0.0 

AEAS 

B 

nozzle entrance area/nozzle 
throat area (used In calcula- 
ting the empirical submergence 
loss for the nozzle). 

none 

1.0 

AVELS 

- 

motor average L* 

in. 

0.0 

KN0Z 

SEND 

C 

nozzle type flag 

-1 steel nozzle 

«0 all other types nozzle 

none 

0 
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of * 


2.4 Trajectory Directive and Input 


The TRAJECT0RY directive allows the user to input either the altitude or 
exhaust ambient pressure as a function of time during a motor firing In the $TRAJ 
namelist set. This data allows the SPP computer program to calculate the de- 
livered performance at the actual motor conditions. If this Input Is not speci- 
fied only vacuum performance can be calculated.' 

The format of the TRAJECT0RY directive and description of the Input is 
as follows: 

The word TRAJECT0RY must start in card column 1 . (Only the first 4 char- 
acters are checked). The card following the TRAJECTORY directive should contain 
$TRAJ In card columns 2 through 6. The cards following the $TRAJ card may be 
in any order upto the $END card which specifies the end of the trajectory Input 
data. All cards in the namelist input set must start in card column 2 or greater. 
The input items are 


Item 


Description 

Units 

Assumed 

VfllmeM 

$TRAJ 


namelist name 



PAM BI(1) 

m 

ambient pressure table 

PSIA 

0.0 

H(l) 

=5 

altitude table 

FT. 


TTRJU) 

m 

times corresponding to the 
ambient pressure or altitu- 
de tables. Must be mono- 
tonically increasing in time. 

sec. 

0.0 

NTJPT 


number of entries to be used 

In the Input tables. NTJPT* 1. 

If NTJPT* 1, then the following 
is done: H(2)-H(l); PAMBI(2)* 
PAM BI(1); TTRJ(2)*lxlO*. 

none 

0 

IPH 

m 

0 - pressure table Is to be 

used 

none 

0 


1 - altitude table is to be 
used 

SEND end of STRAJ namelist set 

Note: The PAMBI and H arrays have been equlvalenced to each other in order to 
save core storage. Hence, If both array names are Input, the last Input 
values will be used In accordance with the IPH flag. 
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2 . 5 Problem Directive and Input 

The PR0BIEM directive and $PR0B namelist input allows the user to 
specify which calculation modules are to be executed and which losses are to be 
considered. The occurence of the PR0BLEM card and associated data in the input 
stream causes the selected loss modules to be executed in the sequence shown 
in Table 2.1. 


The format of the PR0BLEM directive and description of the input is as 

follows: 

The word PROBLEM must start in card column 1. (All characters are 
checked). The card following the PROBLEM card should contain $PR0B in card 
columns 2 through 6. Ihe cards following the $PR0B card may be in any order 
up to the SEND card which specifies the end of the problem input data. All cards 
in the namelist input set must start in card columns 2 or greater. The input items 


are 

Item 


Description Units 

$PROB 


namelist name 

ODE 

SI 

r these input iters call out the cal- none 

CAL 


culation -nodules which are to ex- 
ecuted or specify ho w the dntu wMch 
would have generated by those , 

modules is tc be obtained . ihe 

JfDK 

■ < 

allowable values for these variables 
are: 

none 

i 0.0 either not con* leered or 

TDiP 


1 not executed, 

j 1.0 module la to bo executed. none 

THL 

' 

SB 

2.0 data previously generated by 

module ;s :o be road in or 
logical uni: IMP. 

none 



^3.0 values which would pe gen- 

erated uy modules are to be 
read in. 


iREO « 0 kinetic st* rt conditions for none 

0DK are to be read in via 
the SPECIES directive. 

1 kinetic start conditions tor 
0DK are to be computed 
from a 'restricted equilibrium 
ODE calculation', (see 
ECRAT in $ODE NAMELIST) 


Assumed 
Value (s) 


0 . 


0 . 


0 . 


0 . 


0 . 


0 


2-18 


Item. 


Description 


Units 


ETAId) 


Input values of the loss mechanism none 
efficiencies fractions. (Used only 
if module option *3, is selected) 


i. Mechanism 

1 finite rate kinetics 

2 2 dimenslonal-2 phase 
flow 

3 combustion efficiency 

4 submergence 

5 2 dimensional divergence 
loss 


INF 


Unit number on which previously none 

generated data Is to be read In on. 

(See Section 2.11) 

SEND 

s 

End of namelist set SPR0B 


Assumed 

Valuefa) 

1.0 
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2.6 


0DE INPUT DATA (ALL PROBLEMS SPECIFYING 0DE=1 .0 OR IREQ=1) 


The 0DE input data described here Is as defined In NASA SP-273, Refer- 
ence 1 # except namelists Input $INPT2 and $RK1NP have been combined Into a 
single list named $0DE. Only the rocket option (RKT In the namelist) is described 
here since all the other options have been removed from the 0DE calculation module. 

The 0DE input data consists of the following input groups: 

1. REACTANTS directive card, followed oy up to 15 data cards, 

followed by a blank card, specifying reactants. 

2. 0MIT and INSERT directives to omit or insert species for equilibrium/ 

frozen calculations. 

3. NAMELISTS directive card followed by input namelist $0DE 

specifying input case data. 

2.6.1 REACTANTS CARDS 

This set of cards is required for all 0DE problems. The first card in the 
set contains the word REACTANTS punched in card columns 1 to 9. The last card 
In the set Is blank. Tn between the first and last cards may be any number of cards 
up to a maximum of 15, one for each reactant species being considered. The cards 
for each reactant must give the chemical formula and the relative amount of the 
reectaat. For some problems, enthalpy values are required. The format and con- 
tents of the cards are summarised *n Table 2-G. A list of rome REACTANTS cards 
is given in Table 2-7. 

Relative amount s of reactants . - The relative amounts of reactants may be 
specified in several ways. They may be specified in toi r.s of molts, mole fraction, 
or mole percent (by keypunching M in card column d 3) or in terms of weight, weight 
fraction, or weight percent (blank in column >3). 

Relative amounts of total :\tvi to total oxidants can also oc input. Fnr this 
situation, each reactant must be specified as a ;uei or cn oxiciser by keypunch- 
ing an F or O, respectively, in col urn r. : 2 of the REACTANTS card. The amounts 
given on the REACTANTS cards are relouve o total fuel or total oxidant rather 
than total reactant. 


TABLE 2-6. REACTANTS CARDS 


Order 

Contents 

Format 

Card columns 

First 

REACiANTS 

3A4 

Hot 

Any 

One card for each reactant specie* 
(maximum 1$). Each card contains: 




(1) Atomic symbols and formula num- 
bers (maximum 5 sets) 1 

5(A2. F7. 5) 

Ito 43 


(2) Relative weight* or number of 
moles 

F7.S 

4t to 32 


(3) Blank If (2) is relative weight or 

M if (2) Is number of moles 

At 

33 


(4) Enthalpy or Internal energy 1 , 
rat mole 

Ft.3 

3410 02 


(S) r»Mt. N L. or C for solid, 
liquid or gas. respectively 

At 

II 


(6) Temperature associated with 
enthalpy in (4) 

FIS 

44 to 71 


(7) F if fuel or 0 if oxidant 

At 

72 


(8) Density in g cm 3 (optional) 

Ft. 3 

73 to 0) 

Last 

Blank 




a Procr:';w will calculate th* enthulpv or i«nermt energ;.* (4) fur sotcies la 


the THERMO data at the temperature (6) if .'.cro« ire punched in card 
columns 3? and 38. (See section Reactant enthalpy fur additional in* 
formation.) 

^Relative weight of fuel in total fuels or oxidant in tola! oxidants Alt 
reactants must be given either all in relative wem hts or a!! in number 
<4 moles. 


L 
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TABLE 2-7 LISTING CF SAMPLE REACTANTS CARDS 


4 


REACTANTS 



H 2. 

100. 

0. 6298,15 F 

N .7806810 .209795AR, 00*662 

100. 

•7.2021646298,15 0 


REACTANTS 

N i, H 

4, CLl. 0 4. 

72,06 

•70730, 

S298,15 

F 

Cl, H 

1.669550 .031256S ,008415 

18.56 

-2999,0821296,15 

F 

ALl. 


9,00 

♦c.o 

8296.15 

F 

*6i. 0 

1. 

.20 

•143700. S296.15 

F 

N 2. 0 

1. 

.16 

•66317.4 L296.15 

F 


reactants 
« 2, 

00 

100. 

0. 

6296,15 F 

0 2. 

00 

100.0 

0.0 

6298,15 0 


REACTANTS 

N 2, Hi, C 2, 

50,0 

12734,8 

1296,15 

F ,766 

N ?, M 4. 

50,0 

12050. 

12*8,16 

F 1.003 

F 2, 

100, 

•3050.8921 85,2* 

0 1,5* 

• 

reactants 

III, 

• 

100. 

C. 

$296,15 

r 

F 2. 

100, 

•3030. 892L 65,2* 

0 1,54 


REACTANTS 

N 2, M 4, 

60. 

12100, 

L298.ll 

F 

1.003 

6E1. 

20. 

0.0 

8296.15 

F 

1.65 

M2, 02. 

100. 

•44860. 

L296.15 

0 

1.407 


"Luted above ari six examples. Each example must end with a blank card. 
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There are four options in the $0DE namelist for indicating relative 
amounts of total fue.! to total oxidant as follows: 

1. Oxidant to fuel weight ratio (0F is true). 

2. Equivalence ratio (ERATI0 is true). 

3. Fuel percent by weight (FPCT is true). 

4. Fuel to air or fuel to oxidant weight ratio (FA is true). 

For each option, except 0DE with 0DK^1, and TD2P^1, the values are given in 
the 0FSKED array of $0DE (described in Section 2.6.4), For 0DE=1, 0D1^1, and 
TD2P^1, the MIX array is used as described in Reference 1. 

Reactant enthalpy . Assigned values for the total reactant are calculated 
automatically by the program from the enthalpies of the individual reactants. 
Values for the individual reactants are either keypunched on the REACTANTS cards 
or calculated from the THERM0 data as follows: 


Enthalpies are taken from the REACTANTS cards unless zeros are punched 
in card columns 37 and 38. For each REACTANTS card with the "00” code, an en- 
thalpy will be calculated for the species from the THERM0 data for the tempera- 
ture given in card columns 64 to 71. 


When the program is calculating the individual reactant enthalpy for values 
from the THERM0 data, the following two conditions are required: 

1. The reactant must also be one of the species in the set of THERM0 data. For 
example, NH^g) is in the set of THEPM0 data but NH^(d is not. Therefore, 
if NH,(g) is used os a reactant its enthalpy could be calculated automatically, 
but that of NH^(: ) could not be. 


2. The temperature T must be in the range T Jow , 


/ X • 4 ^ l 


to T^ ^ Is the temperature range of the THERM0 data. 


high 


a i.i) who*** 


low 


For many cises it may be desirable to modify the enthalpy entered on the 
REACTANT cards. Th is car* be done by using the DELH input array. The DELH entry 
will be added to the system enthalpy as computed by 0DE from the reactants cards 
(see above). For example, overall system enthalpy of the propellant can be input 
through the reactants cards and the energy added or extracted due to conditioning 
the propellant can be input by the DELH entry. An alternate method would be to 
input zero enthalpy un the Reactants cards and input enthalpy by the DELH entry. 
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The second and less important reason is that if one knows that one or 
several particular condensed species will be present among the final equilibrium 
compositions for the first assigned point, then a small amount of computer 
time can be saved by using an INSERT card* Those condensed species whose 
chemical formulas are included on an INSERT card will be considered by the 
program during the initial iterations for the first assigned point. If the INSERT 
card were not. used, only gaseous species would be considered during the initial 
iterations* Howevar, after convergence, the program would automatically insert 
the appropriate condensed species and reconverge. Therefore, it usually is 
immaterial whether or not INSERT cards are used. For all other assigned points 
the inclusion of condensed species is handled automatically by the program. 


NOTE: For all metalized propellants It Is strongly suggested that the appropriate 
condensed phases be included in the insert cards to Insure convergence 
in the equilibrium calculation. 




2 . 6.2 


0MIT and INSERT Cards 


0MIT and INSERT cards are optional. They contain the names of particular 
species in the library of Thermodynamic Data for the specific purposes discussed 
below. Each card contains the word 0MIT (in card columns 1-4) or INSERT (in 
card columns 1-6) and the names of from 1 to 4 species starting in columns 16, 

31, 46, and 61. The names must be exactly the same as they appear in the 
THERM0 data. 

2. 6.2.1 0MIT Cards 

These Coids list species to be omitted from the THERM0 data. If 0MIT 
cards are not used, the program will consider as possible species all those 
species in the THERM0 data which are consistent with the chemical system 
being considered. Occasionally it may be desired to specifically omit one or 
more species from considerations as possible species. This may be accomplished 
by means of 0MIT cards. 

2. 6.2. 2 INSERT Cards 

These cards contain the names of condensed species only. They have 
been included as options for two reasons. 

The first and more important reason for including the INSERT card option is 
that, in rare instances, it is impossible to obtain convergence for assigned enthalpy 
problems (HP or RICH without the use of an INSERT card. This occurs when, by 
considering gases only, the temperature becomes extremely low. In these cases, 
the use of an INSERT card containing the name of the required condensed species 
can eliminate this hind of convergence difficulty. When this difficulty occurs, 
the following message is printed by the program: "LOW TEMPERATURE IMPLIES 
CONDENSED SPECIES SHOULD HAVE BEEN INCLUDED ON AN INSERT CARD". 
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2.6.3 $0DE NAMELIST INPUT 


The 0DE subprogram contains namelist input sections $0DE. The Name- 
list $0DE must be preceded by a card with NAMELISTS punched in card columns 
1-9 • 


The $0DE Namelist is required for inputs to the $PR0B namelist set if: 

0DE = 1.0, 

0DK = 1.0, 

or 

IREQ = 1 . 

If both 0DE=1, and IREQ=1, are specified in the $PR0B namelist, two sets of 0DE 
data will be require d. 

The variables input by the $0DE namelist are listed in Table 2-8. Addi- 
tional information about some of these variables follows: 

Pressure units . - The program assumes the pressure in the P schedule to be in 
units of atmospheres unless PSIA - .TRUE, 

Relative amounts of fuel(s) and oxidizer (s) . - These quantities may be specified 
by assigning 1 to 15 vlaues for either o/f, %F , f/a, or r. If no value is assigned 
for any of these, the program assumes the relative amounts of fuel(s) and oxidlzer(s) 
to be those specified on the REACTANTS cards. (See discussion in REACTANTS 
Cards, Section 2.6.1). 

RKT problem . - Only one value for chamber pressure, P, is to be input. 

Print out will be given for the chamber pressure condition (l.e. stagna- 
tion) and the throat condition. Print out may be requested at other conditions by 
use of the PCP schedule and the SUBAR and SUPAR schedules. The values of AS UP 
and ASUB from the $GEOM namelist set will automatically be Inserted into the SUPAR 
and SUBAR print arrays. Hi us any input into these arrays will override the values 
input ir the $GE0M namelist set for the 0DE calculation only. 

The prograrr will calculate both equilibrium and frozen performance un- 
less FR0Z=F or EQL=F are input. If FR0Z=F, only equilibrium performance will 
be calculated. If EQL-F, only frozen performance will be calculated. 
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TABLE 2-8. VARIABLES IN $0DE NAMELIST 


Variable 

NO * Of 

entries 

Type 

Value 

before 

read 

Definition and comments 

RKT 

1 

L 

True 

Rocket problem 3 

p(d 

25 

R 

0 

Assigned pressures: stagnation 
pressures for rocket problems: 
values in atm unless PSIA, or SI= 

. T. , (see below) 

SI 

i 

■ 

False 

a Values in P array are in N/m a 

PSIA 

1 


False 

3 Values in P array are in psia units 

xpc) 

50 

R 

1. 

Not currently used 

0 F 

1 

L 

False 

Oxidant to fuel weight ratios are 
to be input a 

ERATI0 

1 

L 

False 

Equivalence ratios are to be input 3 

FPCT 

f 

i 

L 

False 

Percent fuel by weight are to be 
input 3 

FA 

I 

L 

False 

Fuel to air weight ratios are to be 
input 3 

0FSKED&) 

50 

R 

0 

For a RocI et problem, and ODK~ 

1 or TD2P=1 , 0FSKED will be used 
rather than MIX (see Reference 1 ). 
Relative amounts of total oxidant 
to total fuel are input as defined 
by OF, ERATI0, FPCT, or FA. 

DELH(I) 

5C 

R 

0 

This value will be added to the 
system enthalpy input thru the re- 
actants cards. Units are BTU/ 

# if PSIA- .T. , otherwise cal/gram. 

I0NS 

1 

L 

False 

Consider ionic species 3 

PCP(I) 

50 

R 

0 

Compute and print solutions at these 
ratios of chamber pressure to pres- 
sure (entries must be < 1.) 

SUBAR(I) 

50 

R 

0 

Compute and print solutions at these 
values of subsonic area ratios (en- 
tries must ^ 1 .) 

SUPARC) 

50 

R 

0 

Compute and print solutions at these 
values of supersonic area ratio (en- 
treis must ^1.) 


a If variable Is set to be true. 

b Set variable false U these calculottons are not desired. 
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Table 2-8 (cont'd) 


Variable 

No. of 
entries 

Type 

Value 

before 

read 

Definition and comments 

ECRAT 

1 

R 

0 

Subsonic area ratio to start 0 DK 
calculations with computed equil- 
ibrium conditions. The SUBAR in- 
put table must inplude an entry 
equal to ECRAT. 

EQL 

1 

L 

True 

Calculate rocket performance as- 
suming equilibrium composition 
during expansion* 5 . 

FR0Z 

1 

L 

True d 

Calculate rocket performance as- 
suming frozen composition during 
expansion* 5 . 

LISTS? 

1 

L 

False 

List names and dates of all species 
residing on thermodynamic data 
used 3 . 

KASE 

1 

I 

0 

Optional assigned number associated 
with case. 

-0LPT 

1 

L 

TRUE 

For cases with condensed phases, 
the equilibrium conditions at the 
onset of solidification will be cal- 
culated. 

NLPRNT 

1 

L 

FALSE 

Write out the $0DE namelist data 
set using a namelist write. 3 


°If variable is set to be true. 


b 


c 






b 


Set variable false if ‘here calculations are not desired 
Required only for the restricted equilibrium option, l.e 
Assumed falcc for a restricted equilibrium calculation. 


• » 


IREQ=1 in $PR0B namelist. 
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2 7 Grain Design and Ballistics Input. 

The input for this module Is described in three parts. Section 2.7 1 de- 
scribes the $BAL NAMELIST input, while Sections 2.7.2 and 2.7.3 escr e 
formated grain design and formated ballistics input, respectively. 

2.7.1 $BAL NAMELIST 

In the following list of input quantities an * preceding the variable name 
means that it is internally transmitted from the 

$PR0B NAMELIST 'at ^ * are used to denote variables that 

. *- ~ ~ 

will over-ride internally transmitted values. 


Variable 

♦PC 

♦CSTART 

♦TFT 

♦CPG 

♦CPM 

♦GAM 

♦XMG 

♦XMM 

*FCP 

♦XMV 

♦XKG 

*XB 

**RN(1) 

**RN(2) 

**RN(3) 

DTIME 

KBAL 

MODE 

**NNOZ 


Description 

Estimate of initial chamber pressure 
Theoretical O* 

Theoretical fleme temperature 
Specific heat of the gas only 
Specific heat of combustion products 
Ratio of specific heats 
Moles of gas per lOOgm of products 
Moles of products per lOOgm of products 
Moles of condensed phase per 100 gm of 
products 
Gas viscosity 
Gas therm :.l conductivity 
Nozzle erosion parameter, * 

Nozzle entrance radius 
Initial nozzle throat radius 
Nozzle exit radius 
Not used 

Logical unit number for internal informa- 
tion transfer (set equal to 8 in MC module) 

Integer used for control of Internal inform*- 
tlce transfer. (Set equal to 0 in MC module) 

Number of nozzles 


Units 

lb/in 2 

ft/sec 

°K 

Btu/mol/°R 

Btu/mol/°R 

mol/ lOOgm 
mol/ lOOgm 
mol/lOOgm 

lb/in/ see 
Btu/in/secAR 

in 

in 

In 
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2.7.2 Grain Geometry Inputs 



The data cards required to define the initial grain geometry are divided in- 
to two categories: (1) Initial data cards and (2) Figure cards. This input must be 
preceded by a title card. This title card is not used, however. The title card de- 
scribed in Section 2.1 is used instead. 

2.7.2. 1 Initial Data Cards 

INITIAL is punched on the cards starting in column 1. These cards con- 
tain data which define the region in a coordinate system where the surfaces will 
be located plus data which control the limitations and X, Y and bum increments. 

The burn interval may be divided into a maximum of four sub- intervals; e.g., BO, 

Bl, B2, B3, B4, where the Kth sub-interval (the interval defined by B^) 

would be divided into NB^ parts. This may be desired where, for example, it is 
obvious that a complicated shape will bum into a simple shape, warranting a 
change lu rr.esh size to save computer time. 

The subscripting for X, NX and Y, NY is also identical* and is related to 
the burn intervals cj follows: for each burn sub- Interval, a complete set of X and 
Y values must be p r epared, thus the reason for the double subscript on these 
quantities. The procedure for the first subscript is the same as for the bum and 
burn increments and the second subscript corresponds to the burn intervdl. To 
illustrate, suppose the burn interval is divided into two sub-intervals BO to Bl to 
Bl with NB1 and NB2 parts; thus, two sets each of X and Y data are required, as 
follows: 

InteA'al BO tc Bl: 

X01, XI 1, NX1 1 , X21, NX 2 l , X31, NX31 uptoX41 f NX41* 

YOl, YU, NYU, Y21, . ., Y41, NYU 
Interval Bl to B2: 

X02, x:2, NXi:, X 22 , NX?2 X12, NX4 2 

Y02, VI 2, NY12 Y42, NY.2 

♦EXCEPTION: When the only burn interval is from BO to Bl it is possible to have :p 
to sixteen changes in the X interval. The subscripting -hi the X input would bo as 
follows: 


From one to nina interval changes: from ten to sixteen changes (continuing) 



X01 XI 1 NXI 1 
X21 NX21 
X31 NX31 
X41 NX4i 
X51 NX51 


X01 NX61 
X71 NX71 
X81 NX81 
X91 NX91 


X03 NX£3 X43 NX63 

X13 NX33 X53 NX73 

X23 NX43 X63 NX83 

X33 NX53 
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where X01 » beginning X 
XI 1 = ending X 

NX11 = number of Intervals in between X01 and XI 1 
X21 = next ending X 

NX21 * number of intervals in between XI 1 and X21 
etc. 

If too large a mesh size is specified a plot of surface area, or pressure, 
versus bum distance, or time, will reveal a saw tooth effect. 



To avoid this situation the maximum X and Y Increments should be based on the 
following criteria: the maximum X and Y increments, in order that a step or saw 
tooth effect in the plot of surface area or pressure versus bum distance or time 
is not obtained, are: 


max AX * aB esc A 
max * /\B sec A 

where A is the angle between the X axis and the line tangent to the point of in- 
tersection of the actual void and the representation figure that constitutes most 
to the surface area in the region. 

Y 



2-31 



It is not necessary for the number of increments of X and Y corresponding 
to the first bum interval to equal the number of increments of X and Y correspond- 
ing to the second burn interval and also for any more burn intervals. The restric- 
tions on X and Y are that the initial and final values for X are equal and the initial 
and final values for Y are equal. To illustrate, suppose the following case: 

For the Interval BO to Bl: 

X01 = 0.0 
Xll = 5.0 
NX11 = 10.0 
X21 = 10.0 
NX21 =6.0 

Tor the interval Bl to B2: 

X02 = 0.0 
X12 = 10.0 
NX12 = 10.0 

In this case X01 = X02 and X21 = X12. The same idea follows with Y. 

The number of bum intervals must be a whole number, and the number of X and 
Y increments must be whole numbers, and the X and Y increments must be even whole 
numbers. 

The burn interval, or sum of burn sub- intervals, constitutes the total pro- 
pellant web. The number of bums, NB, determines the frequency of geometry out- 
put; thus, NB 20 will yield answers for each 5 of web burned. This Is not to bo 
confused vath NX or NY, wHch define the caiculatlonal rresh at a given time or 
web burned, or with separate conuol of ballistics output. 

The computational time will be proportional to the product of the number of 
X un i Y increments, in other words, if both «re doubled, the time will increase by 
roughl y a factor of A . 

2.7. 2.2 Figure Cards 

The first card for each surface contains the names of the figure starting ir. 
column 1. The cards immediately following contain the data defining the figure 
•md how it will burn. The first card contains Information of whether the fiqure is 
to consist of propellant or void. If the figure is to consist of propellant, a CR is 
punched in columns *9-20 and the figure is considered to burn Inward; otherwise 
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the figure Is considered void and bums outward. 

The burning of the figure Is controlled on a BURN card. A figure may burn 
In any combination of the following types of burning. (This card is used only if 
a figure is to be burned in a special manner.) 


Type 

Bum 

Variables 

Name 

Columns 

Value 

Columns 

Bum 

Considered 

Action 

Nonbum ing 

NB 


Bum = 0.0 

No burning occurs in 
the figure. 

Delayed bum 

DB 

a 

Bum = B-a 

The figure is delayed 

In burning the speci- 
fied amount (a). 

Rated burn 

RB 

b 

Bum - (B)(b) 

The figure is burned 
b times as fast (or 
slow) as the normal 
bum. 

Corner rounding 

CR 

c 


The figure (prism, 
cone or cylinder) Is 
backed-up the distance 
c and then burned out 
again the distance. 


The data defining the figure is contained on the DATA cards. Punched in 
columns 1-4 of these cards is the word DATA. These cards will contain the control 
characters XI, X2, X3, Yl, Y2, Y3, Zl, Z2, Z3, R, Rl, R2, H. The X, Y, and Z 
characters define points in an X, Y and Z coordinate system. R is the radius and H 
is for prism height. In the case of a cone, where two radii are required, Rl and 
R2 are used. When R1-R2; e.g., in a cylinder, the R is used, which will cause 
the same value to be stored in Rl and R2, 


The figures the program will handle and their inputs are: 
Cone (right) Two ooints (the center of two bases) 

Two radii (the radius of the two bases) 
(Rl is radius at point 1 and R2 is radius 
at point 2) 

Cylinder Two points (the center of the two bases) 

Ught) 

One radius (the radius of the cylinder) 


XI. 

Yl. 

Zl 

X2, 

Y2, 

Z2 

Rl. 

R2 


XI. 

Yl. 

Zl 

X2, 

Y2, 

Z2 


R 
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Prism (right 
triangular) 


Throe points (corners of bisecting 
plane) 


XI, Yl, Zl 
X2, Y2, Z2 
X3, Y3, Z3 

Height of prism H 

Sphere One point (center of sphere) XI, Yl, Zl 

Radius of sphere R 

2.7.2. 3 Propellant Periphery and Exposed Insulation 

The complete description of the motor geometry requires a series of cards 
in the form of a table of outer propellant radius versus axial position. The X con- 
vention should be consistent with that describing the configurations above, but 
the number of entries need not be consistent with NX. Only the extremities in X 
need be consistent with those on the initial *INIT) cards. Up to 100 entries are 
allowed; however, far fewer will generally suffice. 

Insulation area exposed as the web burns is not computed, but may be in- 
put in the form of a table. (See Section 2,7.3, T-ble 2-10) 

2.7.3 Ballistics Inputs 
2.7,3, 1 Input Scheme 

Following the last card containing groin design Information, but preceding 
the ballistics input, a card containing the word COMPUTE, starting in card column 
1 , must appear. 

The ballistics input is controlled from subroutine TABIN, but is actually 
read in subroutine DA TAIN, in the following sequence: 

1 


2 


3 
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All of the data tables (1) are read first, with a limit of 50 tables and 3000^ words 
of storage, both about three times more than normal input for one case. The tab- 
>r input is assumed to continue until a READ IN completion code (loc. 4*-l) is 
given. The input in the remaining two sets (2) and (3) is then Interpreted in the 
program. Integer words (2) are logic controls for the calculation, and data words 
(3) include coefficients, constants, and some problem logic control words, to be 
described in detail in subsection 2.7.3 .2. 

The input is transferred to working storage in the main program (4) if non- 
zero (nz) values are given. As a result, only changes in the problem are needed 
in successive cases. Also, only input pertinent to the case and options to be 
used is required, so unneeded input can be completely ignored. Since zero input 
is ignored, no "filling N is required, but values in the cards do no harm if desired 
for column alignmert, etc. One limitation in the method is that options can not 
be "turned off" by inputting a zero in a later card. The node points are derived from 
the NX values in the grain design input. Points defining radial slots Include one 
point at the end of each grain segment and one (only one) in the slot. The number 
of NX Increments is given as -2.0 for the special case of a radial slot to key the 
logic. Once accomplished, the number is made positive and the grain design 
program operates as u3ual. 

2. 7.3.2 Details of Input 

In this subsection, the mechanics o. r card format and table input are con- 
sidered; then the parameters and requirements are specified. 

D&ossi 

Three card Input formats are used by the program, 
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4 (12. 14, 6#, 12, 28A 

Logic, first of set 

6 (12, 14, fiX), 12 

Logic, remaining cards in set 

6 (12, F10), 12 

Data 


For either logic (Integer) words or data cards, the Input is Interpreted as (loc for 
location) . 

1-2 3-12 13-14 

loc, Input loc — 

value 

12, 14, 6XorF10I2 

The Input cards are read In sets, and a set Is assumed to continue until 
a completion sign is encountered . Two completion signs are recognized: 

(1) A blank card 

(It) A negative location, l.e. , location « -1 In any of the 12 fields 
on the card 

Note that the completion sign for the group of tabular inputs is different In mean- 
ing and form { ^4^ ^ in 12, 14), 

Rules governing the interpretation of a card set include: 

Any number of cards may form a set 
Any order of locations is acceptable 
Zero values of input are ignored 

If locations are repeated, the last value of input is used (overwrites) 

If a location is blank but data is given, a sequenced list is assumed, 
so the following two cards arc equivalent: 

1 10. 2 11. 3 12. -1 

1 10 . 11 . 12 . *1 

The sequencing doer not follow from card to card, so the first location on the 

next card must be given* 

Tabular Data 

Data tables are assumed of the form z»f (x, y). Tables are termed rec- 
tangular If the number of x arguments in the rows is the same for each of the col- 
umns of y arguments. A non -rectangular table, then, has a varying number of x 
arguments in the columns. 

The card sets to input a table are given in Table 2-9 with the location 
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Table 2-9 Input of Tabular Data 


Input for Form Z = f (X, Y) 


Input Type 

Input 

Description 

Integer card set 

Required 

Control constants 

Data card set 

Required If NY > 1 Omit If NY=1 

Y values 

Data card sets 

1 set If rectangular NY sets If 
nonrecta ngular 

X values 

Data card sets 

NY sets required 

Z values 


Integer Cards 4 (12, 14, 6X), 12, 28A; 6(12, 14, 6X) 12 


Location 

Input 

Description 

1 

Required 

Table number (see Table 2- 1 0) 

3 

Optional 

Key to LOG-LOG interpolation (+) 

4 

Optional 

Table type non -recta ngular = (+) (Note that (-) 

Is key to end table inputs) 

5 

Optional 

Number of Y (columns) values, required if NY > 1 

6 

Optional 

Number of X (rows) values, required if NX > 1 

7 to 25 

Optional 

Number of X values in 2nd 19th Y columns, 

non-rectangular tables 


Data Card Sets 6 (12, F10), 12 


Location 

Input 

Description 

1 

Any 

Optional 

Required 

List of X, Y, or Z values sequenced from Location 1 
Completion sign 


tyOTES : 1. A table can degenerate to one point In X, Y, or both. 


2. Interpolation Is linear with no extrapolation: values at boundary 
of table are returned. 

3. A missing table causes Z=0 and one (only one) error message. 
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of the words in the sets. The logic (integer) card set is always required and must 
include a table number and a completion sign. The number of y columns (NY) in 
location 5 and x rows (NX) in location 6 must be given if more than one. Location 
3 is used to key log-log interpolation if desired, and a (+) value in location 4 keys 
a nonrectangular table. 

Usually, one set of y values, one set of x values, and NY sets of z values 
compose a table. If only one column is specified, the set of y values is omitted. 

It a non-rectangular table is specified, the number of x values in each column 
must be given, starting in location 7, and an x argument set is required for each 
y argument. 

Assignment of Input Tables 

The numbers of input data tables are assigned in the subroutine MAIN3 due 
to compiled-in constants in the call statement: 

CALL TAB (X. Y, Z, NO) 

The correspondence ox Input data and the table numbers, NO is given in Table 
2-10 with the FORTRAN symbol and assumed arguments of the tables. In Table 
2-10, data that is designated as "optional" is needed only if the related option 
is specified. Data that is designated as required is almost always needed since 
zero values of the parameter (Z) are rarely meaningful; the possible exception is 
the list of times of detailed output where no input results in printout at all time 
steps. 

Integer Inputs 

A series of integer words may be input to control program options, as shown 
in Table 2-il. Loc (4)^NUM(4)--1 is the required input to indicate the end of data 
tables. The none number of each radial slot is required input. Other input In Table 
2-11 is used only if the particular option is desired. 


Data Input Words 

Up to 99 data (real) word locations are provided to complete the input 
constants and control options in the program, as listed in Table 2-12. Note 
that one data set is used, regardless of the number of cards, as discussed earlier. 
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Table 2-10 Input Table Assignments 


HI 

Function 

Input 

Z Symbol 

X_ 

Y 

1 

Propellant strand burn- 
ing rate at conditioning 
temp. 

Required 

RO 

PC) 

- 

2 

Pressure exponent 

Required 

SLOPE 

PC) 

- 

3 

Motor geometry (grain 
outer radius) 

Required 

R(l) 

X123 

- 

4 

Insulation area 

Optional 

AI 

TIME 

- 

5 

Time increments 

Required 

DT 

TIME 

- 

6 

Time of detailed output 

Required 

TLIST 

TIME 

- 

7 

Throat area table 

Optional 

AP (NN); 
KODE=5 

TIME 

- 

8 

Slot bum rate correction 

Optional 

DR 

WEB (I) 

X(I) 

13 

Slot diameter 

Optional 

DSLOT 

TIME 

AFI(L) 

14 

Slot bum-back correction 

Optional 

DX 

WEB (I) 

AFKL+10) 


Table 2-11 Inteaer Input Parameters 


Location 

Symbol 

Input 

Description 

Num (4) 

- 

Required 

(=-l) indicator of end of tables 

Nam (5) 

Ign 

Optional 

Key (+) for ignition dt option 

Num (11-20) 


Optional 

Node number of radial slots 

Num (21-30) 


Optional 

Propellant number for following 
segment 
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Table 2-12 Data Input Parameters 


Location 


Units 

Input 

Description 

Z(l) 

TSTOP 

sec. 

Required 

Calculation run (expected bum) 
time 

2(2) 

TRISE 

sec. 

Optional 

Time of rapid pressure rise 
(Ignition rise), an aid to con- 
vergence 

Z(3) 

com 

— 

Required 

Convergence ratio (.00005-. 005) 

?.<:>) 

PAM 3 

lb/in 8 

Required 

Ambient pressure 

2(7) 

KODE 


Required 

Nozzle material 

1. carbon phenolic 

2. ATJ Graphite 

3. pyrolytic graphite 

4. non-eroding 

5. a table of erosion rate Is to 
be Input 

2(8) 

KODL 1 

— 

Required 

1. If fully Insulated motor 

2. If uninsulated 

Z(9) 

KODE 2 


Required 

Binder type 

1 . PBA Binder 

2 . HTPB Binder 

3 . CTPB Binder 

4. NC Binder 

5. PU Binder 

6. PGA/NC Binder 

7. NF Binder 

Z(10) 

ALUM 

per cent 

i 

Requires 

Aluminum concentration 

Z(ll) 

MMO 

gm/mole 

Required 

Molecular wt. metal oxide 

2(12) 

TAMB 

°K 

Required 

Ambient temperature 

2(13) 

FACTOR 

-- 

Required 

Loop gain ratio 

2(14) 

PSTART 

lb/ln ? 

Required 

Trial value of motor pressure 

Z(1 5) 

ETAF 

fractional 

Required 

Nozzle efficiency; 1.000 in 
the Integrated mode 

2(16) 

RHO 

lb/ln 3 

Required 

Propellant density 

2(17) 

RHOI 

lb/ln 3 

Optional 

Insulation density 

Z(lfl) 

RI 

in/sec 

Optional 

Insulation ablation rate 

Z(19) 

ANOZ 

In 2 

Optional 

Projected nozzle surface area 
for reradlatlon to propellant 

Z(2U) 

ETAIN 

fractional 

Required 

Recovery efficiency for radial 
slots 
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Talkie 2-12 Data Input Parameters (Continued) 


Location 

Symbol 

Units 

Input 

Description 

Z(21) 

Cl 

in/sec 

Required 

Erosive burning constant 

Z{22) 

C2 

— 

Required 

Erosive burning exponent 

Z(23) 

RHOAMB 

lb/ln® 

Required 

Ambient gas density 

Z(24) 

CS 

BTU/lb°R 

Required 

Heat capacity of propellant 

Z(2S) 

TSUP 

°K 

Required 

Propellant surface temperature 
(850 for AP-composite, 520 for 
double-base) 

Z(26) 

ALPHA 

ln a /sec 

Required 

Propellant thermal diffuslvity 

Z(41)-Z(50) 

AF 


Optional 

Key word for radial slot ig- 
nition, 10 forward faces, 10 
aft faces (value = ignite) 

Z(71)-Z(80) 

A FI 

in 

Optional 

Table argument value for up to 

10 radial slot ignition diameter 

Inputs (key to input table number 

13). Ignited diameter becomes re- 
quired Input if slot face is burning. 

Z(91)-Z(100) 



Optional 

Key word to suppress pressure 
rise down the port. Requires 
definition of slots and segments 
to function. 


Termination of Ballistics Data Set 

Following the last card containing ballistics input a card containing the 
words END 0F PR0GRAM, beginning in card column 1, must be inserted. 










2.8 0DK INPUT DAT A 


0DK input data is required only in 0DK = 1 as specified in the $PR0B 
namelist set. 

2.8.1 SPECIES 

Species used by the computer program are determined in several possi- 
ble ways, depending upon the problem type. Methods used to determine 
chemical species for each problem type are discussed below. 

0D K f, 1 , IREO " 0 

For 0DK problems, species names and concentrations must be input, 
see Section 2.8.1 . 1 

0DK - 1 , IREQ = 1 

For 0DE (restricted) - 0DK problems, the initial start conditions for 
the kinetic expansion are obtained from a restricted equilibrium calculation. 

The species list generated by the equilibrium calculation generally contains 
many more species than the 40 species for which the <2fc)K subprogram is dimen- 
sioned. Therefore a selection process is required to interface the 0DE calcu- 
lated equilibrium start conditions with the GfDK kinetic expansion calculations.' 
This selection is performed using the following rules: 

Kuie 1 If a species appears in a reaction, it is selected for the 
kinetic calculation. 

Rule 2 If a species is specified using INERTS directive it is 
selected for the kinetic calculations. 

Rule 3 If any species has a mole fraction greater than an input 
criterion, it is selected for the kinetic calculation. 

Species which arc selected but which do not appear in a reaction are treated 
as inert and listed as such on the output list of selected species . 
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2. 8.1.1 (*DK OPTION FOR INPUT OF INITIAL SPECIES CONCENTRATIONS 
(APPLIES TO THE 0DK ONLY PROBLEM) 

This input begins with a single card with SPECIES in columns 1 through 
7 and with either MASS FRACTIONS or M<2<LE FRACTIONS in columns 9 through 22. 

If the identifier for mass or mole fractions is omitted, mass fractions are assumed. 
Up to 40 species cards may be input. Only those species specified by input 
species cards will be considered for an <2fc)K problem . The order of the input 
species cards is independent of the order in which the species appear on the 
master Thermodynamic Data file. 

A chemical species is Identified symbolically by 12 alphanumeric char- 
acters and must correspond identically with the species name as it appears on 
the Thermodynamic Data file. A complete list of the current species names are 
listed in Table 2-4. The species symbol may not contain the characters * or =. 

Col Function 

1-10 Not used 

11-22 Species symbol (left -Justified) 

23-30 Not used 

31-60 Value of initial species concentration (if zero 

must be input as 0.0) free field F or E format 

61-80 User Identification if desired 

2.8.2 REACTIONS 

Chemical reactions must be input for an 0DK problem. 

Up to 50 reactions with an implied third body and a total 150 reactions 
may be input to the program. Only one card per reaction, and only one reac- 
tion per card is permitted. Cards specifying third body reactions must precede 
cards specifying all other reactions. Species names appearing in the symbolic 
reaction set mur correspond identically with the species names as they appear 
in the master Thermodynamic Data (see Table 2-4). A card listing for a sample 
reaction set is presented in Table 2-13. 
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Table 2-13 Sample Aluminized Propellant Reaction Set 
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The symbolic reaction set contains directive cards and reaction/data 
cards in groups as outlined below: 

REACTIONS Directive for start of symbolic reaction 

input 

# Reactions with Implied third body species 

END TBR REAX Directive for end of third body reactions 

• 

# All other reactions 


LAST REAX Directive for end of reactions 

INERTS Specified inert species 

THIRD BGtoY REAX RATE RATI0S Directive for start of third body reaction 

rate ratios 

• Third body reaction rate ratios 

LAST CARD Directive for end of REACTIONS input 

The content and format of each type of card is defined as follows: 

2.8.2. 1 The symbolic reaction set begins with a ^ard containing the word 
REACTIONS in columns 1 through 9, Other columns on this card can be used 
for comments . 

2 .3 .2. 2 Each card defining a reaction is divided into five fields , separated 
by commas. Each field contains: 

field 1 a reaction 

field 2 A - followed by a value of A 

field 3 N - followed by the value of N 

field 4 B - followed by the value of B« the 

activation energy (Kcal/mole) 

field 5 available for comments 

The general form of a reaction is: 

Nj * Symbol j + N 2 * Symbol 2 + . . . « N ft *Symbol a + N fe *Symbol b 4 . . . 

where the left hand side represents reactants and the right hand side represents 
products. The reaction can be either endothermic or exothermic. 


rate parameters 
for the reactions 
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The multipliers, N, must be integers and represent stoichiometric 
coefficients. If no stoichiometric coefficient is given, the value 1 is assumed. 
The dimensioning currently in the program requires that: 

N + N 2 + ... <10 

and 


N + N. + . • • <10 

a o 

The chemical species (denoted by the word " symbol *' in the above gene- 
ral form) can contain up to 12 characters, each of which must match a species 
name contained in the thermodynamic data (see Table 2-2, card 3). Examples: 

Reaction Interpretation 

NA++CL-=flACL Na* + Cl’ « NaCl 

R+2-H?-2=B0 B ++ + (f~ = B0 

BE+2f2*<2(H-=BE0H0H Be ++ + 20H* = Be«2)H) 2 

The value assigned to A, N, B defines the forward (i.e. left to right) reaction 
rate, k, as 

k = A • T* N • e -< l000B ^ 
in units of cc, °K, mole, sec. 

All three reaction rate parameters must be input. The numeric value of each 
parameter may be specified in either I, F, or E format. If E format is used, 
the E must appear before the exponent. s 

2.8.2. 3 The reacaons with an implied third body must precede other types of 
reactions, and mujt be followed by the directive (columns 1 through 12): 

Col 1 | 

END TBR REAX 

all reactions prior to the above directive will have a third body term added to 
each side of the reaction, e.g.. 
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H2 = H + H, ... 

END TBR REAX 

defines the chemical reaction 
h 2 + m =H + H + M 

where M is a generalized third body. Specific third body effects may be 
included by lnputing specific third body reaction rate ratios as outlined in 
2.8.4. Cards encountered after the END TBR REAX directive card do not have 
a third body term added. 

2.8.2. 4 All other reactions are input next, exactly as. described under 2. 8. 2. 2. 

2.8.2 .5 After the last reaction has been defined, a card with LAST REAX in 
columns 1 through 9 is input. 

2.8.2 .6 Reaction rate data for 14 dissociation-recombination (implied third 
body) reactions and 59 binary exchange reactions are listed in Table 2-13 for 
propulsion systems containing elements C, N, O, H, Cl and Al. Cards can 
be abstracted from Table 2-13 for input to the computer program. For the 
implied third body reactions, the third bcdy for which the rate applies is indi- 
cated in parenthesis in the comment field (M represents a "generalized" third 
body). 

2.8.3 INERT SPECIES OPTION 

Inert species (i.e. species not appearing in reactions) can be included 
in the input by input of a card with INERTS in columns 1 through 6 followed by 
a list of inert species names. The species names must each be followed by a 
comma and each name must be written exactly as in the master thermodynamic 
data. The last comma must be followed by the word END. See Table 2-14 for 
an example. The species list can continue cn to the next card, but a species 
name cannot overlap onto the next card. 
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2.8.4 THIRD BODY REACTION RATE RATIOS 


As described above in Section 2. 8. 2. 2 for the j** 1 reaction only one 
reaction rate, k^ where 

k - »t' n i ,- b / st 


can be input. For three body recombination reactions such as 


H + 0H + Wl = H 2 0 + M a 


the rate of reaction is in general different for each species, M^, depending 
upon the efficiency of the species, M^, as a third body collision partner. The 
temperature dependence of a recombination rate is approximately independent 
of the third body; i.e. , for the i th third body and j th reaction; 


K tJ = A tJ T" N j e' B J /,T . 

th th 

The third body efficiency of the i species for the j reaction is then defined as 


Hi a A VA. 
iJ ij J 

Thus rrij. i r, the ratio of the reaction rate with species M. as the third body to 
the reaction rate input on the reaction card described in Section 2 .8.2. 2. 

If reaction rate ratios, m.^, are to be input id tlie dissociation- 
recombination reactions, a card with TH'RD B0DY REAX RATE RATI0S in columns 
1 through 27 must be input next. If this card is deleted from the input, the 
program assumes all rr^ . -• 1. If this card is included in the input, it must be 
followed cither by i card with ALi !QUAL 1.0 in columns 1 through 13 (which 
sets all m^ - 1) or by SPECIES cards as described below. 

The m can be input using a card with the word SPECIES In columns 1 
*J 

through 7. This word is followed by the name of the 1 species followed by a 
comma, followed by the values m^ in F format, each followed by a comma. 
These values can be continued onto succeeding cards. Note that the 
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TABLE 2- 14 LISTING OF SAMPLE REACTIONS CARDS FOR AN 0,/H, PROPELLANT 

REACTIONS 0-H HAY 3-4 197? JANNAE PSWG 

H ♦ OH « H?0 . A*7,5E?3 » Ne?,6 • P=0.« (AR) NO* 

0 ♦ H « OH * A«4,oElR » N»l« *.R“0*» (AR) NO. « 

ft ♦ 0 a 0? * A«1 »2E) 7 « NbJ, t 0«0** *AR> NO* 

H * H a H? t Ab 6.4E17 * Noj, * B»0.« (AR) NO* ' 

END TOR REAx 

H? ♦ ftH e H * H20 * Ao?.)9E13» NjiO, * R"5,1B* RAULCH NO. ? 

OH * OH a 0 ♦ H?0 « As5,75Ej?» N«0» * B»*7Bo» BAULCH NO. 

H ♦ ftH B 0 ♦ H? « Ae7*33El?« Ne<), * fl«7,30<>» RAULCH NO, 

0 ♦ OH « H ♦ 0? * Ab|*3F43 , nbq* * Bsq* * . BAULCH NO* ? 

last REAX 
JNERTS N>*AR*ENO 
THIRO HOOY REAX RATE RATIOS 
SPECIES AR»)**1*«1*.1*« 

SPECIES H?»5.«5,*5,j4,» 

SPECIES H?0|?0» *3**5, »?a, * 

SPECIES 0?.S..5..4.5.l*5, 

SPECIES N? *4, *4, *4, *1,5* 

SPECIES H*l?.Sil2.5*)2.r.,?5.» 

SPECIES 0»l?.5*l?*5»l?«5»?3*» 

SPECIES 0H*i?.S,i2,s,i2 t 5«25,* 

LAST CARO 
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*\J M 


i 



values depend on the order of input of the reaction cards, i.e. the J tfl reaction 
th 

is defined by the j card input after the REACTIONS card. 

Table 2-1 1 gives a sample input for a hydrogen/oxygen system using 
third body reaction rate ratios. In this example the three body recombination 
rates are input with argon as the third body. The rate with H 2 as a third body 
is a factor of five larger than with Ar as a third body for the first three reactions 
and a factor of four larger for the fourth (hydrogen recombination) reaction. 

At this point in the data input deck, a card with LAST CARD in columns 
1 through 9 must be input. 

2.8.5 SCfDK NAMELIST INPUT 

$<#DK Namelist input specifies the conditions for the kinetic expansion 
calculation. The input is read in subroutine 0DKINP and consists of the follow- 
ing groups of data as outlined below: 

2.8.5 .1 Inlet, Throat, and Expansion Nozzle Geometry 
2 . 8 „ 5 . 2 Integration Control 

2.8. 5.3 Print Control 

2 .8. 5.4 Species Selection and Mass/Mole Fraction Check 

2 . 8 . 5 . 5 (SDK Problem Input 

2.8.5. 1 INLET, THROAT. AND EXPANSION NOZZLE GEOMETRY 

The nozzle geometry is defined in Figuro 2-2. At the tangent point 
where the nozzle is attached to the throat section either a cone, parabola, 
circular arc, nozzle defined by a wall table may be input. Ttoo radii of curva- 
ture for the throat wail (upstream and downstream) are allowed in ODK (TD2P 
allows only a single radius oi cutvature). 

Certain of the 3<2fDK Namelist items are communicated from the SGZ0M 
Namelist data set. These items will be prefixed with a plus (+) sign. Input 
of any of these data items in the $0DK Namelist set will override the values 
input in the $GE0M Namelist in this calculation module only. 
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Figure 2-2 ODK Nozzle Geometry 


Item Name 

+ RSTAR 

+ RWTU 


+ RWTD 

+ THETAI 

+ RI 



Units 

SI Units 

throat radius 

inches 

meters 

upstream normalized 
wall throat radius of 
curvature* 

none 

none 

downstream normalized 
wall throat radius of 
curvature** 

none 

none 

inlet angle 

degrees 

degrees 

normalized inlet wall 

none 

none 


radius of curvature 


For a Conical Nozzle Option 


Item Name 


Input Quantity 


Units 


+ IWALL 


requests cone option none 


+ THETA 
+ EP 


EPS 


cone angle degrees 

nozzle expansion ratio 
(£400). If EP is input, the 
Qi DK expansion will be 
through EP. 

set equal to EP none 


Note: Parabola (IW ALL-2) and circular arc (IWAIX-3) options are available 
only through use of the $GEOM Namelist. 


♦The uunsomc analysis requires that a value of RWTU £ .5 be input. RWTU - 
RW YD if data comes from $GEGfM. 

* * If a corn r expansion (i . * * . Prandtl-Meyer fan) is desired, a value of RWTD = 
.05 is recommended. Experience has shown that values smaller than this give 
the same result but are computationally less efficient. 
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For a Contoured Nozzle Option 


Item Name 


Units 


+IWALL 


requests contour option 


THETA 


wall angle at tangent point degrees 


wall angle at exit 


degrees 


NRZS 


number of points in 
table, £20 


none 


+ PWRS(2) 


normalized radial wail 
coordinate table* 


none 


PWZS{2) 


normalized axial wall 
coordinate table* 


none 


RZN0RM 


optional normalization 
factor for input wall 
coordinate tables 


none 


2. 8. 5. 2 INTEGRATION CONTROL 


The program begins its calculations using an implicit integration 
method to integrate the fluid dynamic and chemical relaxation equations. 

For near equilibrium flow an implicit method is used since it is inherently 
stable. However, once the flow has become sufficient frozen, explicit 
methods become numerically stable and can be used more efficiencly. The 
program uses temperature as a freeze criterion in order to switch from an 
implicit to an explicit method. This criterion is T < TEXPLI where TEXPLI 
is assumed 0°R but may be input. The program always prints the axial posi- 
tion at which the switch occurs, For high area ratio nozzles this procedure 
can be especially useful. 


The integration routine controls the step size such that the relative 
error in the dependent variable Increments are less than a prescribed fraction, 
DEL. Only doubling or halving of the step size is permitted, and on option, 
either all the variables may be considered (JF*0), or only the fluid dynamic 
variables (IF=1) may be considered. 


♦PWRS(l) and PWZS(l) are internally computed at the coordinate* of the contour 
attachment point. 
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When the flow becomes supersonic and the area defined fluid dynamic 
equations are used, an additional check on continuity is applied in the form 




< C0NDEL 


where C0NDEL is an input relative criterion. 


The step size is held between the two input bounds HMIN and HMAX. 
Fixed step cases may be run by setting input values for HI, HMAX, HMIN all 


Item Name 


InDut Quantity 

Units 

Assumed 

Value(s) 

HI 

= 

initial step size 

none 

.01 

HMAX 

- 

upper bound on step size 

none 

0.10001 

HMIN 

= 

lower bound on step size 

none 

.005 

DEL 

= 

fractional incremental 
error 

none 

.001 

TEXPLI 


temperature below which 
explicit integration will 
start 

°R 


CGfNDEL 

= 

relative error criterion for 
continuity check for super- 
sonic flow 

none 

lxlO' 6 

IF 

- 0 

all variables considered 
for step size control 

none 

0 


a i only fluid dynamic variables none 

considered for step size 
control 


2.8.5. 3 PRINT C ONTROL 

Output from the Kinetic Expansion Calculation consists of complete 
output for each print station selected. The end point of the nozzle is always 
printed. Print stations are selected from one of the following options: 
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Item Name 


Function 

Assumed 

ValueM 

JPRNT 

= -2 

print throat and input area ratios 

-1 


= -1 

print at internally set area ratios for conical 
nozzle.* Print at input wall contour points 
for contoured nozzles 



= 0 

print at every integration step 



= +1 

print every ND3rd step up to the throat and 
then nominal area ratios 



= +2 

print every ND3rd step over entire nozzle 


If JPRNT is +1 or +2# 

the following must be input: 


Item Name 


Function 


ND1 

= 

first integration step to be select** lor 
print 

- 

ND2 

= 

last integration step to be selected ter 
print 

— 

ND3 


print every ND3id step between ND1 end 

ND2 . 

- 

If JPRNT is -2, the following must be input: 


Item Name 


Function 


+ ARPRNT(1) 


requested area ratios for print# oust be 
monofonic increasing and greater then 1.0 
(usually entries are the same es those 
used in SUPAR of $0DE) . 


+ N JPRNT 

s» 

number oi area ratios requested for print 

S ICO, 


An extended print option may be 

selected as follows: 


Item Name 

Value 

Function 


IDYSCI 

= 0 

no extended print requested 

0 


* 1 

extended print option selected (not suggested) 


For JPRNT and a conical nozzle (i.e. IWALL * 1), the internally aet area 

ratios are: 

ARPRNT(1 ) = 2| 3*4 # • • • # 39 #40/42# • • • # 58 # 50 # 64 #•••116# 120 # 126# • • • # IW « 21b# 



220 400,412,444. 
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2. 8. 5. 4 SPECIES SELECTION AND MOLE/MASS FRACTION CHECK 


In order to interface 0DE equilibrium calculated start conditions with a 
the kinetic expansion calculations, special consideration must be made for 
inert species (those not appearing in the reaction set) . Inerts may be selected 
explicitly by use of the INERTS directive or by use of a relative selection 
criterion. The INERTS directive is described in Section 2 .8. 2. 5. 

The relative selection criterion is described below: 

Item Name Function 

EPSEL = all species which do not appear explicitly 

in the reaction set but whose mole fractions 
are greater than the input value for EPSEL, 
will be retained for the kinetic expansion. 
Species selected under this criterion are 
treated as inert. The program assumes 
EPSEL = 1.0E-5, unless input. 

In some instances it may be desirable to use input species concentra- 
tions which do not sum to unity. Species concentrations, either input or from 
equilibrium start conditions, are summed and the sum checked as described 
below: 

Item Name Function 

XMFTST - input species concentrations are summed 

and checked versus unity using this input 
criterion. If 

j 1 - L species concentrations] < XMFTST 

then the test is passed. The species con- 
centrations will then be normalized such that 

£ species concent rations * 1 . 

'.he progiam assumes XMFTST = 1.0E-3, unless input. 

If the test is not passed, an error message wi.'l be given and the run 
terminated. 
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2.8.5. 5 0DK PROBLEM INPUT 


This input is required when an 0DK calculation is requested without 
restricted equilibrium start conditions being requested. A kinetic expansion 
from input arbitrary start conditions is to be computed. In addition to the input 
items described in Section 2.8.5* an GfDK problem requires input of those items 
described in Sections 2.8.1 and 2.8.2. 


Item Name 

Input Quantity 

Units 

PC 

chamber pressure 

PSIA 

T 

initial temperature 

°R 

V 

initial gas velocity 

ft/sec 

JPFLAG * 0 

pressure table calculateo 
internally 

none 

= 1 

pressure table input 

none 

ECRAT 

initial contraction ratio 

none 

For JPFLAG = 0 option* the following must be input: 


Item Name 

Input Quantity 

Units 

PI 

initial pressure 

PSIA 

PESTAR 

throat pressure 

PSIA 

For JPFLAG « 1 option* the following must be input: 


Item Name 

Input Quantity 

Units 

PTB(l) 

normalized pressure table 
entries" 

none 

ZTB(l) 

normalized pressure table 
coordinates** 

none 

NTB 

number of pressure table 
entries* £127 

none 

Z 

initial axial position 

none 


"normalized to input chamber pressure* PC 
""normalized to input throat radius* RSTAR 
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Two-D. ! menslonal-Two Phase Flow (TD2P) Module Input 


Data required to execute the TD2P module is read in under the namelist 
name $TD2. Certain of the input items described below are either communicated to 
the TD2P module by the GE0M, BAL # and/or 0DE modules or are preset In the TD2P 
module. Hence, these data items do not necessarily need to be input in the $TD2 
namelist set. However, any values that are input will override the communicated 
or preset values. 

The input data items to the TD2P module are as follows: 

Description Units Assumed 

Value (s) 

namelist name 


The following items are communicated from the 0DE or BAL module 
if 0DE or 3AL=1 . or 2. has been input in the $PR0B namelist. 


PR 

ss 

Prandtl Number 

none 0 . 0 

GMG0 

s 

chamber gas viscosity coefficient 

lbm/ft.sec. 

CAPN 

= 

viscosity temperature exponent 

none 

CPL 

= 

liquid particle heat capacity (T p >T p ) 

ft 8 /sec? °R 

TPM 

r: 

T , particle solidification temperate 

°R 

DHM 


m 

latent heat of melting for the particles 

if/sec? 

EPM 


expansion ratio at which solidification of 
the particles begins in an 0DE type calcu- 
lStion 

none 

UGM 


gas velocity at the 0DE solidification 
paint 

ft/sec 

TPS 

2? 

gas temperature at some point after soli- 
dification has been completed. 

°R 

UGS 

K 

gas velocity corresponding to the value 
of TPS 

ft/sec 

PC 

- 

chamber pressure 

psia 

TGO 

s. 

chamber temperature 

°R 

WPWGT 

= 

ratio of particle to gas weight flow 

none 

XMLW 

e 

molecular weight of the condensed phase 

lb/lb-mole 101.96 


On the following geometric input data, the items proceeded by 
an astcrix (*) will have been communicated by the GE0M or BAL 
module if they were input in the $GE0M namelist set or calcula- 
ted by the ballistics (BAL) module. 


Item. 

$TD2 
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Inlet and Throat Parameters (see Figure 2-3) 


tyem 


Description 

Units 

Assumed 

Value(s) 

DZI 

s 

Az, particle trajectory integration step 
sice 

none 

.002 

DZMIN 

= 

^ z min # in * et step slze P arameter 

none 

.002 

NILP 


N^ # number of initial line points 

none 

20 

*RRT 

= 

R , throat radius of curvature, A value 

RS 1 is required 

none 

- 

*RT 

= 

r*, throat radius 

ft 

- 

SAUR(l) 


first estimates of x , u , or, 8, and v 
for the special throat expansion. Re- 
quired only if 8 

none 

-•15, 1 # ,5, 
.3,-1 

THFD 

s 

Qf, faring angle ( 9 ^> 0 ^ = >no faring) 

degrees 

5. 

*THID 

ss 

R, , inlet angle 

degrees 


THIW 

s 

0 . , intersection of initial line and wall 

"tAf 

degrees 

12 . 

THJD 

s 

w 

R. # angle defining the zone farthest down- 
stream 

degrees 

9. 

VAR(l) 

m 

first estimates of x , u , a, 8 , and y for 
the zone farthest upstream 

none 

•3,0,0, 

,l,.l 

ZAX 

s 

2 axis' lntnrsection initial line and axis 

none 

- 

l\ 

■3 

number of upstream zones 

none 

3. 

ZJ 

sz 

# number of downstream zones 

none 

2 . 

NTBL 

s 

number of points in the subsonlc-tran~ 
sonic region to be output for input to 
the TBL module 

none 

15 

XITBL 

B 

starting normalized axial position for the 
output to TBL 

none 

- 

*XLSTAR 

s 

motor L*, If input in $GE0M, It will be 

in 

- 


used unless It is calculated in the Bal- 
listics module. In which case, the 
average L* from the ballistics module 
will be used, A value input here will 
over-ride previous values. 


Characteristics Mesh Control Data 

DL 

s 

&i, maximum LRC mesh width 

none 

.2 

DTWI 

M! 

&R , maximum flow angle change along 
the wall 

degrees 

3. 

DR 

sz 

Ar, maximum RRC mesh width 

none 

.2 
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Va/iable 

Input Name 

r i 

RRT 

r* 

i RT 

« f 

THFD 

9 i 

THID 

e i 

THIW 

„ w 

J 

ThJD 

*1 * 

Z axls 

ZAX 

S J 

THJW 



Item 


Description 

Units 

Assumed 

Vaiuefe) 

EW 

=s 

c , end of nozzle wall criterion 
w 

none 

.001 

IMAX 

s 

i , maximum number of iterations per 

point 

none 

5 

NI 

5S 

n 1# select each n,^ LRC for print 
n 2 , print each n 2 point on selected 

none 

1000 

N2 

ZK 

none 

1 

Nozzle Wall Contour Data 

•IWALL 


Option flag 

0 = > tabular input (not allowed in $GE0M) 

1 * >cone 

2 * > parabola 

3 * > circular arc 

4 = >wall contour (spline) option 

none 



If the wall Is to be Input in tabular form (rWALL*0) 

• 

• 


PW(1) 


(rj, Zj) , wall coordinates 
l = 1,2, . n points, viz.: 

PW(l)«r 1( 2j, r 2 , z 2 r n , z n , 0,0 

none 




Note: 





a) always mark the end of the table with 
two zeros 





b) n * 74 is required 





c) always set r.^l+RRTU-cos THIW) and 
z.-RRT sin THIW 



•EPS 

at 

E, nozzle expansion ratio 

none 

- 


If a ccnc, parabola, or circular arc contour Is to be specified 


then: 




•THJW 

a 

6, , attachment angle for the contour; 

J w 

e.g. for a cone, the conical half angle 

degrees 


•BPS 

at 

c, nozzle expansion ratio (cone only). 

none 

- 

•RWMAX 

at 

r , nozzle exit radius (parabola or arc 

m. 

none 

- 

•ZWMAX 

a 

z , nozzle length from throat to exit 
(parabola or arc only). 

none 

- 
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Description 


Units 


Assumed 
Valu 



IT0T 

SMP 

R(l) 


number of particle groups to be considered none 
particle density lbm/ft 3 

r , the radius of each of n particles Is ft 


WPWT(l) 


r , the radius of each of n particles Is ft 

to^be Input so that r p < r p . . .< r p 

set r = 0. n < 10 Is retired. n 
p n+l 

w AIw # particle weight flow fractions none 
p j p j 

corresponding to each of the above particle 
radii, r 


5 

250. 

0.0 


Note: The mean particle diameter and distribution will be calculated from the 
following data if R(l) 3 0.0. The equation 


d p - XK • PC PEXP . 
*(l+24'XD.r*) 


(£)XIEXP (] -tXP(XL*XLSTAR)) 


is used to calculate the mean particle diameter in microns. 


PEXP 

AL 

XD 

XILXP 


constant in the & calculation 
P 

constant hi che calculation 
P 

constant in the d calculation 

constant in the d_ calculation 
p 


microns 

none 
. -1 


in 4 
none 


0.454 

.33334 

-.004 

.045 

.33334 


XILXP » constant in the d p calculation none .33334 

In the abovo formula: r is the concentration of condensed phase tmoles/100 gr products) 

calculated from WPWGT and, PC is in psla, r* is the throat radius in feet. 

SIG * geometric standard deviation used in none 1.9 

the particle siz e distribution 

SEND end of $TD2 namelist set 

*PWRS(1) and PWZS(l) are internally computed at the coordinates of the contour 
atta hment point. 
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2.10 


Data required to execute the Turbulent Boundary Layer (TBL) Module is 
read in under the namelist name $TBL. Certain of the input items described be- 
low are either communicated to the TBL module by the 0DC and/or TD2P modules 
or are preset in the TBL module. Hence these data items do not necessarily need 
to be input to the TBL module. However, any value that is read in will override 
the assumed or ..ommunicated value. 


Ihft input data items to the TBL module are as follows. 


Item 


Ds&sfleiigg 

Units 

Assumed 

valMfife) 

$TBL 


namelist name 




The following items are communicated from the 0DE module if 



0DE« 

1 or 2. 



TO 

» 

Stagnation temperature (from 0DE-chamb- 
er temperature) 

°R* 

- 

PO 

* 

Stagnation pressure (from 0DE-chamber 
pressure) 

psf* 

- 

PR 

S 

Prandtl numl>er 

none 

- 

ZMUO 

m 

Valuo of viscosity, ju , at the stagna- 
tion temperature 

lbm/ft-sec.* - 

ZMVIS 

■ 

Exponent in the viscosity-temperature 
equation 

ZMU0*(T/T0) ZMVI5 

none 

‘ 


The following items are communicated from the TD2P module 
if TD2P • 1 or 2. 


PO 

• 

Stagnation pressure (from TD2P-the value 
of stagnation pressure at the initial line) 

psf* 

GAMO 

* 

Stagnation value of specific heat ratio 

none 

REAR 

m 

Value of gas constant 

ftVse<^°R* - 

DCTAB 

m 

Number of points In the x, y and Mach 
number tables 

none 

SCALE 

a 

If it is desired, non-dimensional values 
oi x and y may be input and then multi- 
plied by a single scale factor, (from 
TD2P, SCALE*!*) 

FT* 1.0 

TO 

m 

Stagnation temperature (from TD2P the 
value of the stagnation temperature at 
the initial line). 
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Item 


Description 

Units 

Assumed 

Value(s) 

WD0T 


Mass flow rate 

lbm/sec 

- 

XITAB(l) 


IXTAB values of x, axial distance , in 
monotonically Increasing order (from 
TD2P-XITABfl)-Xj/r*) 

IXTAB values of y, nozzle radius or con- 
tour height in consecutive order corres- 
ponding to each x of XI TAB. (from TD2P- 
YITABQhrj/r*) 

none* 

- 

VITAB(!) 


none* 


ZMTAB(l) 


IXTAB values of free stream Mach number 
in consecutive order corresponding to the 
contour points x i# y^ 

none 



The following items are preset or must be Input to the TBL 



module 

• 



MZETA 

s 

Velocity power law exponent, n 

none 

7 

IPRINT 

s 

0 Prints omy at input intervals 

1 Prints at all calculated points 

none 

0 

ICTAB 

* 

0 Indicates a constant specific heat 
calculation 

none 

0 



* 2 ? 3 Number of points in C versus 

T table p 



I TWTAB 


-l T w * adiabatic wall temperature 

0 T w - constant (see TWTAB) 

1 IXTAP values of T will be innut 

(see TWTAB) W 

none 

-1 

ZNSTAN 

B 

Interaction exponent n, in Stanton number 
relation 

none 

.1 

DXMAX 

s 

Maximum length of stop size* If DXMAX 
is not Input DXMAX i» computed to be 

1/100 the length of the contour 

ft* 

• Olx 

EPS 2 


Geometry indicator 

0, Two-dimensional Planar flow 

none 

1.0 



1,0 Axisymmetric flow 



FJ 

s 

Conversion ! actor between thermal and 
and work units 

ft-lbf* 

TJTTT 

778.2 

G 


Proportionality constant in equation F* 
(m/g)a 

ft-lbm* 
sec 11 ibf 

32.174 

TCTAB(l) 


ICTAB values of temperature in increa- 
sing order, defining the C p -T table. 

The maximum T should exceed and 
the minimum T should be below the low- 
est temperature likely to be encountered 

In the system 

°R* 



-> t 4 

-"Ot 



Item_ 


Description 

Units 

Assumed 

Value(s) 

h CPTAB(l) 

ST 

ICTAB values of C in consecutive order 
corresponding to the values in the temp- 
erature table. 

BTU*/ 

lbm°R 

- 

TWTAB(l) 


If ITWTAB = -1 omit this input 

If ITWTAB = 0 input one value of wall 
temperature 

If ITWTAB = 1 input IXTAB values of wall 
temperature corresponding 
to the wall locations of the 
x and y tables 

°R 


T0LCFA 

= 

Tolerance in Cf-CfRQ interaction loop 

none 

.0001 

T0LZET 

= 

Tolerance in zeta interaction loop 

none 

.0003 

T0LZME 

= 

Tolerance used in gas property eval- 
uation loops 

none 

*— < 
o 

1 

■vi 

THETAI 

mmm 

Initial value of momentum thickness. If 
THETAI is input as a negative value, the 
sonic point start procedure will be actuated 

• 


PHII 

= 

Initial value of energy thickness. This 


- 


need not be input If the sonic point start 
procedure is used 


Edge condition option flag 

IFEDGE = 0 edge conditions are computed from the none 0 


PITAB(l) 

Mach number table 

“ 1 tables of P^ , T /T c , C , V,, 

are to be input as Below 

> 2 tables of P , T . C , V , o arc to 
be input as e beiOw p c 

~ IXTAB values of P /P_. or P rt versus X (see 


TITAB(l) 

IFEDGE) e c e 

= IXTAB values of T /T or T. versus X(see 


+ CPITAB(1) 

irEDGE) e c 

~ IXTAB values of C versus X 

BTU/lhm* - 

Pc 

J R 

VITAB(l) 

- IXTAB values of versus X 

ft/sec* 

R0ITAB(1) 

= IXTAB values of versus X 

Jbm/ft 3 * 

SEND 




’TNOTE-Any consistent set of units in accordance with the values of FJ and G may be 
used. The units shown are those conslstant with the values communicated from the 
0DE and TD2P modules and by the TBL module back to the SPP control routines. 

*NOTE-If a C table is not input into the program (see CPTAB and ICTAB) , the pro- 
gram %;ill attempt to edit in a C table from the values of Cp e vs. T e . However, 
the edited values may not (especially in the cold wall case) cover the entire range 
of values required for the mn. The safest procedure is to generate by hand a CPTAB 
table which includes several values of C p for both the upper and lower temperature 
extremes. 



s 


! 

1 

i 


\ 

i 

i 

j 

i 
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2.11 Linkage 


The linkage between the five basic computational modules in the SPP com- 
puter progrrm was designed to give the user of the code flexibility in selecting 
which of the modules need be used In any one computer run. While all of the 
modules may be used in any one computer run, this is not usually the most efficient 
approach since input mistakes or the use of the program in parametric studies 
could require the repeated successful execution of several of the modules. Thus, 
for the sake of efficiency, and overall flexibility, all of the linkage data between 
modules (except for the restricted equilibrium interface between the 0DE and 
0DK modules) is written out on logical unit 8 as formated data. This data may then 
be manipulated in a number of fashions depending on the computer hardware and 
utility routines which are available. On 6000 series CDC equipment which the 
program was developed for, logical unit 8 was equivalenced to the formated punch 
file. Hence, at the completion of a run, the linkage data would be punched out 
and made available for subsequent runs (described below), eliminating the need for 
repeating an already successful solution. 

This type of linkage has two major advantages. The first is obviously a 
savings in computer run time on series of cases wherein the output from one or 
more modules is unchanged. The second is the ability to use other analyses (if 
appropriate) to ’‘fool” the SPP code into accepting them as an integral part of its 
own computational procedure. This second approach allows the engineering user 
the flexibility to substitute the output from any other code or analysis into the 
procedure defined in the SPP coot as long as it conforms the proper format. 

Data is read in on subsequent runs by the SPP progiam on logical unit 3 
if the calculation module equals option Is selected in the $PR0B namelist. 
However, this logical unit may be changed to any value by input in the $PR0B 
namelist. This feature allows the punched output linkage data to be read In as 
part ci the normal card input stream au long as the linkage data is In the follow- 
ing order: 


$PR0B 


module option = 2, 

$END 


0DE 

data* 

BAL 

data* 

0DK 

data* 

TD2P 

data* 

TBL 

data* 

the rest of the input data 


♦this data is only needed If the module 
option *2 Is selected in the $PR0B namelist. 

The data does not need to conform to the above order a l long as the data resides 
on logical unit 3 or any other logical unit which is not assumed to be the card 
input stream. 

Certain complications can arise when previously calculated (punched) 
data is used out of the ordinary calculation sequence presented in Table 2-1 
of Volume in or Figure 2-2 of Volume II. For example, if the 0DE module is ex- 
ecuted while the ballistics module (BAL) data is read in from cards, the chamber 
pressure from the 0DE module is used in the TD2P module instead of the average 
pressure calculated by the BAL calculation. While nominally small, the difference 
between these two pressures can cause shght differences ir the performance re- 
sults calculated by these modules. In any event, this deficiency v.ill bo cor- 
rected in any further versions of the SPP code. 

The punched (or formated output) from each of the computational modules 
conforms to the same basic fonr.ot. The f irst card image contains the name of 
the moduie # whilo the last image contains the word END. The data between these 
curds consists of sequenced "packets” of information which contain the linkage 
data. A packet may consist of one or ir.oie data cards, cl! of which have the 
same sequence number. Table 2-15 contain.’ rhe basic format used in generating 
the linkage data. Tables 2- 16 through 2-20 vicscribe th* specific data transmit- 
ted between modules. 
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Table 2-15 Basic Format of Linkage Data 


Card 

Format 

Description 

first 

all others 

last 

A4 

5X, 5E14.6, 15 

A3 

name of module 

data and packet sequence number 

the word END 


Note that In the following five (5) tables the first and last cards have been 
omitted and that the five fields refered to are the five floating point data fields 
described in Table 2-15. 










Table 2-16 0DE Linkage Data 


"Packet" 

Field 

Description 

1 

1 

- reference viscosity - Ibm/ft. sec. 


2 

ce - viscosity exponent in u a M re f (T/T 


3 

T c - chamber temperature - °R 


4 

? c - chamber pressure - atm. 


5 

p c - chamber viscosity - lbf • sec/ft 3 

2 

1 

* c - chamber gas thermal conductivity - lbf/sec-°R 


2 

y __ - chamber equilibrium ratio of specific heats 
eq 


3 

B - erosion parameter 


4 

c th ~ theoretical c*=P c A*/m - ft/sec 


S 

c - chamber equilibrium specific heat - BTU/lbm °R 
p c 

3* 

1 

c melt area ratio at solidification 


2 

u melt ” velocit y dt solidification * ft/sec 


3 

T me i f " temperature at solidification - °R 


4 

y - fraction of condensed phase to gas phase 


5 

* H melt * hedl ^ ’ Tielt l ng °* condensed phase (ftVsecVlbm) 

4* 

1 

c - specific heat of the liquid - (ftVsec^/lbm °R 

p i 


2 

c - specific heat of t^e solid - iftVsec^l/lbm °K 

P. 4 


3 

M - molecular weight of the condensed phase 

cp 


4 

T It * temperature at the exit plane - °R 

t 

I 

5 

U exit ” ve l° c i t V at the exit plane - ft/sec 

5 

i 

M, v - molecular weight J the gas only 

’ gas 


2 

M - ino*e. dor weight of the mix 

nix 


3 

* effect iv molecular weight of the mixture 

w eff 


4 

Pr - chamber Prunutl number 
c 


5 

I S p - theoretical I CJ> at the exit plane 

6 

1 

I.„ - frozen I__ at the exit plane 

spj sp 


*Not used if no condensed phase is present. 


j 
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Table 2-17 Ballistics Linkage Data 


Packet Card in 
"Packet* 

Field 

Description 

all 



1 

1 

t - time - sec 


2 

P c - chamber pressure - PSIA 


3 

P - maximum pressure to this time - PSIA 

max 


4 

R* - throat r.;dius - In 


5 

L* - motor L* - in 

2 

1 

rf * - fractional c* efficiency 


2 

m - mass flux - slugs/sec 


3 

m./m - ratio of insulaticm to propellant mass flow 
^ P 


4 * 

P ci - head end chamber pressure - PSIA 


5 

fmdt - mass expended to this time - slugs 

3 

1* 

F - thrust - lbf 


* not currently used 
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Table 2-18 0DK Linkage Data 


Packet 


UN* 


Field 


N+UEND 



Description 


€ or - area ratio corresponding to I in field 2 

Kt) Sp 

1 - "Restricted" equilibrium I 

s PrE s P 

- area ratio corresponding to I in field 4 
I - "Restricted" equilibrium I 


c K - area ratio corresponding to 1 in field 2 
r - "kinetic" I as calculated by flfDK 











Table 2-19TD2P Linkage Data 


"Packet" 


Field 



Description 


- discharge coefficient 

R - gas constant, ftVsecV’R 

% - specific heat ratio as calculated by TD2P 

- specific heat, ftVsecV^R 

m - mass flux, lbm/sec (for single nozzle in mul- 
tiple nozzle cases) 

P - stagnation pressure at the TD2P wall initial 
point, psi 

T stagnation temperature at the TD2P wall Ini- 
tial line point, °R 

fj/r* - normalized nozzle wall radius 

z./r* - normalized axial distance 

M, - Mach number cm wall streamline 

^ b,7 TD2P 1 


f, TD2P, 


able 2 -20 IBL Linkage Data 


Packet 


Description 



- wnll radius - ft 


1 * decrement due to turbulent boundary layer 


•An ^xlsyr.i-vtric nozr.le been assumed here since output will be for A/ A* 
ec '. g.. 1 where A/A* (r *r* v \ 



2-7? 













3 . OUTPUT DESCRIPTION 


As a result of the many types of analyses performed by the computer pro- 
gram considerable output can be generated. Since the amount, and nature, of the 
output will vary over a wide range, as a function of th»» options selected, the out- 
put descriptions will be keyed to the individual overlays (types of analyses). The 
sample cate output in section 5 is a valuable adjunct to the output description, 
and should be liberally referred to. 

3.1 Overlay O (Program control and summary) Output 

Since all of the program Input cards do not appear in formated output 
statements, it is recommended that a card image listing of the input cards be 
obtained, via control cards, at the initiation of the run. 

The first output from the program is the case title. This is followed by 
an acknowledgement of the program directives that have been considered (GEOM, 
TRAJECTORY, etc.). After acknowledging the TRAJECTORY directive the trajectory 
information is output in tabular form. 

Following the PROBLEM directive, the proqram enters into execution of the 
various overlays as directed by the SPROB input. The output produced during the 

execution of These overlays is discussed in subsequent ‘ubsectlons. 

Upon completion of each of the "'verljys fron ,0 t-' *j.Q, the following 
outputs utc produced. 

QVEFLAY 3 . If BAL « 1 or 2, tV calculated values of F, (psi), L* (in.), r* (in.). 
df #/ dt (in/ sec). M uiugsAer), and f j # are primec »s .* function of time in sec- 
onds. * tifcle of the average, xirr.u^. md nirr. urn A * .--' u of the above quanti- 
ties is ulnv jjtput. T^e average values ' otj jteu tine., on the total calcu- 
lated bum iim» unless the following message iu printed out below the table: 

VALVES AT TIMES GREAT r *R THAN HAVE BEEN EXCLUDED 

FROM ThE AVERAGE CALCULATION 

Ihc excluded point* (should there be any) represent points during tail-off that might 
produce spurious , # values that unduly affect the calculate average f| c# value. 
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OVERLAY 4 , If QDK * 1 , or 2, a table of the one dimensional kinetic specific 


impulse* ISPqpk* and the kinetic ISP efficiency factor, are printed out as 

a function of nozzle area ratio. 

OVERLAY 5 . If TD2P * 1, or 2, a table of the calculated two dimensional, two* 
phase, perfect gas, specific impulse, ISP^^p* and the i>D, two phase, 1 
efficiency factor, iB P rlr *ted out as a function of nozzle area ratio. 

If RSD01)*O in $ the SGEOM namelist, or if Overlay 3 (Ballistics program) 
has been executed and indicates that the throat erodes, then the nozzle expansion 
ratio will vary as a function of time. 

In such cases and *sp are tlme dependent. The program cal- 

culates and prints out the time averaged values of these quantities. 

The program also calculates and outputs a value of t0 ** used * n 

the overall efficiency calculation. When there is throat erosion the use of this 
modified efficiency corrects the predicted overall delivered specific impulse for 
the fact that I vaiies with expansion ratio (see subroutine SUMRY). 

OVERLAY a . If TBL * i, or 2, a table of the incremental loss in I , aI 

5 P *1 TOT 

is output as c function of nozzle area ratio. 

After ail of the individual loss calculations have been completed a sum- 
mary page is printed out. Figure 3-1 reproduces a typical output summary page. 

A series r l descriptive notes are given below. T^ese ci*«n be correlated with the 
circled numbers on Figure 3-1. 

Not* 4 : 1. These values are calculated whe n the relevant parameter in the $PROB 
input equals 1 or 2. 

2. The calculated 2-12. two ise efficiency Mn includes the so-called 

erosion loss. The magnitude of the cr oon loss (in terms of an efficiency) is 
equul to thi ratio os (correct* \ :o, ihmat erosion) to (time averaged' . 

iaee summary output from ?D2P module). 

3, A ‘•calculated 1 * value of combustion efficiency *s obtained only if 
BAJ. ^ ! , or l and 7D2P ■ 1 , or 2 in the SPROB namelist. The "calculated** corn- 
bur: tirm efficiency is equal to the nozzle discharge coefficient (calculated by 
TD2P) times the “empirical** combustion efficiency. 
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4. This efficiency cannot be theoretically calculated with the present 
program. 

5. The calculated 2D 2 phase efficiency includes this effect. 

6. Empirical efficiency factors are always calculated. (Except for 
empirical combustion efficiency which is calculated only when BAL * 1, or 2, 
and TD2P * 1, or 2). 

7. Includes the two-phase efficiency only. If BAL ft 1, or 2 this efficiency 
will be calculated only if AVELS in $GE0M Namelist is ft 0. 

8. If die real nozzle geometry has been altered as a result of the restric- 
tions imposed by the TD2P program, the output value of this efficiency will be 
erroneous when the nozzle attachment angle (supersonic) implied by the input 
nozzle contour differs from the actual nozzle attachment angle. In such cases the 
correct value may easily be calculated by hand. 

9. These values are equal to unity unless the relevant $PROB parameter 
equals 3, and a non-unity value has been input. 

10. For each of the individual loss mechanisms, this column give* die 
efficiency factor selected for use in tho overall efficiency calculation. The 
selection criteria is as follows, in order of precedence: 

1. Input values. 

2. Calculated vaiues. 

3. empirical values 

11. The product of the selected efficiencies is not. in general, equal to 
the product of the efficiencies calculated by any of the methods, as it may 
contain efficiencies calculated by different methods. 

12. This column is i key to the .ethod selected for use in the overall 
efficiency calculation. The meanings of the ietters arc as follows: 

I * Input 

C - Calculated 

E - Empirical 

N - Not considered (a value of unity 
is used in such cases) 

13. The values of delivered vacuum specific impulse are equal to the 
theoretical specific impulse at the initial (t«0) nozzle expansion ratio multiplied 


W 
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Figure 3-1. Sample Summary Page From Program Output 




by the appropriate product of efficiencies, and then decremented by the corres- 
ponding turbulent boundary layer loss. 

14* The empirical delivered vacuum specific Impulse does not Include the 
so-called erosion loss-see note 2 for how to Include this loss In the empirical 
result. 

15. The average ambient pressure Is the average of all points In the tra- 
jectory table. The average chamber pressure Is obtained from the summary table 
following the ballistics calculation (or Is Input If a ballistics calculation Is not 
called for). The average expansion ratio Is an average based on the whole calcu- 
lated bum time. The discharge coefficient Is obtained from the TD2F program. 

If BAL * 1, or 2, and TD2P=4 or 2, In the $PROB Namelist, a page entitled 
"Thrust - Time History" Is output at the end of a run. The table Is basically self- 
explanatory, however the following Items should be noted: 

1. The mass flow (MOOT) and the Integrated mass flow (MDOT-IN1), 
Include the corresponding total mass flows (l.e. both propellant and In- 
sulation). 

|j 2. Due to Inaccuracies in the input grain geometry the total Integrated 

mass flow may not equal the known propellant (and burned insulation) 
weight. In such cases the calculated total Impulse (F-INI) should be ad- 
justed by the ratio of calculated to known propellant weights. 

3. The average I gp Is the average delivered I gp a;.d will In general differ 

from the "I delivered to the average ambient pressure" which appears on 
sp 

the aforementioned summary page. 

4. The average MDOT may d^fer from the corresponding value In the Bal- 
listics Ave.-Max.-Min.. Table as a result of different time bases for the 
two averages. 

5. Negative I-DEL values appearing during tailoff do not unduly affect the 
calculated average, AVE. ISP - F-INT/MDOT-INT. 

3.2 Overlay 2 (ODE) Outnut 

The first line of output says that one zone Is being calculated (multi-zone 
calculations are not permitted in the present application). The input reactants 
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cards are then reproduced, followed by the INSKRT/OMIT cards (when applicable). 
The NAMELIST Input is not printed out in the program's present form. The word 
NAMELISTS and the print out "NO $ODE VALUE GIVEN. . . " should be neglected. 

A list of the species being considered in the calculation is then printed. 

Each species in the list is preceded by some identification such as J12/65. The 
J refers to JANNAF data^ The letter L refers to unpublished data calculated at 
the NASA Lewis Research Center. The number refers to the month and the year the 
data was published or calculated. 

Following the list of species is the current value of O/F. This is follow- 
ed by a listing of the enthalpies or internal energies of the total fuel and oxidant, 
and of the total reactant. Following this is a list of the kilogram-atom per kilo- 
gram of each element in the total fuel and oxidant, and in the total reactant. 

The next output yields information regarding the iteration procedure used 
in obtaining equilibrium solutions. User's Interested in interpreting this information 
should refer to Reference 1 . 

The primary output of this program is basically self-explanatory. If 
EQL - .TRUE, in the $ODE iv r jut the results of an equilibrium rocket performance 
calculation will be printed out. This output contains the following: 

1 . Heading 

2. Chamoer pressure 

3. Proportion of oxidant to fuel 

4. Thermodynamic mixture properties and derivatives 

5. Rocket performance information 

6. Composition oi the combustion products 

Items 4. 3, and t* are output at conditions corresponding to the chamber, 
throat and requested area ratios (SUPAR and SUBAR arrays). They are also output 
(first exit column) at the area ratio e: which particle solidification begins (if 
applicable). 

The aforementioned tables are followed by a listing of the trace species 
which were considered In the calculation, but which do not appear In the mole 
fraction 

A table of the n ass fraction of condensibles is then printed out. These 
values correspond to the same relative locations given in the main table. 


The molecular weights of the condensibles and gas only. In the chamber , 
are printed. The value of the erosion parameter, 6 , used in the Ballistics Module 
throat erosion correlations, is then indicated. 

The next page of output gives the results of the transport property calcula- 
tions. Viscosity, MU, frozen thermal conductivity, K, and frozen Prandtl number, 
PR are given for the chamber and at the throat and exit plane. The power law vis- 
cosity parameters calculated for use in the TD2P program are then printed, followed 
by a list of the species considered in the transport property calculations. 

If FR0Z = .TRUE, in the $ODE input a second series of output is generated 
for the results of a frozen rocket performance calculation. The output is identical 
to that described for the equilibrium solution, with the following exceptions: 

1. There ere no table entries at the area ratio corresponding to particle 

solidification. 

2. Mole fractions are frozen and hence do not va.y with location. 

3. The information concerning the condensibles and transport properties 

is not calculated or printed out. 

3.3 Overlay 3 (Grain Deslg- and Ballistics) Output 

The first output from this overlay consists of the current values of the 
variables in the $BAL NAMELIST. Values input directly via the $BAL NAMELIST 
over-ride Interne 11 y generated values. 

The next ser.es of output is a repeat of the Information contained In the 
formated Input cardt . 

The main outpu’- from the Motor Performance Module is of two distinct types: 
(see sample case) 

1. Results of the grain design calculations 

2. Results of the ballistics calculations 

Grain design results are printed at t * 0, and then each time the number of 
bums becomes equal to NB. Thus, it is possible that more than one ballistic dis- 
play can appear between the grain design displays. The first column of me grain 
design output is the station number. The number of stations is NX + 1; station NX 
appears on a succeeding page, and station NX+1 is the plenum aft of the grain. 
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or the nozzle entrance. For a submerged nozzle, the nozzle entrance Is artificially 
located at the aft end of the grain. The second column Is the axial location of the 
station In inches. The next two columns are web positions between successive NB 
increments; the first Is the beginning web, the second Is the ending web, In Inches, 
and the difference Is equal to B10/NB10. Variations down the port stem from un- 
even burning. The following two columns are port areas, in sq. inches, corres- 
ponding to these two web positions. The next column is the average local grain 
perimeter in inches. The next column Is the local volume of propellant in the In- 
terval in cubic inches. The next column Is the local bum area In sq. Inches, 
and following that Is a running sum of bum area down the port. The last column is 
a running sum of local void volumes. 

In the case of end-burning grains, where the burning follows the X-direction, 
results which seemingly vary with Y are artificial and should not be taken literally 
without interpretation. Port areas will be zero, opening to the motor area across 
the end-burning face; this Interface will propagate upstream with time. Bum 
volume and local surface area also will appear .o propagate upstream, being zero 
both upstream and downstream of the interfacial region. (Web position is really 
longitudinal even though the output may make it appear to be radial). 

In the set of ballistics oriented output the first column is the station num- 
ber, followed by a loop counter. This is followed by local values of gas density 
(lb/in 3 ), mass flux contribution (lb/ln 2 -sec), gas velocity (in/sec), pressure 
(ib/ln*), port area (in*), weight flow contribution (lb/sec), burning rate (in/sec), 
cumulative web (In.), insulation area (in 8 ), bum area (In 8 ) and the axial position 
of the station (In.). (Again, In the case of end-burner, values appearing upstream 
of the burning surface are artificial). 

The sequence of output appearing after the columns are time (sec.), vacuum 
thrust (lb.), ambient thrust (lb.), delivered specific impulse (sec.), motor stag- 
nation pressure (ib/in 8 ), nozzle exit pressure (lb/ln # ), nozzle exit velocity (ln/sec.). 
Insulation area exposed (in 3 ), burn area (in 8 ), throat area (ln # ), head-end pres- 
sure integral Ob-sec/in 8 ; ignore for end-burners), aft stagnation pressure integral 
Ob-see/ln r ), insulation weight expended (lb.), propellant weight expended Ob.) 
und total impulse Ob-sec.). Ail integrals are up to the particular Ur* m question. 
This summary output is Identified by the throat nede number, NX ♦ 2. 

If this module is run in a stand-alone mode the above motor performance values 



H 
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will be the only ones calculated, hence, nozzle efficiency, ambient pressure, 
etc. must be input. If a grain design and ballistics calculation is combined with 
a nozzle performance calculation these thrust and Impulse values should be ig- 
nored in favor of those appearing in the THRUST- TIME HISTORY (see section 3.1). 

The unlabled output (6 columns) which may appear from time to time Is 
debug print which was never removed (see subroutine MAIN3). The statement "Ref- 
erence To Uncoirelated Table Number" also is debug print and states that the de- 
signated table would have been used had the associated option been exercised. It 
is not an indication that something is wrong. 

3.4 Overlay 4 (QDK) Output 

If IREQ * 1 in the $PROB NAMELIST (and ODK*l) the first output related to 
the kinetics solution is not from Overlay 4, but is a special equilibrium solution 
referred to as the "RESTRICTED EQUILIBRIUM OPTION. M The format of the equili- 
brium output has been previously described in section 3.?. 

The equilibrium output is followed by the values cf temperature, pressure 
and velocity at the equilibrium (starting) contraction ratio (ECRAT). The species 
mole fractions at the initial station are then listed. 

If IREQ * 0 the above output is not obtained. The output described below 
is obtained in all cases. 

The first output from Overlay 4 is a listing of the input reaction and third 
body reaction rate ratio cards, as read. This is followed by a table of all the 
species contained In the input reaction set. A reaction Uoio, containing infor- 
mation from the input reaction cards is then printed. 

The next output is a table of all of the species that will be considered in 
the kinetic expansion rotation. The number of the species, and its initial mole 
flection is also givon in the table. 

A formated table of reaction rate ratios (or the message - ALL REACTION 
RATE RATIOS INPUT AS 1.0) is then printed out. This is followed by a me.- sage 
relating the type of nozzle geometry option which was selected and a table of the 
relevant wall geometry. The two values alpha and alpha-bar are the particle to 
gas weight flow ratio and the mass fraction of the condensed phase, respectively. 




The next output describes the flow properties, nozzle geometry and chemical 
composition at the initial expansion ratio. This output is self-explanatory. The 
same quantities ate then printed out at a series of axial stations. The number «nd 
location of these outputs being controlled by the print flag (JPRNT). 
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When the temperature becomes equal to, or less than, the solidification 
temperature of the metal oxide being considered (if any) a message is printed out 
to Indicate the beginning of solidification. Likewise, when solidification is com- 
plete, a message so indicating is printed. 

The note EP not reached means that the final requested expansion ratio was 
not reached before the solution terminated. In some case 1 ? this is due to an input 
error which causes abnormal termination. In many cases, however, the message is 
printed somewhat erroneously. In such cases the program actually has effectively 
reached the final e;q>ansion ratio but due to round-off error the message is still 
printed. 


3.5 Overlay 5 (TD2P) Output 

The case title card is printed followed by a list of quantities which were 
either internally transmitted or input in the $TD2P NAMELIST. 

The following quantities are calculated in subroutine AGP. Only the math- 
ematical symbol and units for each quantity is given here. See subroutine AGP for 
definitions. 


Variable 

Math. Symbol 

Units 

G(I) 

V 

- 

Evl) 


- 

CPG 

C 

pg 

ftVsec 5 / 5 R 

G 

V 

- 

KCAP 

R 

fr/secV’R 

CPS 

0 

* 

ftVsecv^R 

HPS 

h 

p.M 

ftVsecV^R 

HPL 

h, 

Pi 

(tVsec* 

HPO 

V 

ftVsec 8 

PR 

Pr 

- 

RH0GO 

r go 

lbm/ft* 

B 

B 

- 

BS 

B s 

- 

GSB 

>s 
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The specific Impulse, ISP(IDOL), and exit velocity, UE(IDOL), calculated 
from a one dimensional zero lag analysis are output at the nozzle exit expansion 
ratio. The mean particle size and particle size distribution (input, or internally 
generated) are then printed out. 


The next series of output comes from the one -dimensional inlei and transonic 
throat portion of the solution (see subroutine PARTIL). The output appears in the 
following sequence: 

1. The values k and m^ for gas-particle equilibrium, 

2. Initial and final conditions for the one -dimensional inlet integration 
(see subroutine 0NED). 

3. Converged values forx Q , u Q , a, B, yand f 1# f^, f^, f^, f^ for each 
transonic flow zone (see subroutine JAMES). 


4. 


5. 

6 . 


7. 

8 . 


The corrected estimate for k, corrected values for x Q# u Q# a, ft, and 
y in the transonic zone containing the initial supersonic data line, cor- 
rected values for m„ and rrr . 

Items 2,3, and 4 are iterated twice. 

The gas-particle flow properties P . e # u , r, z, h , o. # u^ 

9 9*9 p j P J p j 

and v along the initial supersonic data line. 

p j 

The nozzle weight flow, WDOT. 

The nozzle discharge coefficient, CD. 


A message Indicating the selected wall geometr/ opticn is then printed, to- 
gether with the input (or internally transmitted) nozzle geometry. A finely spaced 
Internally calculated, table of the wall coordinates to be used in the solution is 
also output. 

The results of the* two dimensional, two phase, supersonic solution are 
printed out along left running characteristics. A header is printed, for identifi- 
cation purposes, above the output for each left running characteristic. The out- 
put may be identified as follows: 


How One: 




Item 

Header 

Mining 

Units 

L«R.C. number 

IRC 

left running characteristic number 

none 

ldent. number 

10 

type of point (see below) 

none 
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Item 

Header 

Meaning 

Units 

r 

R 

r position coordinate 

none 

z 

Z 

z position coordinate 

none 

M 

MACH 

Mach numb*. 1 * 

none 

T 

Q 

TG 

gas temperature 

°R 

V 

g 

VG 

gas velocity (scalar) 

ft/sec 

•« 

THETA-G 

streamline angle 

degrees 

v\ 

TG/TGO 

ratio of gas temperature to chamber 
temperature 

none 

'A, 

PG/PGO 

ratio of gas pressure to chamber 
pressure 

none 

V 

DG/DGO 

ratio of gas density to chamber 
density 

none 

n 

V „ /„ 

- Pj c g 

SDK/DG 

ratio of total particle density to gas 
density 

none 

C F 

CF 

thrust coefficient 

none 

l S p 

isr 

specific impulse 

see. 

iteration no. 

IT 

number of Iterations required 

none 


Rov c Tv yj. throu g h 


A iow is printed ter each particle size present at the point In question. 


Item. 

Header 

Meaning 

Units 

k 

K 

^article rise number 

none 

t . 

REK 

particle Reynolds number 

none 

V 

K 

VPK 

particle velocity (scalar) 

ft/soc 

K 

THETA-K 

particle streamline angle 

degrees 

\ 

% \ 

TPK 

particle temperature 

*R 

DFK 'DG 

ratio of particle density to gas 
density 

none 
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ritdifa 




m 



Item Header Meaning Units 

p /p DPK/DPO ratio of particle density to chamber none 

p k p o particle density 


r RPK particle radius ft. 

p k 

The table below relates the printed identification number to the type of 


point calculated: 

ID 

1 

2 

3 

4 

5 

6 
7 


Type of Point Calculated 

initial line point 
interior point 
axis point 

Kth particle boundary point 
wall point 

point Inserted on the previous L.R.C. 
point inserted on the R.R.C. 


The quantity IP printed out at the end of a TD2P solution is the number of 
mesh points on the last left running characteristic. 


3.6 Overlay G (TBL) Output 

The case title is printed, followed by a list of the current values of the $TBL 
NAMELIST variables. Values input in $TBL over-ride internally transmitted values. 

A table of the input (or Internally transmitted irem T r ;2P) wall geometry and 
trail streamline mach numbet is then printed out. 

At every printout station the program then prints t groups of quantities. 
(Note: DELF and DELI were added to the TBL output and do net conform to their 
respective group headings. These two quantities arr; not output until the wall 
slope becomes positive.) 

Contour Properties 

X Axial Distance 

XL Arc length of contour up to X 

Y Radius or height of contour 

DY/DX Slope of contour 
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Flow Properties 

M 

Mach number 

TE 

Static temperature 

TW 

Wall temperature 

TAW 

Adiabatic wall temperature 

DM/DX 

Mach number gradient 

UE 

Velocity of fluid 

PE 

Static pressure 

Boundary Laver 

DELTA 

Velocity thickness 

DELTB 

Temperature thickness 

DELT* 

Displacement thickness 

THETA 

Momentum thickness 

PHI 

Energy thickness 

H 

Shaoe factor 6*/0 

DEL F 

Thrust decrement due to boundary layer effects 


Heat Transfer 

HG Heat transfer coeffirient 

QW Local rate of heat transfer to wall 

SUMQ Integrated heat transfer rate to point X 

FORCE Drag force in axial or X direction 

LAT. F Force normal to X direction for two-dimensional planar flov/ 

DEL I I ^ decrement due to bounc ary layer effects 

Internal Integrals 

ZETA. II, 12. 13. 14, I?. 16. 17. IIP, UP, I3P 


Coefficients 


CF 

CH 

RETH 

REXL 

REPH 

RED* 


Skin friction coefficient 
Stanton number 

Reynolds number based on momentum thickness 
Reynolds number based on arc length 
Reynolds number based on energy thicknesj 
Reynolds number based on displacement thickness 
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In addition. If the sonic point start procedure Is used, the Initial values 
of momentum thickness, THETAI, and energy thickness. THH, are printed out. 

After the last printout station the value of the throat radius corrected for 
displacement thlcknoss, and a table of the normalized wall contour points cor- 
rected for displacement thickness are output. 

a. 6.1 Dimensions of Variables 

are six physical dimensions used In this program. 

M * Mass 
F » Force 
T * Temperature 
S ■ Time 
L ■ Length 
K - Heat 

Any consistent set of units may be used with this program. If Input guan 
titles are In units consistent with those below, the output quantities will have 
the units indicated below. 


REAR 

- L x F/M x T 

C i 

- Dimensionless 

c_ 

* H/M x T 


* Dimensionless 

P 

n 

- L x F/H 

) 

* L 

G 

« M x L/F x SP 

+ 

* 1 

TE 

* T 

5 

- 1 

PE 

- F,l» 

p 

» L 


■ M 'S x L 

6 * 

* l 

s 

* M/L 

QW 

* H.V x S 

UE 

« L/S 


* H/S axisynmetric flow 

Re 

■ Dimensionless 

SUN'O 

* H/S x l planar flow 

X 

■ L 

hG 

* H/V x S x T 

Y 

• L 

FORCE 

- F(F/L for 2D flow* 

XL 

« L 

LAT.F 

* F/l 


The built in values of FJ and G require 

use of the Foot, Poundp 


BTU, Poundp, Degree Ranklne, system of units unless H ai d G are input. 
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For example, using the above conventions, and using the built in values 

of rj uftd G 

M - lb-m 
f = Ib-f 
T = °R 
S * aec 
L * ft 
H -- 31 U 


R = ft lb, /lb 9 R 
* m 

c f 

* dimensionless 

C = BTU/lb S R 

p m 

C H 

* dimensionless 

FJ ■ ft Ibf/BTU 

6 

“ft. 

G = (lb m /lb f ) ft/sec 2 

9 

-ft. 

TE = °R 

A 

• ft. 

PE * lb f /ft* 

* 

6 

» ft. 

u * lb /sec it 

ro 

6* 

-ft. 

o - lb At' 
n 

<J 

» BTU/ft* sec 

UE * ft/soc 

Q 

- BTU/sec (cxlsymmetric flow) 

Re = dimensionless 

0 

- BTU/sec ft (planar How) 

X = ft 

h 

g 

- BTU/ft* sec °R 

V - ft 

FORCE 

- Lb. 

XI s. L 

LAT.P 

- Lb, 


Note: When this program is used as an integrated module in the SPP Program 
certain constraints are placed on the units by data communicated from other modules* 
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4. RECOMMENDED PRACTICES AND GUIDELINES 


This section has been included to provide a place for some of those comments 
that do not seem to belong anywhere else, but should be recorded. Thus, the follow- 
ing is a potpourri containing recom endatlons on using the SPP program, miscel- 
laneous notes, and guidelines for selecting values for certain of the input quantities. 

The input oriented comments are presented in the order In which the variables 
appear in Section 2. Amplifying comments, of foe type contained herein, for foe 
grain design and ballistics inputs have been integrated into the text of Sections 
2.7 and S. 

The best recommendation that can be made Is to encourage prospective user's 
of this program to carefully study this manual in order to learn foe capabilities. 
Intricacies, end limitations of foe program. Time spent in this way will be repaid 
several fold when foe user tries to run cases. 

Until foe user has achieved a high confidence level In running foe program. 

It is suggested that the modules be executed one or two at a time. The output can 
then be examined to see if it appears to be reasonable, and that foe case that has 
been run corresponds to the one that the user desired. Th*s approach can save foe 
user considerable computer time during his initial attempts to run the program. 

SGEOM NAMELIST 

Many of the inputs from this Namelist arc autom.n'cully transferred to foe 
relevant calculation modules. These geometry inputs can tx.- over-ridden by direct 
input to a given module, but remember that such a geometry change is known only 
to that module. 

AS UP It is suggested that some extra values of ASUP be packed 

around the exit plane, average expansion ratio, and any 
regions of particular interest, since interpolation is per- 
formed in this table. Otherwise, e reasonable number of 
roughly equally spa ced values will suffice. 

ASy§ These need to be Input only if subsonic information la 

desired (unless an ODJC solution is callad for, than 
ASUB(l) must he input equal to ECRA7). 



RCURV TD2P requires that the throat be circular with a single radius 

of curvature. Different upstream and downstream radii of 
curvature are allowed in ODK. If the throat is not circular, 
RCURV must be obtained by a best fit to the specified throat 
geometry. Values of RCURV < 1.5 should be avoided; the 
TD2P transonic analysis will fail for lower v-’ues. 

RSD0T These values are over-ridden by the results of the ballistics 

TBURN module if BAL * 1 or 2 in the $PR03 Namelist. 

THETA1 the transonic method in TD2P will fail if 0^ > 45° . 

THETA If the real nozzle shape cannot be input due to program re- 

strictions U is better not to try to match THETA with its 
actual value. Instead, one should try to select THETA to 
minimize the differences between the real and input nozzle 
shapes (see comment on THJW under $TD2 Namelist}. It 
should be understood, however, that if THETA is not equal 
to its actual value the empirical 2-D efficiency calculated 
by the program will be erroneous. 

» Remember that this table starts at 2, so the total number of 

points, NWS, is equal to the number of input points +1. 
Check the wall derivatives tables printed out in ODK and 
TD2P for wiggles. Don*t try to input too many points, or 
get spacing to close. If table is bad (not due to input 
mistake) it is usually better to throw out points rather than 
add them. 

AVELS If RAI>1 or 2, AVELS must be input In order to calculate an 

empirical two-phase efficiency. 

$PRQB NAMELIST Setting <8DE, BAL, etc. equal to 3 oilows the relevant loss 

efficiency to be input in the ETAIfl) table. This allows the 
user to bypass both the calculated and empirical options 
of the program in favor of a value obtained from a different 
analysis (if the user has reason to believe the present 
analysis is inadequate for his problem), or for any other 
reason. 

IREQ When this Has is input as 1 the initial data for the kinetics 

calculation is automatically obtained from a “restricted 
equilibrium “ solution. Unless the user has a strong reason 
for wanting to input his own Initial concentrations, etc. the 
use of !REQ*0 is not recommended. Also, If !REQ"0, the 
nrogram cannot calculate r ^ unless it is "fooled** by user 
prepared punched cards containing the I f p values that 

are normally punched following execution if a restricted 
equilibrium solution. 
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For 0DE*1 # ODKjft and TD2P/1 , the MIX array can be 
used as described in Reference 1. 

For the present application the ODE program has been 
modified so that it can execute only rocket problems. 

Thus, only P(l) will be considered. 

In its current application the ODE program has been modi- 
fied so that only one zone can be considered. 

The user is cautioned that ODK execution times become 
very large as pressure increases be/ond about 600 psi. 
Fortunately, at high pressure, the kinetics loss tends to 
be quite small (a few tenths of a percent of loss) and the 
empirical relation for the kinetics efficiency can be utilized 
without qualms. 

The ODK running times are not insignificant even at moderate pressures 
(P < 600 psi). The size of the present reaction set (Table 2-13) and the speed of 
the aluminum reactions are such that execution times on # .he order of 2-5 minutes 
(CDC6600) can be expected at moderate pressures. This, ODK should not be ex- 
ecuted indescrim inately. Its use is warranted only when the utmost accuracy 
(that the program is capable of) is desired, or when it is suspected that the kinetics 
loss may be relatively large. 

A completely satisfactory set of reactions and reaction rate data for alumin- 
ized solid rocket propellants does not exist at this point in time. The set given 
in Table 2-13 represents the results of the limited study that was possible within 
the constraints of the current effort. Should the user be Interested in modifying 
this reaction set it should be pointed out that execution time is linearly propor- 
tional to the number of reactions and proportional to the number of species to a 
power between 2 and 3. 

With the previous comment In mind the user should use OMIT cards (section 
2.6.2. 1) to eliminate as many of the unimportant (from a performance standpoint) 
trace species as possible in the ODE restricted equilibrium solution. 


SODE NAMELIST 

m 

OFSKED(I) 

DELH(I) 

ODK USAGE INPUT 


KWTD 


sa 


The kinetics program allows for different upstream and 
downstream throat radii of curvature (the TD2P program does 
not). If a combined solution including both TD2P and ODK 
modules is to be executed (and In reality RWTDjflWTU) it 
is probably not worth changing to a different nozzle shape 
for the ODK solution unless the downstream expansion is 
so rapid that the single radius (less extreme expansion) 
approximation would lead to serious underprediction of 
the kinetics loss. 
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RZN0RM 


HI 

HMAX 

HMIN 

DEL 

TEXPLI 

IF 

JPRNT 


IDYSCI 


EPSEL 


$TD2 NAMELIST 

DZI DL N 
DZMIN DTWI 


SAUR 

THFD 

ZI 

ZI 

VAR 

HFD 


DR 

EW 

IMAX 

N1 

N2 


NILP 


RRT 


THIW I 
THJD 
THJW I 


ZAX 


PWRS and PV.ZS are already normalized If they have come from 
$GEOM. Thus, RZNORM is not relevant unless ODK is run 
in a stand alone mode. 

The nominal assumed values should work well in most cases. 
User’s should not change these values without first becoming 
quite familiar with the ODK program; they can have a very 
large impact on exectuion time. 

It is recommended that these values always be input as 0 
because of the speed of the aluminum reactions. 

It is recommended that the -2 option be used. The other 
options (except - 1 ) can yield voluminous output. 

The user is warned not to change the assumed value, turn- 
ing on this option uses alot of trees. 

To eliminate species for the ODK solution the use of OMIT 
cards is somewhat preferable to EPSEL. With OMIT cards 
the user has direct control over which species are to be 
eliminated. 


The nominal assumed values for these variables should 
suffice in almost all cases. The user should be quite 
familiar with the TD2P program (and be unable to run a 
case otherwise) before venturing to change these values. 


More than 20 points should not be required, but one should 
probably not use less than IS. 

Should be greater than, or equal to, 1.5 to avoid problems 
in transonic solution. 

THIW should be selected such that the Mach number at the 
wall initial line point is about 1.1. THJD should be set 
to about 3/4 of THIW. THJW must be greater than THIW 
(the initial line must be on the circular portion of the throat). 
See comment on THETA under $GEOM Namelist. 

An estimate for this quantity may be made as follows: 

1 . 2 , 1 - cos «. where 0 , ■ attachment 

ZAX t (R, sin f S n exit angle 

This estimate is based or a first order source flow approxi 
motion. 
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NTBL 

XITHL 


IT0T 


XK 

PEXP 

XL 

XD 

$TBL NAMELIST 

IPRINT 

TH£TAI I 
PHII 


MZETA 

ZNSTAN 

T0LCFA 

T0LZET 

T0LZMC 

4 


The nominal value should suffice. 

The results of the TBL program should not be too sensitive 
to this value (for further information see Reference 4). 

A value of 3 is recommended for IT0T. A larger value , 
while possibly a little more accurate. Increases the risk 
of Impingement. 

Since the particle size correlation is empirical the user 
Is preceding at his own risk if the^e values are changed. 


Should be left * 0. IPRINT = 1 leads to voluminous output. 

The solution should not be too sensitive to the value Input 
These thicknesses decrease as unit Reynolds number In- 
creases. 

The user should be familiar with the TBL program and/or 
boundary layer theory before trying to change the nominal 
values. 
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5. SAMPLE CASES 


The Extended Delta Motor^ 0 ^ is representative of modem finocyl grain 
designs whose complexity warrants the use of a three dimensional grain design 
program. A sketch of this motor is shown in Figure 5-1. 

The input and portions of the output from two sample cases based on this 
motor are presented in this section. The sample cases serve to facilitate program 
checkout, while serving as adjuncts to the input and output descriptions. 

The first sample case utilizes all of the SSP program modules except for 
ODK. Thus, ODE, Grain Design and Ballistics, TD2P, and TBL solutions were 
obtained, as were complete summaries and a thrust-time history. The second 
sample case represents an ODK (and restricted equilibrium) solution for the Ex- 
tended Delta. A complete summary performance page based on the combined first 
and second sample case was not available for inclusion in this manual. 

The sample cases are preceded by card image listings of the input cards 
used to run them. The following discussion of the Grain Design and Ballistics 
input for the first sample case is presented to further enhance the user's under- 
standing of this module's relatively difficult input requirements. 

The first group of cards present the grain geometry input. BIO- 11. 92 denotes 
the web of 11.92 inches. NB10*2G denotes the number of bum increments; its 
significance is the number of new geometries called by the ballistics calculations 
(returns to HERCULES) without interpolation. SXOO denotes taking advantage of 
the symmetry of the motor, and 16 represents the fact that the diametric cross- 
section may be divided Into 16 pie-segments. Such may bo done with an 6-point 
star, each star segment being divisible In half. The next card defines the longi- 
tudinal grain mesh; the grain goes from X01«0 to Xil*45.72 Inches, and is herein 
divided into a mesh of NXl 1*100 Increments. Thus, ballistics answers will be 
output at 100 grain station. . NX roust be an even whole number, cannot exceed 
118, and for gn end-burner must be selected such that the mesh dimension (0.4572 
Inches in this case) is larger than the expected product RO*DT (bum rate times input 
eelculational time increment). The next card defines the lateral grain mesh; the 
grain goes from a lateral distance (radius in this case) of Y01-0 to Y1 1-20.00 
inches, and is herein divided Into a mesh of NY11-40 Increments. Note that Yll 









1. selected .0 be Iftgft then ft. mftft r.dte.l ftl. »«.!, «■“«» «»' ”» 

„U1 be inside ft. .«« In c. =1 «**■ **• »>« *» ***** ““ 

be pnpoctlonel ft ft. .codec. HX11 time. NKll. * 100 x 60 B«tl« we. ft. 1 -- 
g.st tound ft be neeeeser, ft .void eew-ftoth pce..ote-tl»e me. ontpnt. In 

complex grains . 

Ibe significance ol ft. licet oceup el needs Is ft. definition ft . oeometty 
w hlch is completely filled wlft propellent «ft divided inft . efteftetlonel meeb. 
Tb. following group, ft Cftd. define fte bom sftf.ee by sob.ftc.lo- ft volm.s 


and codes. 

Tbe first cylinder group (OKU) subtracts out s cylinder wblcb fills fte 
central port ft fte moftc. and Is tangent ft fte »•« POlnft a. soon a. fte, bftln 
at about 1/4 of the grain length. The centers of the circles bounding thls cyllnder 
are located at (0. 0. 0) and (46.00. 0. 0): (X10. Y10, Z10) and OHO. V20. ). 

respectively. The cylinder may extend beyond the boundary, as done here there- 
in the 46.00 exceeds the 45.72. Hie last card sets the radius R00=1.8B In. 

The second cylinder group (CYU) subtracts out the cylinder In which the 
subr irged nozzle throat region Is snuggled. Ae coordinates of the detlnlng 
circles are centered at (34.22. 0. 0) and (46.00. 0. 0). It may overlap the pre- 
vious cylinder, as Is the case here. A prior card inputs the rounded comer seen 
In the drawing. This is the BURN card with CR00 and 0.75 in. as the radius of 
the corner round. The radius of the cylinder is ROO-4.55 in. 

The next group labeled CONE subtracts out the conical volume that paral- 
lels the submerged nozzle exit cone. The surface of this cone, as the previous 
cylinder, is tangent tc star points. The coordinates of the ^ 

centered at (45.72. 0. 0) and (39.40. 0. 0). The radii, respectively, are RI0-7.30 
in. and R20^4.55 In. Note that a cylinder Is really a special case of a cone. 

The next cylinder group (CVL1) subtracts out the aft volume and allows the 
cRd-f.ce to bum. The coordinates of the defining circle, are centered at (45.72. 

0 0) and (46.00. 0. 0). The radius of the cylinder Is R00-12.00. (It might be 

noted that, If the aft end were inhibited, the second card would have to be a 
BURN card, with NB00 (no bum) as the proper analogy to CR00). 

(As a further note, if the aft end were. Instead. a division of segments of 
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a segmented motor, then in addition to this cylinder with a no-bum card, there 
must also be a division of mesh in the 1NIT cards. The first segment would go 
from X01 to Xll as here, with NX11 mesh. The first slot would go from XU to X21 
with NX21 set to -2 to key HERCULES for a segmented grain. The next segment 
would go from X21 to X31 with NX31 mesh, etc. The reason for the no-bum cylinder 
fitting the slot is that the HERCULES scheme needed to be modified to properly 
compute the bum-back in the slot region. Other inputs required for a slot are 
associated with Table 2-12. The new scheme was verified for the Lockheed Pro- 
pulsion Company 156-5 motor). 

The next two groups subtract out the star volumes. The first of these, 
labeled CYLI, deals with the circular cut at the front end A the star cavity. The 
second, labeled PR1S (prism), deals with the remainder. 

The cylinder group, CYLI, attempts to lay a pancake cylinder inside a star 
so that the curved surface of the cylinder approximates the curvature of the grain. 
The comer round card takes into consideration the 0.25 radius of the star valley. 
Thus, the fitting cylinder does not have a sharp edge, but a rounded edge fitting 
into the star valley. Because the axis of the cylinder is toward the perpendicular 
to the paper, the X coordinates are constant. However, the Y and Z coordinates 
are offset from the plane of the paper and the centerline of the motor as appropriate 
to lay the cylinder inside the star as closely as possible. 

The prism group, PRIS, begins at the end of tho curved surface and extends 
down the grain for a height H00-27.80 inches. Three coordinates a~e required 
to define the prism, which lays inside the star. Again, the X-values are constant; 
the Yand Z values define the triangle at the constant X position, and the H00 ex- 
tends this shape all the way down the grain where it overlaps the cylinders and 
cone at the aft end. A comer rojnd card ac~~unts for the curvature of the star 
valley. This completes the grain definition except for the outei periphery. 

The next group of cards tabulates a series of pressures followed by cor- 
responding strand burning rates. Five *ntrle* are used in this case; the maximum 
allowed is 99, so this presents no limitations. The number 5 appearing at location 
6 on the first card below the logic card must equal the number of entries. The 
bum rate table is the first table at location 1. Location 3 specifies one rate for 
each pressure. It is good practice to input burning rates over a wide pressure 


range to account for exponent variabilities, and to Include a pressure higher than 
the expected maximum* The following group of cards is a corresponding group 
for pressure exponent, or slope, which is the second table at location 1. 

The next group of cards tabulates the axial posit-ons of the propellant 
and the corresponding propellant radii. This is the third table at location 1. 

Again, 9S entries are allowed; this case used 36 entries, and this number of 
entires must appear at location 6 on the first card below the logic card. The 
initial and final axial positions, zero and 45.72, must correspond to the zero and 
45.72 appearing in the INIT cards (X10*0; Xll*45.72). No other consistency is 
required. The fourth table at location 1 is insulation input. The input is values 
of time (first card) and corresponding insulation areas exposed (second card) . 

Time Increments are the fifth table at location 1. In selecting the cal- 
culational time increment, the product of this time increment and the number of 
bums (NB in the grain design input) should approximate the expected bum time. 
This increment will Impact computer time, but is of secondary importance as com- 
pared to NX and NY. The next entry is the desired intervals for output. This 
latter time interval can be larger except where detailed study or plotting of results 
are desired; it has a trivial effect on computer time. If is is desired to have 
variable time increments in the course of the calculations, the tables may be ex- 
tended by inputting additional times and time increments associated with those 
times {c.g., finer increments for boost-sustaln transition). 

The last group of cards ere Table 2- 12 Inputs, m the order of and with the 
units given therein. The 4, 1 is a required key to end table loading. 

The first entry, TSTOP, should bo selected to be somewhat greater than the ex- 
pected burn ti: o. If TSTOP is inr.ched before the computer time reaches DTIME, 
the computations will terminate, so there is the possibility that the tail-off will 
not be computed if TSTOP is too short. The input value of nozzle efficiency in 
(15) should be 1.000 for an integrated mode calculation, and can be something else 
for a stiind-alone mode calculation. The projected nozzle area in (19) is defined 
by the annulus formed by the outer nozzle material shell and the open throat. 

Values of ambient density (23) and thermal dlffuslvity (26) are, respectively, 
4.43x10 "lb/ln 3 and 3.0xl0" 4 ln 2 /sec. The user is cautioned to remember all 
-1 completion entries. 

Other tables do not appear because they are not used. 
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Sample Case Number 1 
Input Listing 
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1 ,»•,*;<. • f 1 .*rt4655f,rt* .6',*)6F*"6 .1* T 7«9*.60 .1416587*0(1 .2*l72§r-oS 

2 .*).*. -.r. . | .Vl** r *M »62 t *‘* i r •** ,2lc63|F*0O .|76p28€«O0 .4*6871F«o5 

1 ,2>* K ' ' ■ *r» .l9T n *« r **,» ,*#kFM*.,j • 53*0*8* *o* *3Q639|f**Q *2507227*0* •8*?*?J7-|S 

y 4 .•**') .'5*6* ),44»5 59 J 9, J 1*,049 .62946 .2?^’ .?T**9 .367l* 0.0888 8.88 5 

1 ,6i4’i* r .04 .liM65* #t ? , *106?. '.04 ,lT?625F.po .1237057.00 .2517257.0* 

2 ,'V. *'••••: ,46r94’r.04 ,)02i*9f. f? % 4)?36»r.ft* ,21*5307*00 .1531147*00 ,*6*0MC.§5 


* 4 |.'*i> 

, ' ■*«»? 1,1%S 

4075.4 

*C»‘.o 

18.436 

. 6n*77 

• 1897* 

.71*91 

•170*8 

0.8888 

0.99 5 

• 


*.<*: ,55*6577 

• 04 . 

1566/ 94 *f 

? .490009' «o* 

,•.70*611*50 

•1041*57 

•08 • 

2517257-8* 


T « , ,i^**»* 

. 4 * » 1 • 1.74)1 

4*65.1 

6451.1 

22.998 

.77431 

.14999 

.19357 

0.88880 

1.3259 

212*70 5 


L® r t r * 

f M*r m 

t* 

• 3 

TW7T4-G 

TG/ir.O 

P8/P60 

OG/GGO 

<0«/08 

CP 

1*8 IT 

* 

• ill « 

VP* 


Tm* T4.« 

TP* 

PP6/P6 


0P6/0P8 


PP4 


• 1 , 4 > • } 

.*VTTO 1.1674 

4331. 7 

4.s9,) 

5.9*5 

.86117 

.1640* 

.*|4,4 

.57524 

P.9999 

0.69 1 



,4 j?407f • o* ,4**K^«,01 ,537639F*C4 .14*7657*00 .1588741*00 .25|?2S(*85 
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I 


2 [6T57f5F.Pl ,40*6**r.e4 .31?3e?2.01 ,5419982*04 # l795?7f.O0 ,1931312.00 .4*697)2.05 

3 .21 96902.62 .365260C.P4 • 1 1 1 • 1 *0 1 .S^042(*0« .2470252. PO .2667442.00 .O074T3E-0* 

* 4 ,0*77i .7*2*5 1.35,#i 5125.1 «?06,» 12.99* ,8*7o3 ,2823* ,3*9?* .69160 0*0000 0*00 5 

1 .lOpfliOF.fl .4937372*04 .101535**02 .5163942*04 .1682772*09 .1354272*00 .2517252-05 

2 .694245F.fl .459*692*04 .7959522*01 .5249642*04 .2U37)F*00 .17*9152*00 .4669712*0* 

3 .??2415F*02 .4057*02*04 .4505962*01 .5334452*0* .3111542*00 *2*04132*00 .0674732-0* 

5 4 .9*074 .70002 1 .40.32 4977.0 5710.9 16.9*3 .82256 ,2lj40 .25944 .30937 0.0000 0.00 * 

1 .2612127.01 ,5300672*04 .1375042*0? .50*7092*04 .1731452*00 *1160502*00 .2517252-0* 

2 .7024ME.01 ,4917052*04 ,1115722*02 .5137442*04 ,2160242*06 .1456202*00 ,4669712-0* 

0 4 1,27039 .04704 1.6292 4017.1 41»«,2 20.696 .79613 .17337 .21777 ,17215 0.0000 0.00 S 

1 .2090062*01 .5*»07032*o4 .l7o5022*02 .4930632*04 .1721522*00 •*?*0O5C-0l .2*17252-0* 

• 5 1.22506 ,9j9fl 1.774* 4657.0 6620.4 22.907 .76900 .14422 .10735 0.00000 1.33*9 jW.)# 5 

tOC 10 0 2 MOCM t© VO tH€t*-8 TO/TOO PG/POO 00/000 505/00 CF ISP I? 

K R£* VP* THCT6-5 TP* DP5/00 DPK/OP0 RP6 

9 1 ,47444 .61404 1.1596 5349.2 4631.3 5.364 .00250 .36952 .41060 .5560? 0.0100 0.00 0 

1 , 1 705202*41 ,434242F*«4 .464577E.Pl .5304262*64 .1460002*06 .1500262*00 .2517252-0* 

2 .669250F.41 .4475132*64 ,286J36F*01 . 5425612*06 .1745062*06 .1090342*00 .4669712-0* 

3 .2154602.4? .3644062*04 .1067132.01 .5405272*64 ,236?29C*00 .2560022*60 .0674732-0* 

« 4 .47)67 .77060 1.370? 5092.9 5375,6 13.69) .04171 ,2*063 .29139 .69799 0.0000 0*00 5 

1 .1406462.41 .5426152*44 ,1607002*0? .5)64142*44 .1606622*00 .1294402*00 .25172*2-0* 

2 .470324^.41 .4602342*04 .0671262*01 .52?122F*44 .2137112*00 .1640392*00 .4669712-0* 

1 ,?1 S642F .4? .4133362*64 .5209442*01 .5346702*04 .3155662*00 .?4???1E*60 .0674732-0* 

6 4 .OOiTO .04)39 1.6)67 4931.6 6034,1 17,838 .8l§)9 .19930 .24452 .36210 0.0000 0.00 5 

1 .1640732.01 .5417432*44 .1464572*0? .502*912*04 .1743262*06 .1167532*60 .2517252-0* 

2 ,6761*4fOI .5421512*04 .1265002*02 .5101662*04 ,2170622*00 .1304122*00 ,4669712-0* 

6 4 1 ,*6760 .96665 1.670 4750,0 6329.6 21.063 .70649 .15894 .26209 .17325 0.0000 0.00 * 

1 .?603?5r*6l .5827162*64 .1«191S2*02 .400191F-O* .1732552*60 .6097302-61 .2517252-05 

9 5 1.>M*5 1.66467 1.7966 46?8,V 6692.6 22.690 .70503 .13042 ,)0094 0.06000 1.3457 21*.00 5 

tl»e tO • 2 MACH 76 VO THETA. 0 70/760 PO/PGP 00/061 OK/OO CF ISP IT 

w R2« VP* TM2T6.6 IPO RPK/Cft OPO/OP6 OP* 

16 1 ,6P1*5 ,6?696 1.1*27 6347.7 4607,0 4,069 .80361 .37365 .42277 .§4226 0.0000 0.00 0 

1 .176674F.fl ,433?f 3f *64 .l64665F.pl .*361067*04 ,1439*6F*00 .1501322*00 .2517252-05 

? ,6447l7F.fl .4662362*64 .2507402.6) ,54lP60r*04 ,1767*52*00 .1075762*00 .4669712-05 

1 ,2117717. f? ,3637642*64 .906l62r.Pl .5409532*04 .2275662*00 .2496732*00 .007473C-0* 

16 4 .07*61 .0)063 1,4643 5661.3 6466.6 14.321 .6|640 .23586 .20106 .76356 0.0000 0.00 9 

1 . 1632**2. *1 .S|f05?r**4 .ll**vOf.r? .*|3563f*64 .1606*72.00 .1237022*68 .29)7252-0* 

P .*458!5F*M ,47*49*F • fi 4 ,614206C*61 .616.1627*64 .2149502*00 .1574032*00 .4669712-8* 

3 ,p*86 19F .6? ,42f7*2-*.04 .59*6167. r) ,5?0362f*04 .31964)2*08 .?34!782*88 .0074732-8* 

16 4 1 ,f 6*55 .06466 1,4*4. 4891,0 4645.4 10,678 ,80035 .16626 ,23045 ,3641 > 6.8080 8,80 * 

1 . |0*?4*r .6) ,45/1*47 *64 . 1 4*04 4C .0? .4960847*04 .1741662*00 .1046842*00 .25) 7252-0* 

? .*408*47. n .5123422*64 .1292267.0? . 1666252*04 .2)96127*00 .13)5002-96 ,4069712-0* 


1 * 4 1.1174(1 .6*467 1.7101 4711,7 *446,4 22. *17 ,77076 .1*778 .10*08 .1725? 8.0088 0.80 * 

1 . 10)6467 .*) .56664*2*64 .)6?400C.«? ,4840872*04 .»??6l7f.O0 .0567482-81 .25i?2$C-0* 


10 

5 

’ ,/*! *0 1.67069 

1 .0256 

4599.0 

6765.? 

22.998 

•7692? 

.13279 


1.3993 2)7.42 

9 

tor 

to 

• 2 

HACK 

76 

VG 

TH270-0 

70/766 

P0/PGO 

00/0*0 60K/0O 

CF 150 

IT 



« or- 

VP* 


THftA.6 


TP* 

0P4/0O 

or*/O00 

OP* 


11 

1 

, 4*44 7 ,* 4?45 

1.14*7 

5364.1 

4406.0 

4.27T 

.60*68 

.3772* 

.4264% ,5)198 

0.0600 9.00 

0 


1 ,1*39*67.61 .43)41 7f *04 . 9248102.61 .5166842*64 .14/4*81 *60 .1570462*66 . 2517252-6* 
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*4 


1 

2 


.475544F-M 
• ?65?4SF • dO 


,938306F*«4 

•9ft54)9E*ft4 


.177B??E*02 

«18707*E*O2 


•35p925F*0* 

•*16700E*06 


^iS7643E*(it, 

,l55387E*0d 


•513190E-02* 

.3038462-02 


•2S1723E-0S 

•4fr697lE*09 


4 *.67ni 1*. 15109 3.293* 3016,0 0006.? J3.2?7 .49846 

1 .9332! 6r-M ,97p6ft«fc>0* .1698737*0? .324«4SF«P* 


•00690 .01384 .11903 6,0000 0.00 S 

•1196?8E»09 .430051E-02 .ISl72SC-oS 


3 4.09] 7* l*.5R69p 3.427? 29|3,3 1*110.8 13.566 .40148 


•00650 .01349 0.00000 1.0320 293.99 S 


L«f |0 o z 

« ft£< 


"*CM TO VG THE74-0 TG/tGO 
VPK THE74.K TPK 


PG/PGC 0G/DG0 SOK/OG CF ISP IT 

OPK/OG OPK/OPO PPK 


45 3 4.00060 3.8)884 

1 ,?30604F*fiO 

2 .ll9O0ir.M 

3 .4)P<iWF.*l 


2.1971 4039,3 7632.0 

73*074E*04 0. 

h9P?30C*O4 0, 

6253l6E»n4 ft. 


0.000 .66757 

•421 152F*04 
.434749E*04 
•452649£*p4 


.03372 .05051 .77933 0,0000 0.00 5 

.169iA8f.oo .222032E-01 ,2$l7?5r-05 

,231 9ft5E*00 .304263E-61 .466971E-0S 

.378377E.0O .496619C-01 .99747JE-09 


43 4 2.*0*<*C> 7.86683 2.0441 3341.8 8986,? 18.455 .55230 

1 .664279E./-1 ,877696F.n4 .17?760E*0? .402099f*P* 

2 .?40766F,fifl ,045355E*n4 ,160?o n F*ft? .4107oOE 4 O4 

3 .lr6?34r.rt| ,7T173?E**4 .*6f»655E*o2 ,*)67 q«E*04 


.00668 ,01?09 .7349] 0.0000 0.08 9 

.!56636E*00 «49]956E*02 *2Sl72SE-oS 

.2101237*00 .659947E-0? ,466971C-aS 

•308] 52E*AC • 1 1S628C-A1 .947473E-CS 


43 4 4. 14**4 9.81607 3.1177 3152.1 9566,9 16,848 .37096 

1 .439814F.M ,94**B«F«*4 .l7?95PF»n2 .3475727**4 

? .24987PF.O0 ,91*124^**4 ,185943£.*2 .4**7ftOE*r»6 

45 4 4.75879 1 ; .474?3 3,1]** 2998,1 0910,4 14,922 .49549 

1 .M54»9f*M «9R]?|1E**4 ,1683fl3E.r2 ( }ff7Str* 4 

43 5 5. *7748 1“.05«*8 3. *57! 28*7.3 1*153.7 lSv#8T «477?3 


.88625 .01200 .31308 8.0000 0.00 4 

, 157702F* 00 *49] 464E**? .2Sl??5E-09 

,155379E*oo .484420C-02 .460971E-Of 

.80644 ,01350 ,11930 0.0000 0.00 S 

.1196987*00 .4191582-02 .2S1725C-0S 

.86613 .01284 0.60808 1.8374 *94.76 3 


45 4 2.05P7J 8, *7024 2.*7i9 3312.0 9f>32.2 18.755 ,i»U6 

1 .-»189i*F-*] ,8*?13*f**4 .178??*r*»!? .18l4l)r.H 

2 .248944F, rn .04908*F»*6 .16Q718F.0? 

3 .10G9^3»**1 .7784|6E**4 . 16734*F.n? 


♦3 

4 

4,??o<7 4.70)36 

3.130? 3171.9 

0Of8.4 14.955 ,51761 



1 .4l/‘7l 7F.' 1 

.9444P7F*'4 . 

1?17,5Ft? .300574 #. *4 



2 ,244r*4*.*p 

.4|45pOF**4 

1 84 ] 2 Of «*2 ,4l6Tp*«,84 

49 

4 

4,0404* « *.7*69* 

3.3*11 24T-..4 

*4-9.7 ) 4 ,8*4 ,49[ ,; 3 



1 ,7l4079F«r| 

.903035E*** 

I61453C8? .12p773F*84 

*5 

3 

5.!*?4* 11.12526 

3.4645 2064. 0 

1ft 192,* .47318 

L6f 

10 

• 2 

M.fM ?G 

7C. 7**ET/^0 TG/ T C.ft 



n »E 4 

VP« 

Tm. ‘ 4.«r IP, 

46 

3 

* . * 8 * i> f ^^.*,>7? 


7723,! P.CP'. , ♦> 4 9 * 8 



1 ,2]4»' l ^f.»* 

.T*4«**f*p4 r. 

t 4)**“'f.«4 



? , 1 .<• ) 

,*VO?.ir'.-4 *. 

,* Jf “*Pf •' 4 



7 ,78ir64l*.-l 

.O370*«f*f4 *, 

. 444|7i» ♦ .4 


46 4 7,-1777 ‘>.?9344 ?,09A* 3285, « »o71, T 19.4^ .547-8 

1 ,6i.J3’M.-l ,ft*A?OF.*4 .177**4#.*? ,39?4)!t«'h 

7 ,7'57qsf. ftW .8?4?7«C*^4 ,!’I24*0,*? ,4]8 ? *r*»'* 

7 .7«3l t '*».P* ,Vfo«A7f**4 ,)60l47r.r>? ,414" «r.* 

46 4 *,77?7J | * . *46* | 3,1**» 3110.1 1«,H§ ,*]^2 

\ 1 . 4*1*1 79F • -4 ,lA4Sv?».r? ,1*i579f.«.4 

2 , 9 ! * 4 3 Qr * * 4 .101472#. f ? .4J6TCC*-'* 

46 4 4.4**** **.«‘:|7* 3, VIS 7948,9 * - P £ 9 . f 14.6,8 ,447> r 

1 ,907*011.84 . !5*4^^.e? ,3htM1*^ 

40 5 8.10143 11,4*7)4 1,6805 2843,4 1*??5,5 J3.676 .4*499 


.99616 .01126 .74430 0.9008 0.08 S 

. I5733if.ro .468243E-0? .2S172SF-QS 

,21231 1E*00 .621081E-62 .466971E-0S 

,J74937f**0 .1O9661E-01 .947473C-09 

,90632 .rtinl ,3103# 0,0000 0*00 S 

.:5*44lf*r8 .47?702f.*2 .25172SF-0S 

.133037F O' .404896E-O2 .*4697lC-eS 

.9*434 ,c3*f . 12884 6.9088 0.08 5 

.I2«*34f.'0 .4981722— 6? .Hl7?Sc«0S 

,4PM« .0 226 0.60408 ».«61t 293,48 S 

•’3/00* 2 ft /06ft 800/96 CF ISP IT 

&047r0 ^00/000 0#H 

.0300* . t 4S4? .79010 9,9900 0.09 S 

,!7*7./ *** ,20)904*-*! .29l??3F-09 

.?>-04n. ' . 7I09337-* } ,4669Tlr-PS 

,99V a .«»699tr-0| .8674 73C-0S 

.004’) .01452 .7422% 0.0466 0*09 5 

,!V »/*•** .4>0766f-«? .9tlT?SC-0S 

,7!i4.*r.- .594 7890 —8? .666t7lf-AS 

. 30964 *.-r .1O4«14C-0| ,6674?3E-«S 

.9*0Vi ,c»l49 ,lp)4i P.0466 6.66 S 

.05*6662-*? ^S|7?Sf.fS 
,104?’! .44*«7K. r ? ,4669 ’If - qS 

.00^5 ,i7?9l *.9966 0.O9 9 

.122087C. ; .>447666*# - ,#fl 7232-of 

.00694 ,ni 70 0.00690 *.6*96 296.2? « 



44 4 3.*»332 6,51733 2.9766 3299.2 4113.1 19.243 .53949 .00930 .00903 .759 . 0.0900 0*00 9 

1 .468253F-61 .69p435F*04 .l798T*F*02 ,30fl9l6F*04 ,1379677*00 .*035917-0? .2Sl?2S£-0S 

2 .??3O40F«ft6 ,B565*3F*ft* ,17?866F.02 ,*167O0F*04 .2lS?9?F»A0 «M99R2f.02 ,*0&971F-0S 

3 .400U0F*60 ,783A4^F*04 .1689787*02 .*167007*04 ,3861?47*00 .0974737-0* 

49 , 4.41041 U. 39939 3.1535 3110.1 9624.4 19.109 .91534 .00398 .01160 ,2004? 0.0009 0.00 5 

1 .*4367or«M ,95?f 717**4 . *66991F*6? .3367137*64 ,1466647*00 «44?7l7E-6? .2317237-03 

2 .?7I036F*60 .923122004 .179257002 ,6167007*04 .1426327*00 .420094E-62 .4469717-oS 


40 4 3.00047 11.27710 3.40?9 2920.9 10090.7 J4.96? .40279 .00S7i .01102 .12450 0.0000 0.00 9 

1 .A43800F.01 .940472F«04 .1365187*62 .3190387*04 , 1244992*00 .lt?429E-02 .2317257-0* 


49 

5 


3.741*5 11.49360 

3.3277 

2027.3 

1*252,1 

15.464 .46727 

.••533 .01141 

9.00000 1.0310 296.94 5 


10 


* 2 

MACH 

TO 

VG 

TMET4-G TG/TG6 

*o/*oo oo/roo 

505/09 C * IS* 17 



5 

RE* 

V*K 


TMET4-3 

T*« 

o*x/oo 

OP 5/0*0 0*3 

47 

3 


*,66060 3,6*609 

2.7763 

3937.6 

7003.4 

o.oro .65077 

.02999 ,04*47 

.79704 0.0009 0.00 3 



1 

,1Q05*«r.oo 

. 73**10E*04 0. 


.416700C*04 

.1710342*00 

♦lOA 3147-01 .2317257-0* 



7 

,967961F*60 

.767064F *04 0 , 


.*270387*04 

.2363747*00 

•7547137-61 .4469717-0* 



3 

,350*»0?F*6l 

.941007F 

•04 6, 


,449019E*04 

.3904347*00 

.4207477-61 .097473F-0* 

47 

* 


3.i*«99 0.7?9l* 

2.943* 

3234,0 

4140,6 

19,437 .53463 

•00494 ,06924 

.79402 0.0000 0*0) 9 



1 

.932275F.M 

.89*19'F**4 

101937E**? .3833337*04 

•159054E*00 

•379422E-02 .2*17237-0* 


* ,?!07l6r*00 , r 6?5l4r**4 .17*3737*62 .4167007*04 .2192107*00 .310943C-0? .4669717-0* 

3 «8315?3F*00 .7070912*04 .1640137*62 .410700E*04 .3905171*00 .9372002*02 .007473E-0* 

47 4 *.404*17 10.43374 3.14*7 JloO.4 9643.9 14. *45 .31373 .00393 .01130 .20153 0.0000 0.00 3 

1 .730l4?r*M ,9* r 46 ’f*64 .14i70lF«0? .3361927*04 .14*4677*00 . 4333947**2 .2*172*70* 

2. .2480O2F.00 .9F*71?C*04 .l742*7f.o2 .*14700E*p4 ,1374901*00 «4l39l3E*02 ,409*717-0* 

47 4 9.47343 11.53979 3. *295 2097.* 1*049,4 1*.991 .*7949 .005*2 .01131 .12577 4.0000 0.00 9 

1 ,43893?7-*l ,9937 o7C*04 .15417*E*p 2 .3140397*04 .1257471*40 • 399964 C-* 2 .23j725E*0» 

47 9 9.334*2 11.94075 3.34*6 2*1*. 0 1*272.5 15.192 .46320 .00319 .OIK) 0.00000 1.0991 297.39 5 

47 4 4,21412 0.9*326 2.9*34 3212.1 4|66.6 19.5*6 .53000 .04462 .00070 .7490* 0.0000 0*00 3 

1 .*9?7*9F.ri ,99773?F*0* . 1 418 1 7f .1670217*04 ,1906097*00 .13727S7-*? .*51725£*C3 

* ,19709*0.46 ,«69306F*6* .1?48?9F.6* ,*)6700F*44 ,2166*97*0* ,*89961E-*2 ,4469717-0* 

1 .6***420*40 .79l**?-*04 ,l7n66HC*(i2 ,*147oor*e* ,39*1?2C*00 .69P7m-02 .9074737-0® 

*7 * *.49«»6 1 * .94 1 1 7 3.10?’. 3093.0 9671.* 1 3.9*3 ,51110 .00506 .011*7 ,279*3 4,0000 0.00 5 

1 , 7365*57. *| .9303**7**4 .l9Tof6r**? .333564F**4 ,142*&6r*40 «*2**?6l>*2 .25l725c*05 

2 ,*94l*^F*46 .93**1 *£*0* .I7?7*3f*c? t *l*0*6E*04 .l3)94Vf*06 .399224E-*2 .46997ir.o3 

47 4 3.141*0 11.0*613 3.*3|7 ?«74.* »*»73.9 1*,964 .475J9 .00517 ,01007 ,1*904 9.0000 0.00 4 

1 . * p 69*0F •<* l .9971417*04 .19*6437*02 .31*9907*04 ,1240417*09 .3*597SC-*2 .2H72SC-0* 

*7 5 9.4164* 12.7*1*2 3.4379 20*1.9 1*293.* 1*.V*4 .*6360 .09992 .01033 0.09000 1.0593 *90.29 • 

LOT to o 7 «*Cm ’0 *c. TkiT*-5 Ta/Tr.O *0/*96 09/0*9 5(34/00 C f 15* *7 

« 4f« V*ff T**C74-* TO* 0*4/00 0*4/ 0*9 **4 

*0 1 *.60009 *.*14*1 2.1177 1009.6 7403.6 o.OOC .6**4* .92**0 .037*9 .06699 0.9099 0.09 * 

1 ,i?«4!5F«*ti ,7*,?4***.*44 r. , 1 1 6 7pOF • o* a l7?ft6«r*60 . 16O30OC»6i .291 7257 -0* 

7 ,«al677«*** .?»>3?1F‘A4 *. .4>J*»*f*** .2309461*00 «?Jl)40C-6l .*669717-03 

3 .47i56ir.*i .9494917*64 *. .«4?9CS7«*4 .399)**F«09 .JOTOSK-Ol ♦•074737-63 

60 * 4 . 7471 * 9.14135 2.9019 114*.? a/ll.* 19.719 .5?7*9 .96*31 .00*21 .77*4* 9.0990 .49 3 

I .**?07ir.ft1 .9*1141* *64 .|4*5*7l«n7 .3)49*07*04 .1579*7^*06 «j3706?f-*2 .*9l7|5F-«9 

7 . 1 »* 071 V «66 , 07609 U*«** .1M/1?I.«? a *U 79 »f *64 .?17f74r»64 .*03l*77-*2 .4969717-03 

3 . 741*1 l»* 6 t ,79*69«r«A4 .1715217*6? .*14709r*e6 ,3999797*60 ,4499797*62 ••6?4?3f-03 

4* 4 *.*417? 11,7796* 3,7*** 3*7J.9 4764,* 13,39? .9*797 .•6373 .011*7 ,2737* 6,0000 0.90 3 


1 ,4979647. M .9*1 7?if«6* ,l4?5»5f.r7 .131f?»F*o« ,1*19*07*00 .*131717-0? .2S}7?5r-0S 

7 ,79094*7.60 .41*014* •««* .lA0194f.fi* ,*0«25O9*6* .1317*3***0 .404010F.47 .606971F-49 
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4.??657 

17.10066 

3.4734 

7154.7 10151.4 

14,634 

.47755 

.00497 

,0105? 

.12537 

0.0000 0*66 
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i ,34ome-M 


.3117052*04 

.1253732*30 

• 3474352* 

02 .2517252-1 

»9 

46 

5 

4.46151 

I7.571r7 

3,5639 

7796.9 10301.9 

14.462 

•46225 

.00496 

•01074 

0.00000 

1,6633 296.90 

9 

44 

4 

3.14*4$ 

5.16315 

3.0015 

3170.9 9?37.6 

19.767 

.52405 

•00406 

.00779 

.77194 

0.0000 0.00 

9 



1 .4Ul?4r.M 

,9o434lf*«4 .1670567*07 

.3760452*04 

.1574172*00 

•1165452- 

02 .7517752-1 

D9 



? ,1??85flF,nO 

.4716777404 .1745447*0? 

.4167002*04 

•7164932*00 

•4361452* 

02 .4669712-1 

D5 



3 .7141144*00 

.7946417*04 •1773777*07 

.4167002*04 

,3960?9E*00 

•6055442- 

02 .6674732-09 

4§ 

4 

4.74794 

11.53154 

3.714? 

1047.0 9757.5 

13,491 

•50356 

.00550 

•01092 

.27490 

0,0000 0.06 

9 



1 .*474414.*! 

,945?JJf *04 .1444107*0? 

•3745212*04 

,1437632*00 

.4064772- 

02 .251 7252-1 

»* 



? . 7732747*00 

,9374577 

•04 .1650032*02 

.4046692*04 

•1316702*00 

•3713432- 

02 .4669712-09 

44 

4 

5,4*444 

17.14673 

3,4494 

2644,6 10175,1 

14.607 

.47017 

•0046? 

,01025 

•17399 

,6006 6,06 




1 .1905167.01 

•100356C«65 .1514737*0? 

.3094572*04 

• 123 f 

.E*J0 

•3303652- 

2 .2517292-1 

»9 

4» 

5 

5.51473 

17.66740 

3, *647 

2792,9 10309.0 

14,300 

•46159 

• 0049? 

.01065 

6.00006 

,6652 299,22 


44 

5 

5.53797 

17.73743 

3.567J 

7794.3 lo3«9.6 

14.103 

.4616? 

.00*93 

.01066 0.00960 

,6662 299,39 


46 

9 

5.94617 

17.76405 

3.5706 

7791.4 10309.6 

14. *03 

•46140 

.01490 

•01963 4.00000 

,9667 299,49 


4i 

9 

5.49044 

17.7*564 

3.471? 

7790.6 10311.6 

14.103 

•46120 

•00469 

•OlOoO 

0.00006 

,6669 299,49 


4i 

9 

9.99301 

17.79449 

3.572J 

2769,9 10312,4 

14,103 

• 46110 

.00486 

.01059 0,00604 

.6670 299,91 


40 

9 

9.95315 

12.795*5 

3.57J3 

2746.9 10312.9 

14.1o3 

.46109 

•00466 

.01050 

0,00000 

,6671 299,51 


44 

5 

«. 55371 

17.79530 

3.5724 

2799.9 10317.9 

14.103 

•46109 

, 00406 

.01059 0.00000 

,6471 299,51 


44 

9 

4,55175 

17.79543 

3 17?4 

2749.9 !*J|?.9 

14.103 

.46109 

.00466 

.01059 0.00004 

.9971 299.51 
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SUMMARY OUTPUT FROM T02P MOOUl* 


A/A* 

2.ADO0O 

*•**000 

A.onooo 

1 *. 0 ** P * 
2* OfOOO 
3* ,0009* 
3. 93959 


J*PtTD?P» 

?2*.300T* 

2*7.33*7* 

27f.***51 
?•*. 6*929 
?9*. 9*0*3 
299 . 613*0 


•• timc avfraocd values** 

CTAITD2*) • •951*' !SP<T02P>. ?^.T»3 ISF-SEC/U* 

•• FTAIT02P1 eoRRCCTCO F0« THROAT EROSION. .f*7f» 


ETA (TD2PI 

•9*13* 

•93939 

.9*19* 

•93*02. 

.9*397 

•9922* 

•9927* 
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TURBULENT now**** l4V£i LOSS *OOUU CT»U 


rjiTrsoeo rELTA VALf 04110 N CASE 
KfFTo a ¥rtnr|1> »9nr TLF POWFP LAW EXPONENT 

IPPINT • PPTNT AT PVFPV CHCULAfEO P0lN?|«11 00 At INPUT INTE*VALS|«4) 

|IT*4 • NUMBFP OF POJNTA |r. A • »•$• V • M TABLES 

ICTiP • H or POINT* |N CP .*$. T TABLE 4 * 

jTvTiH • wall T€HP. OPTION .. AC I AHA T 1C t •• 1 1 1 CONSTANT <•«) t TABLE <«1) 

IPFOOf • FOftr COnOTTTON OPTION— C*LCULATrOI»#l» INPUT <.$€*!> 

fft « towr S1 U ?A« STAGNATION TEHPEPAlUPf 

Pa « rorr sTpfa* stag* at jon ppcSSUPE 

•APP • ATA«M*Tlf>»* PAT to OF SPECIFIC HEATS 
PPANOTt* stagnation PPAHTTL NUHBE* 

|«|f* a STAGNATION VI*CP«|TV 
} t«v 1 A • r«P 0 * 6 Nf OF VJSCOSI TV«T£HPfBATUPE 1.4* 

ZNSTtfti • AOtivOAPV «AVf« InTEpAcTIO* EXPONENT 

DIM* A • MMfMOP STfP MfF 

THETaI * TN| f ? AL VALUF OF HCmFNTU* TwfCKNESS 

PMIJ ■ T * T T I At VHUE or CNF«6V THICKNESS 

mi ■ WTPwrW.a. •AlsrMPfTPICt*!#* • Pt4Ht«pfr.» 

PPAP • 0#A CONSTANT 

FJ • CnA*VrPStON HI TWEE* THfPMAI. A HO IfOftK UNITS 
• a PPOPOHTIO* ‘LITT CONSTANT |N EQUATION •• FpN/SM 

SCALE • CONTOUR SCALE FaCTOP 


T 

t 

•• 

•I 

ft.N/»»aC«N 

S. 4B9M4<t£»B4 
) •IT-SSSEEaOS 
iiSaimiMt 
4*ti:ftMtE«09 
A.T1TS3ME-41 
l.t04MHC*91 
?»tS7444lE»it 

l.tOOMCei.tA 

momhic*m 

7*A77f 74TE*|| 

T. iimmui 
3,EI74MtC«tt 
N 79*441 4t*U 


1 

X 

•T.i3PT443t>M 

▼ 

S*4437494C«P1 

f 

•4*0i:?31Sf-' , 1 

S.S9934S4E-A8 

J 

•A.fl«*?4*15F<*M 

4. TS4fP|4E»"t 

4 

-A.S4?A4AVF*ri 

4.4P0ATSiE*A2 

S 

-S«53717*SF*f 1 

4.A*Af1A4E*A| 

« 

-i.SU7A?U-M 

4.*9|B4*.*E-M 

T 

•1.9409M7F-M 

3. 737442 3E-aJ 

• 

•1.4*ol073f*M 

3*3t*3P*SfE-A| 

• 

-?.93SAf49e-M 

3 » *2*71 1 TF*M 

14 

11 

-».*iwn63 r *M 

2.4/*34C5F.»>1 

*l # M®*®JT9F-** 

2*7224*SfF-A| 

li 



1| 

•A # 3Ar4*B4r m*f 

».««3*494f-^| 

|4 

•ji.CN4AAAir.P4 

l.TVAf*iMC-*tl 

is 

1? 

»,1AT4?Ptf-'4 

■•A1P71AAL-P? 

J •«Vfi7AQ9E.4t 


|.92T99S0C*A| 

If 

1,MAr4Aff.r» 

l*49ASA33f«A) 

it 

l.|97*?43f«r! 

t.977A?m*->J 

2* 

|,|4|7M7»«*t 

2,^JA*4SSF<PJ 

■»? 


? , 'MFA*9|F. ' 1 

ft 

»,SA(I4a/ IF .“9 

9.1 3?41A4E-A| 

t% 

1 *474441 A* •<* 1 

9.1MT7*V4F.«1 

44 

!.741«A1S*.-’ 

4,4*f74l*E-.«i 

9S 

| # 0||4?*|F»*t 

?.?M?7A?AC.*1 

4* 



ft 

9,|*7A4A)r.r.v 

2.XH449PVF.4; 

44 

F.9AG*4**f* * 


91 

*. VA9**M .P! 


)* 

9,b *4Pf Apr-# 7 

2*H**A|f 14C-AI 

If 


9.4V|4A49C*«| 

.if 


9.**\914«f.«t 

19 

• **JT?r.»! 

2.AV4444SF.A1 

it 

* %*» ’*.»<«•« - *1 

9,n«Au*F.«! 

14 

1*| 9*4AAU.P> 

9.4A4l943r*M 

If 

i.2444-4»»—l 


37 


4,«IJ*9»7F-41 


N*CH NUMtEP 
9,s?meiE*ei 
4,244A9A0E*0? 
T.121HME-M 
I.20I3THI-0I 

1. t?92S«4C>Cl 
I«9SS170#C*«1 
l.4*'2«et'-et 

2, QM760Of*0! 
V 9 Tl4|tME*P| 
3.7EJ?4eH-#1 

s, simiof*«i 

4.?»e««e4i*ei 

•.S*44l0OF-Ot 

1.94S**04E*44 

l.J74TS#0F*O* 
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l.S)9?T«»C»e9 
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i . 74T2494F *00 
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I . ?7*A*eCF*O0 
l,?444?Mf«0t 

i.4Af??#EF*4« 

l.cmiMio* 
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l.f/*4*ME*M 
I.OMITOOEaM 

7.4T2)9Mf*M 

MII0Tl9f*M 


34 

)« 

*4 

41 

4? 

O 

44 

*5 

4* 

4t 

41 

4« 

§4 

§i 

*1 

*3 

*4 

*5 

u 

w 
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4? 

43 
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47 
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44 

74 
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7f 

73 

74 

73 

74 


77 

74 
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•1 
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•3 


•7 

M 

•4 


3.57417405-01 
3.444 74<il5-0l 

3 .4*1 064*5-01 
*.12006?*5-M 
*.20616775-01 

4.77141T9F-4! 

6.141*1305-41 

|,7ftM*4F«<*l 

4*2*71 1015-M 
4 . §4 *46765 -0 1 
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7.4! 2441 05-01 
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*.6*44*725-M 
1 .•4534355*00 
1.051075*5*04 

1. 042344 ?5 *04 
1.1314*75€*90 
1.17*»*?75*00 
1.213*41*5*00 
1 .246101 05*44 

1. 247741147 *44 
l.)*176*7f*44 
1 . 34457444 *44 
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I, 47444444 ♦<»» 

1 .54443?4f*44 
1 .41205*25 *4# 
1 • 454765*5 *•* 

1.74447474*1*4 
1.74344447*44 
1.44|?T74r*^ 
i ,i«4*?*3r**» 
1.437*4^4C*f4 
1.4W1»47*44 
7.33444*1 F *44 
l,MMM?t*44 
7 # |1*4444F*44 
7,14**?47F*74 
7.21441 T7f *44 
7.23443374*44 
7.277344**. 4* 

7. 7734 4*14. p4 
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7 . 2 * 3 * 4417*44 
7.2*13334 4*4# 
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3. *30*1 775-41 
3.0V75674C-01 
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4.44771 902— 4 | 
9.4«43»**C-«1 
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3. U96299C *99 
3.9932199C*99 
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3. 93933995*99 
3.09977995*99 
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3,14413995*99 
3.22343995*99 
3, 25434995*99 
3,20313995*99 
3,30241995*99 
3.31515995*99 
3.31943095*99 
3.31419995*99 
3.33341995*99 
3.34324995*09 
3.35964095*99 
3.3749*095*99 
?. 399*9905*99 
3. *2733495*09 
3.4670*995*09 
3,404*4095*99 
3,»nMl0r*M 
3.52771995*09 
3.5*26*995 *99 
1. 66742004*09 


3 ,5*3490 95 *09 
3.549??O0?*99 
3.5*74500? *00 
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,64*645*0* *5 (.70 • 
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0* 
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ft • -.fftl 331*44 ** • 

XL • .433T4C-41 TE • .*9455 C*«4 

¥ a ,RAft04E*ftf Tip* a ,c9*.TlC*ft4 

Of /OX » -.6T451E*0* TAW ■ .T>9A7*E*ft4 

pwp* a ,1*?6 gE Oft 

UF * .?43AVE*03 

p£ * .SftT2A£*05 

CONTOUR PROPERTIES FLO* ft»ftPEt»TIE$ 

ft a ,,,4*ftf9f*flft * • .T|25TE*01 

XL • ,»?4?S£*OP Tr « •S064ftC*04 

V a # «tssft£*pr 7* * .«W*0¥E*ftA 

Of/O* ■ *.AT45l£*ftf. Tft* « .S966*»E*ft4 

DH/OS « *1 ftftT|E*Oft 

t»£ « .1ftft70£*i»3 

PC « .*A60o€*fA 

contour paoPF*»nes flow properties 

ft • •«RC4?6C*Cf « • ,«*?ftft4E«*91 

XL * ,ia*l2F*ftf If • .*9*40E*ft4 

¥ a .4PftftOC*«ft 7* a .«*9007£ 04 

Of /OX » ••^ME*!'? T4B b .*9ftATE*«4 
Ox/DX s .?319*E*04 

l*£ s .3U44CM3 

PE * .Sft435E*0f 

CONTOUR PROPERTIES flow POOPERTIE^ 

1 a ••S*¥T2E*ft4 X • , 9?0?*£-«l 

ftt • ,>S349E«r* Tf • .S9*??C*f4 

¥ a ,4444fE*4* ?*» * ,«.0fc*4C*ft4 

Of /Ox • -*Pf4ftl€*ftft T*a a ,A94ft*%'*ft4 
pX/OX a .?«3MC*"ft 

"E • .4*3A|C*33 

PC • .:.*>SS4C*0$ 

CONTOUR PaoRfollf* FLOW PoCPtaf JES 

ft • • f 44||ir»^ X a ,Jl?9?f*ftft 

ftl, a .¥1AH7 £.m n * .*440*1*44 

¥ a ,aaaiH*** 7* * ,«*fAA*F.ft4 

Of /OX • •*Pf4MC**« 7ia • *04 

■ *'*79S4E*fta 

t • *4741 JC*01 

Pt a f H4)4(i«4 

CONTOUR »3fPf»TI£5 *LO* P»;»7f«7lCS 

x , * a • iHl’f *PP 

a fft a ,’«F*£*'** 

V v ’* * * 

Of /Oft • •#474«c!*P4 7*« • .«.9*Mt*»4 

m ,4«04f*f*» 
i*f * ,4 MTa£*S3 

•f ■ ,44?4ftC»«: 

C04TM 'i.4* PP'" OPTICS 

a • - 1 . 447F*!** 

«i • ,44fMl*4»* »f * ,M*Np¥|aP« 

V a ,¥¥*4||.M • *..**/?,/% 

Of /0« a •»4»«4|.E**<* **. • *%*>♦ *fa 

* : 


DELTA * *33T33£-ft2 HO • 9. 

PEL TO • .RI200C-03 OH • •# 

NTITP a .4fl0*£-e3 SOHO • 0. 

fHTTA a •327f?E«l3 rORCE • ,214«4C*(I 

PHf a * 734ft 2£«04 LAT.F a ft. 

h a *i 2 rooc»oi * 


ft oo ho A pt later heat transfer 

DELTA a t 4tt4ftC-02 H« • ft, 

OCLTB a ,T0R27C~03 OH • ft* 

DfLT* a •SAf7fC»ft3 SOHO • ft* 

THE7A a ,427{|El3 FORCE * *4|?OoC*ftO 
PHf a *S3ftl«Cf ♦ LAT.F a 9 . 

H a .tNMCHt 


ftOUHOANV LAVER HEAT TRANSFER 

OELTA a .474f3£-02 Hi a « t 

ftcLTB a ,A 2033 E -«3 OH • ft* 

PEL TP a «S0f9SF.f3 SOHO a «, 

T^T* a ,00S0t£-43 FORCE • .4310S£*00 

PHI m •)73ftt€«04 L‘T.T • ft* 

H a .URflEMl 


ROUNDARY LAVER HEAT TRANSFER 

DELTA a »4«t7t£-ft2 hO a ft* 

DELTA a ,377«*E«ft3 OH a ft, 
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“« • 


y a 


•3i??4f »P) 

IP a 

•4iZ4if.it 

«ll»* a 

•i«9m»#4 




• S7T9K*4S 

OTL r 

a 

•itl|*C«tt OCL 1 • .lYittt'Afk 





CONTfttJ* POPPfOTlFS 

FLOY POOPCOTU* 

•OUHOFOY LAYfO Mf AT TH^VtF 

IMCONAL 

ttrciOALi 

cocfficicyti 

X • 

•rt*iM«t i 

Y • 

• IXTlff.M 

pfl TA 

a 

•49SA4f.it Hi a f, 

ZITA a 

•?4t*tf«M 

CF a 

•t72t4C.it 

»L • 

,1»*|tf.*l 

Tf • 

.?rom<A* 

Oft IV 

a 

•)4474f.A9 OY • P. 

IA a 

• liftK.if 

CM * 

•• 

Y • 


TV ■ 

.93)341 **4 

of IT* 

a 

. lt1A7f.il $0Y0 a p. 

!1 a 

• Ulitf.il 

OCT* a 

•19IT4f«AA 

OY/OX • 

• YXpICW 

T*V • 

.43)W*« 

THTfA 

a 

• IflSK-it FOOCf ■ .|70)lf*4S 

IA a 

•IlYttC.iS 

*f*L • 

•!«?1/C*aV 


OY/0* ■ 


PMf a 


* aAT?1C*p 7 LXT.F a p. 

17 a 

• Ultlf.pi 

4 1 a* • 

•rutf-«i 



l*f • 

.PH91*f9 

Y a 


.1i74Tf.il 

IIP a 

•AitSiC.it 

010* a 

• K*44 



Ff ■ 


Oft F 

a 

•if jTTf *4f OCL I • .1 7|f7f*it 





COtfTM* OOfPfoTtfl 

FLM* FYOFfOTtfft 

•OgHOAHY LAYfO MC«T TlAXfffl 

mi 

i 

m 

lYTCiOALt 

COCFFIClfYTI 

X • 


Y V 

.39TIX*FI 

OfLTA 

a 

.494iTf.il NO a A. 

ZCTA • 

•ZAtAlf.ii 

CF • 

• 7T77K.it 

«L « 


TF a 

,TTM4f««» 

OFi TO 

a 

• 1Mi4f.i9 tHI a A. 

IA a 

• 1i9tK*t9 

CM a 

4* 

V • 

«*•)**€ ««* 

TV • 


OF| fa 

a 

• UIXff-M 9«Y0 a p. 

If a 

•tiiiif.il 

•CT* a 

• |S|71f *A4 

OY/OX • 


TAV ■ 

•MSS9f«t4 

THFIA 

a 

.1734 V.P? FgkCC a .)H|X»lS 

IA a 

• t37tH.il 

K«l ♦ 

• UTlAf**? 



OY/Oi ■ 

i499IX*4« 

VMf a 


•44T44C-AT LXT.F a *. 

17 a 

• 1iStiC*M 

Ofa* a 

•?)4tar.it 



a 

.|F|!XM« 

y a 


.1iF90C*«l 

IIP ■ 

•A4t|At*it 

OCO* • 

•ttlfer*** 



•f • 

.)WX*0 

OfL F 

a 

.it094f.it OCL t • *tTAiXC«il 
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C outrun paoPrPTl^S flow Pt>rpE n U c S pounoary La*e« 


x . ,??msE*M 

XL ■ 

V a ,091*fE*''0 

OV/O* ■ ,?S**3f »*)f 


m • . 3 O? 3 F *01 nfl.T* ■ • 45^1 re-c 1 

T r - , >'P<T,t'.04 PF.i TB ■ ,1B<Jft6E-0S 

*W a l«-V??35F* ft 4 m T* ■ .123B?F-ftl 

-. tl0 w ,*;?T35r THfT* ■ .3?365f-ft2 

•j M/0> a .443*9600 P M t ■ .4*7Sftr.C7 

u g a .l^'OE^AS H a »3fi?5ftE*01 

Pf « ,375P9E*ft3 OH F ■ .Ac4|4E*C2 


HO ■ 

o* • 

SUHO ■ 

force <* 

lat.f • 


fOVTOUP P 00 PF 3 T 1 FS *l.OW PPCPrpTlLS P 0 ONP 4 RY LAYER 


X a ,??»76E** 1 

XI a ,X?*36F*«I 
Y a # 00]r7f*rP 
?Y/f)x » ,?®333C*^r 


* a .3S721E01 

TF * ,?7«83F*ft4 

TW a ,5?33^ *ftA 

TAW a ,S? 33 bE«Q 4 
pw/C.x a ,srt03F*<>r 
If a ,1ft '7 ft f *ftS 

pf a .JT'jSRUoS 


ofi.TA • .4S4UE-01 

pfl.Tg a ,jpWft7E»05 
pF» T* a .123*3F-M 
THCTA a .32366E-0Z 
PH! ■ .4B751F-07 

h a .3B?59E*01 

DFi F a .«J7S?E*02 


MO ■ 

0« ti 


CONTOUR PROPERTIES -*LON PROPERTIES ROUNP4RY l A *E® 


X a .fFftlOE**! 
*L a 0’437C««! 
Y a ,ORlhRf*ftft 
OY/D» * .l)Sft7E«ftO 


m a 

T r * 

N a 

TAW • 

r«/n* 
•■C * 
PF a 


,357,’af Ol 
.?7Rtf3f *1-4 
,S?3?5t ♦ '* 

.SP335F04 
.**70?E OP 
.IftXlftf •** 

,37S8Tf,«3 


PFl T 8 a .4S411F.01 
f>Fi Trl = .JP’fTF-CS 
or, r* - .!??A3F-M 

THt Tft a . T,*34fc s -ft? 
pH T « ,4fiT*lF»07 

H » ,J8?S9F*"1 

j£L F « .81 YM>5r* 02 


TRANSFER 

internal 

INTEGRALS 

COEFFICIENTS 

0* 

ft. 

0* 

• 17M9E*ft3 

ZETA ■ 

14 a 

15 a 

16 * 

17 a 

U» a 

•?6241E»00 
.10523E-05 
•loiRir-oi 
.13722F-05 
• 1039ACO0 
.6G215E-62 

CF • 

CH a 
RETH a 
RExL ■ 
REPH • 
RED* ■ 

•27222E-02 

0. 

,|5573E*04 

,15703E*07 

•23457E.01 

,5fS»l£*04 

.17?0SE*01 





TRANSFER 

internal 

integrals 

COEFFICIENTS 

0. 

ft. 

ft. 

» .l7R19E*ft3 

■ 0. 

ZETA • 

14 « 

15 a 

16 ■ 

17 a 
l\P a 

.26241EO© 

,lft523E-ft5 

.1M81E-01 

.13T21EOS 

,lfl390F*00 

•6S2l5E*02 

CF a 

CH a 
RETH a 
RE XL • 
RCRH ■ 
REO* ■ 

.27222E-02 

0. 

,|5573£*04 

,IS703€*07 

.23457E-0I 

,S9StIE«04 

a ,l7o67Eaftl 





Transfer 

INTERNAL 

integrals 

coefficients 

ft. 
ft • 

o • 

* .!7»l9E*ft3 

a ft. 

ZEtA a 

14 a 

15 a 

16 a 

IT a 

UR a 

,?624lE*ftO CF a 
•1OS23E-05 CH a 
,|«ltlE-ft! RETH a 
.IJTJIE-^ PUL * 
,103901*00 HEPH a 
•68214E-02 RED* • 

.272Z2E-0* 

0* 

,155T3C*ft4 

.1S703E*0T 

,234S*E*0i 

,595»1E*04 


hG 

QW 


PEL 1 a .»69UE*0l 
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throat radius corrected for displacement thickness • •itsisaotc^oo 


T*RLE or NORMALIZED CONTOUR POINTS CORRECTED FOR DISPLACEMENT THICKNESS 


OATA POINT 

X 

V 

1 

-•4606556.37*61 

.32900§50E*01 

2 

••3T1Z6336E*(U 

.30§9bl34E»Ol 

3 

-.3418l9?or*Ol 

.269034521*01 

4 

••312351557*01 

.2a9U480E*01 

5 

-.?820A83OF*O* 

, 249235737*01 

6 

-,?53375*9E*o1 

. 22936oloE*Ot 

7 

••?238780BE*C1 

,20949276-*0l 

8 

-.194376847*01 

• 189631 3SE*0l 

9 

..16487376E4M 

,169771 7oE*0l 

1« 

— 13S37AZ6F-0I 

,U991343E*0! 

M 

-.1458SR3SE-PI 

• 1 302 161 3E*01 

12 

-.76347025E*00 

.115US02E.81 

13 

-.468403887*60 

.10S5579PE*01 

>4 

-.1 733S366E* 00 

.I0075o57t*0l 

IS 

• 1 21 7?81 IF *08 

. 1 00351 2BE*0l 

16 

,48S27223€*00 

. 1 65905037*01 

17 

,5508?559E*00 

, 1 n764C93!«01 

18 

•All 321597.00 

,10947H8E*OI 

19 

.0566*8687*66 

. 1 1 0 7*7U5E*0l 

ZO 

.74271*677*00 

.11415S13E-01 

PI 

,64*154667*00 

.116901257*01 

?? 

.8755*7717*00 

.11969947E«01 

P3 

.*41692287*00 

.122483817*01 

?4 

*100782597*01 

.1?S26921E*41 

PS 

•1074p23R r *4l 

t 12BP56S7C*8l 

P4 

,11483120F*<?1 

. 1 3065n 79C«0l 

P7 

,|2f»67Tl9F*0l 

.133fc498«C*0l 

Pi 

,12735?J27'M 

.1 36462687*81 

P9 

,13405458*'M 

.13*284217.01 

1* 

.1448400SF*61 

.I4?l274?7*0l 

31 

.14759846F-M 

,U499f 24F.01 

*2 

,15444*487*01 

, uTsrt.'-.rr.oi 

.13 

.I0l39976f*6l 

,’Sr7P<548F*fll 

“*4 

. |68435)*F *M 

.;5374 ? 7kf.0l 

*5 

.1 F559182MM 

. 14M433H.&1 

•*6 

.1R289lS0F*nl 

.14S>O*l4i.01 

77 

.19e3K?lor.oi 

• 16?93ll> r#oi 

38 

,19i6PJfc8r*M 

, 1 461 4*9rC *0 1 

39 

.?0S9?62M *M 

.149469677.01 

44 

.punreoc-ai 

,i729ojorr.oi 

41 

•???4?33of*oi 

,19844 l 74E.fll 

4? 

.P3151SP1F *«1 

,18‘?2667E.0l 

43 

,744»SPl<!»*f11 

.16417^007.01 

44 

,?043??66r»M 

,194667561 .01 

4S 

.P91T.S520M 

, 2" 641 R 1 PC*M 

46 

,374.T??«.**6l 

.217-27361 .01 

47 

,3Sl043*lf*0l 

.?304«?06f«tl 

48 

*3R?40l2«F*-»i 

, 242900647*81 

49 

.4l626744fM 

,266697177*01 

«A 

.439?lf48C*nl 

,266947517*61 

Rl 

.463944141*41 

.276)2400:. 61 

4? 

.4M79«S9F*M 

.2*1 774651*01 

S3 

*418776467*01 

.36421R59F.6I 

S4 

,44p4440Sr*0l 

.31*624137*61 

SS 

*46746480r*M 

.3310*1441*61 

Si 


.340534787*41 


57 

.6 ne*?97e*m 

.3*88348?E*01 

as 

,6360S33?F*ftl 

•35619973E*01 

69 

,65«70779F*fll 

•36344if CE*01 

an 

.f8^l834C*01 

•370936W4E*01 

61 

,76f0fl7?2f«flt 

, 378989; ?E*01 

•>? 

»7?94l236F*61 

•387C?644E*01 

*3 

• 7519l91l>E«0l 

,395S9lu8C*6! 

A* 

.77R66t82F«0j 

•4643b874E*01 

6^ 

.46*3l565E*0l 

.4l333f77C*0l 

66 

.»3o?2fU8F«ei 

,4222726iE*Ol 

67 

,65639l95F*n 

,43US60?E*01 

68 

.8B?06653F*Al 

c 43956f09E«Ol 

69 

•9A64A3B6E*0i 

.447l9676E*e 

7<* 

• 93l86483f«M 

,4S*tt094SE*0! 

71 

.98?*9Ri9E*Ol 

.462l7«7rE«0 

7? 

•981686088*01 

.4«94?77oE«o: 

73 

•16P91833F*0? 

•4766ni?9E*0 

76 

•1«346717E*62 

.463*81 ,7E*0 

7$ 

• 1 Of 1 76068*02 

c 49C7636lC*0 

76 

.J699p398P*8J 

,*901lffSE*f 1 

77 

• 111M316E**? 

# %oS37n»*C*0 

78 

.l!*41 1 77F *02 

• 5nc*87?C.O 

79 

.1 7278668*6* 

.5?066ir*f«8 

66 

.119942988*62 

,a?793*0PE*8 

61 

,1??776*68*6? 

,53b35736E*0 

6? 

•i?aa7408F*62 

,4*2S9(kO>E*0 

“3 

. l?69683?r*fl? 

, 5%M *4 1 6C*0 

*4 

. l?769im*0? 

,547921ME*0 

8*. 

.1?86*614E*6? 

.S*877*2l£,4 

66 

.;?R7?2B28*62 

.549**3lH*0 

67 


.5*94fl4?F,0 

AM 

, ' ?Ml5«78 •*'* 

,6494?90fE*0 

66 

. ] ’*0!9<4F*f)? 

( 549*9*2 1 £ *0 

9* 

.i?r*fo*ee*e? 

.549 S 5*29C*0 
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SUMMARY OUTPUT prom TRL NOOULr 


A/A* 

OELT* ISPfTBt) 

2.AO009 

.2AAT0 

A.AftOP* 

*4tAA) 

A.ftflOftA 


1A.AAOO0 

•7««A» 

?A»AAQ0A 

».3«432 

3r»***M 

1 .66419 

3' *93857 

! *691*9 
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SUMMARY OF RESUL'S from SPP COMPUTER PROGRAM version S.I DATED 31 JAN7S 
CASE TlTir EMENDED Of L TA VALUATION CASE 
•••• VACUUM PERFORMANCE •**• 


otf 

LOSS 

eta 

ETA 

ETA 

eta 

MET MOO 


MECHANISM 

CALC 

FMRIR 

INPUT 

SELECTED 

SELECTED 


FTNFTIC 

I.aooa 

.9030 

l.POOP 

.9939 

E 


PfWpRM 

,9*?R 

,9H5l 

l.eopo 

. 9478 

C 

> 

COMM EFF 

,«S*9 

,9T JS 

\ .0000 

,9P*3 

C 

J 

Si'RMfRG 

l.ftOFP 

.9921 

l.OOflD 

.9921 

E 

A 

?*o 

l.OOGP 

.97*? 

l.nooc 

1.0006 

C 

P100UC1 *r FT AS 

.9329 

*9?|? 

l.ftooo 

.9198 



TURSUIE^T HOUNOARV L» v £« LO?S DECREMENTS 


mftmOO OELTA ISPaBF'<EC/CRMi 


FAlCOCATFO 

FMRfRfCAl 

fMRUt 

value sctFcrro 


1.68*4 
7 • T3*T 
n.AOCr 

i.*8%* MfnPO SELECTEO c 


DELIVERED vacuum SPECIFIC I^PUlSf alF-SEC/L0M| 


MfYMOO 


ISAlVACin 


* alculatfo 

*-RfRICAl 

TN»UT 

VAI.UC 'Fl'fTtO 


?«7.59?« 

2AJ.V'A6 

oir.otta 


ISR OFCRfMfNf {'uA *0 F*m*us?JL». 10 AM*! fN! PRESSURE 


• VF. 4 «R|FmT P**F„$;jPf IPSfA i .*‘>7’ 

*VF, CMAMMR RPFSSoRF (P$t A! »3l *• ' }* 

•VF, EAPINS|ON RATIO ?* 

PTSfMARRF fOEFF. !.P\M 

rriTA ISP »IPF«<FCAR»»‘ 


TSR OrtlvMJFO TO TMF Avra*f.r »M|fT parSsi'Rf OF .ff9?j* Ml* tS • 2P?*?H8 LtP-ttC/tM 


5-72 


thrust 


T!*c HISTORY 


r 

</M 

F-VAC 

R6**C 

0.00 

30«p9 

13131. 

35. 

2,66 

30 .nr 

13955. 


*,86 

29.91 

13776. 


6 ,0* 

29.7° 

1 3567 , 


6,88 

29,66 

13*39. 


18.88 

29.5* 

13895. 


12.88 

29.43 

176*7, 


1*,00 

29.31 

17960. 


16,88 

29.1* 

13*75. 


16,08 

29.16 

13*99. 


?*,08 

28.96 

13691. 


?2,o* 

28.89 

1 *100. 


?*, 88 

28. *74 

1*559. 


26.88 

28.67 

1*664. 


76,88 

28. 5r 

15771. 


1 *.0* 

28.3« 

1*650. 


32.00 

28.77 

1*696, 


3*.o* 

*••16 

16660, 

*2. 

36,00 

28.8* 

16198, 


36.00 

27.93 

16259. 


48,60 

27.63 

163*6. 


*8,88 

27,7? 

16158, 


44,88 

27,67 

1*869. 


46,08 

27,68 

5861, 


*6,60 

77.60 

1391. 


*8,66 

27,60 

161, 


52.06 

27.60 

5*. 

*3. 


m, fS6» 76?. 6* im-sfc/i#* 

ItlfUMtt rfHRL t*t CUTION TfHf TH*S Cill • 


F-ort 

1 -VOC 

ta*t 

1-071 

13296. 

284,49 

. 7 * 

264.25 

13519. 

2 **, 93 

• 7* 

764,10 

137*?. 

264.86 

.7* 

264,12 

13538, 

264.73 

.77 

263,96 

13**1. 

264.59 

.79 

?63.79 

13057. 

264.42 

.63 

283.59 

1280*. 

284.27 

.66 

?63.*1 

1286*. 

264,15 

.67 

263.26 

13*35, 

284,11 

.65 

263.26 

13*5*. 

264,00 

.66 

263.1* 

136*9. 

263.90 

.65 

263,05 

1*05*. 

283,6* 

.*3 

263,00 

1*51*. 

2*3,77 

.61 

762.96 

1**27, 

263,69 

.60 

262.69 

15224. 

263,62 

.76 

2*7,6* 

156*8, 

263.5* 

.76 

262.75 

15*64, 

783**5 

.75 

262.69 

1601 *. 

283,35 

.75 

262.60 

16156, 

7*3.2* 

.75 

262.5* 

16216. 

283.1* 

.74 

262.39 

16387, 

2*3.03 

.74 

262.24 

16115. 

262.41 

.75 

262,16 

1**26. 

202.67 

.6? 

261.65 

5*36. 

278.41 

2.37 

276.5* 

1346. 

252.6? 

7.6* 

2*5.8? 

117, 

??u 91 

59.46 

16l.*5 

1*. 

2 ?6,*o 

161.93 

*3.*6 


* mi . **oot« *s.?2 utR/sec 


i.ioo iHiu) miimtt 


MOOT 

-oot.int 

*6.775 

0*08 

♦7.57? 

94,35 

*6.369 

190.29 

*7,6*7 

266.30 

*7.222 

381.17 

*6,0*1 

474.44 

*5.19* 

565.67 

*5.396 

656*26 

♦7.431 

749,09 

*7.530 

844,05 

*8.273 

939.81 

*9.677 

I03?.?l 

51.30? 

1138.69 

57.396 

12*2.39 

53.6*3 

13*6.63 

55.196 

1*57.67 

56.063 

1568,98 

56.661 

1681,71 

57.189 

1795,58 

57,423 

1910,20 

57.753 

2*25,3? 

57.111 

21 *0.2* 

57,602 

22*9.95 

16.219 

2320.77 

5.502 

23**.*9 

.7?7 

2350.7? 

• 237 

2351.69 


r*IMT 
0*0 
26815 ,* 
5*0?6.9 
813*9.1 
108 * 80.1 
13*738. 3 
16*603,* 
18627?.! 
21266T.7 
239460,$ 

266567.6 
29*275.6 
3228S2.4 
352192, t 
3622*3,6 
*P*81,0 
***5*1.* 
*76615,5 
5*8589,3 
5* *960, s' 
573*79.6 
605897.0 
636 « 38,5 
6567*2,4 
663*66,6 
664554, * 

664681.6 
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Sample Case Number 2 
Input Listing 


TfTLf e*Te*o?o on.** VALIDATION f k%r 

urn M 

*<114(1 maaumi. 

?.*Ll«** 

44|to fl)a?#4 t 6 t lft'2..«30« 

g>Aj9at, t 

TwC«1*, 

Mfrun* 

TtiFTA I *3* • 

Q'Y*>*14U 

?<!?)• 1.?M4A?.1.t«’.5.3.*«4 c 4* 

9«(>)4l # 37|?m.1.4M99*1.AA94AA.?.3Y*Q4l**. 7374156* 
1,?4«.3»4933 # 4 # 3«M4*4.44?63 f A.553246UT* 

S.MI’ • 

TOv»rcr^R7 

•*04.| 

®AM4t n )4 # 4A. # CW**09«* # l..1. 

TTO Jill ■«)•« I*. • **. .44,. I ftp., 
sit iof«4* 


99f^ L f,^ 

tno'.a 

DM41* MCOoU 
WM 

b*#CT*mts 

H i. *4, Cl 1.1 o 4,f. 

41 1.4 

C *.J10 * 10.404 A H.ftoo 


?e 

• 7ft»»O0 # 

5 

«<.i* 

'4 

ft. ft 

S 


14 

n*o«. 

S 

(H.l> 


N»«r.!$ts 

r*J4T4».3o» 

»«I4«T. 

4«t>T. 
g, osMir. 
c***74T. 

^ • 4 f 1 1 • ^4.045«Mn)4|P«. 
«r*i' 
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Appendix A - Control Card Set Ups 


Control card set ups for two typical cases are shown for both CDC 6000 
(SC0PE 3.3) and Univac 1100 (EXEC 8) seires computers and operating systems. 
In both cases absolute versions of the SPP code and the thermodynamic data files 
are assumed to be on permanent mass storage data files. 

The first case considered is when one or more of the computational modules 
are executed and the linkage data is written (punched) out for subsequent runs. 

The second case considered is when linkage data from the first case is 
used in executing subsequent modules for the same motor. 

The bracket symbol ( | ) is used to explain or comment on the function of 
the control cards. Since the first, and most likely the second, control card of 
each set varies from installation to Installation even for the same computer and 
operating system. The accounting information control card is left off from each 
set and the job specification control card (normally the second control card in 
the set and the first shown here) should be checked for consistancy at each in- 
stallation. 

CASE Is CDC 6000 Series Computers 

J0B, CM 10000 , CL126000 , 1700*, 101000. 

ATTACH, SPP, SPPABS, J get SPP absolute 

ATTACK, TAPE 25, SPPTHERMG j equate thcrn.o c'ata to logical unit 25 
RFL, 126000. j request core to execute SPP program 
SPP. | execute the SPP program 
7/8/9 

program data 
6/ 7/8/9 end of file 

CASE I: Univac 1100 Series Computer. 

@ASG,AX SPP. j assign program file 

&ASG.AX SPPTHF.RM0.* | ussina thermo date 

@USE 25. , SPPTHERM0. * | equate logical unit 25 to thermo data 

@ASG f U LINK DATA. ,F J assign a file for the linkage data 

@USE 8. , L1MCDATA. j equate logical unit 8 to the linkage Jata file 


A-l 


@XQT SPP.SPP j execute the absolute element SPP from the SPP file 
program data 
@FIN 

@@ | remote Job entry end of file 

CASE II: CDC 6000 Series Computers 

J0B , CM10000 , CL126000 , T700* , 101000 . 

ATTACH, SPP, FPPABS. 

ATTACH, TAPE;’5 , SPPTHERM0. 

C0PYBR,INPUT,TAPE3. | copies previously generated linkage data to TAPE3. 
REWIND (TAPE3) 

RFL, 126000. 

SPP. 

7/8/9 

linkage data from Case I 
7/8/9 

program data 
6/7/8/9 

CASE II: Univac 1 100 Series Computers 


@ASG,AX 

SPP. 

@ASG,AX 

SPPTHERM0, + 

@USE 

25. ,SPPTHERM0. + 

@ASG,AX 

LINKDATA.** 

@USE 

3., LINKDATA.** 

@ASG, U 

LINKDATA2. ,F assign and equate further linkage 

@USE 

P . , LINKDATA2 . data to unlt 8 

@XQT 

SPP.SPP 


program data 
@FIN 
@@ 


NOTES: 

* CPU time limit depends on modules to be executed. 

+ Only needed if the 0DE and/or 0DK modules are to be executed. 

** If this data has been punched, the following control cards can be used: 

@ASG,T 3..F4 
©DATA, I 3. 

linkage data from Case I 
@END 
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